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Section  I 


INTRODUCTION 


There  are  three  programs  in  the  set  of  computer  programs  called  GOBACK. 
The  purpose  of  these  programs  is  to  accept  the  coefficients  of  a  curve 
or  surface  equation  as  produced  by  the  linear  programming  formulation, 
and  from  these  to  calculate  the  basic  loft  information  to  produce  the 
ship.  The  use  of  each  of  the  three  programs  is  given  below. 

GOBACK  1  -  Calculates  information  from  a  surface  equation 
which  contains  no  profile  requirements 

GOBACK  2  -  Calculates  information  from  a  surface  equation 
containing  profile  requirements 

GOBACK  3  -  Calculates  information  from  an  equation  describing 
a  curve,  and  can  plot  the  curve  at  the  same  time 

Sections  II  and  III  of  this  Appendix  describe  the  application  of  these 
programs.  Section  II  describes  the  two-dimensional  case  handled  by 
GOBACK  3,  and  Section  III  extends  this  to  the  three-dimensional  case 
handled  by  GOBACK  1  and  GOBACK  2  . 
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Section  II 


TWO  DIMENSIONAL  CASE 


A  typical  curve  which  may  require  analysis  is  shown  in  Fig.  D-l  The 
equations  that  may  be  used  to  describe  the  curve  are  also  shown. 


#  See  text  ("Mathematical  Ship  Lofting  -  Part  1.  Theory," 
Technical  Report  1.0. 0-1)  for  explanation  of  this 
notation. 

Fig.  D-l 
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Lino  (1)  describe  the  order  end  renge  of  Influence  of  the  coefficlente 
for  e  single  spllned  equation.  Lines  (2)  provide  similar  information 
for  e  double-splinad  equation.  Mote  that  the  offset  at  xq  equals  A 
and  the  equation  goes  exactly  through  that  offset.  OOBACX  3  is  capable 
of  providing  offsets,  first  and  second  derivatives  from  these  equations 


Section  III 


SURFACE  EQUATIONS 


G03ACK  1  performs  calculations  from  a  three-dimensional  surface  equation 
such  as  that  shown  in  Fig.  D-2 .  This  surface  equation  is  an  extension 
of  the  two-dimensional  case.  The  surface  equation  is  always  double 
splined  in  the  x  direction  (lengthwise  on  the  ship).  It  may  be  either 
single  or  double  splined  in  the  2  direction  (vertical  direction  on 
the  ship).  The  resulting  equation  for  either  representation  is  shown  for 
a  five-station  by  five -water line  surface  in  Fig.  D-2. 

The  program  has  the  ability  to  extract  the  following  information  from  the 
surface  equation: 

(1)  The  offset,  and  first  and  second  derivatives,  at  any  point 
on  the  surface 

(2)  The  offsets,  and  first  and  second  derivatives,  along  any 
waterline  or  station  at  any  given  increment  or  at  an  exact 
interval,  such  that  straight  lines  joining  the  offsets  will 
have  its  greatest  deviation  from  the  surface  equal  to  a 
given  tolerance. 

(3)  The  Theilheimer  equation,  or  the  standard  cubic  equations 
for  each  interval,  along  any  waterline  or  station. 

(4)  Heights  and  offsets  at  any  interval  along  the  intersection 
of  the  surface  with  any  given  spline  curve. 

(5)  Heights  and  offsets  at  any  increment  along  buttocks  and 
diagonal  planes. 
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Station  Number 

Surface  Equation  (Double  splined  in  x  and  z  direction) 
y  "  A00  +  A01X  +  A02x2  +  A03x3  +  A04(X-X2^ 

+  A10  +  A11X  +  A12x  +  A13x  +  A14(X“X2)#  |  z 

+  (^0  +  A21X  +  A22*2  +  A23x3  +  A24<x_x2^#j  ^ 

+  (a30  +  A31x  +  A32x2  +  A33x3  +  A34(x-x2)3j  z3 

+  |a40  +  A41x  +  A42x3  +  A43x3  +  A44(x-x2)3J  (gmg^3 


Surface  Equation  (Double  splined  in  x  ,  single  splined  in  z  direction) 
y  -  Aqq  +  AQ]x  +  AQ2x  +  AQ3x  +  Aq4(x-x2)^ 

+  A1q  +  Anx  +  A12x2  +  A13x3  +  A14(x-x2)3  J  z 

+  A20  +  A^  +  A22x2  +  A23x3  +  A24(*~x2)^j  z2 

+  A3q  +  A31x  +  A32x2  +  A33x3  +  A34(x-x2)3  j  z3 

+  A40  +  A41x  +  A42x2  +  A43x3  +  A44<*“*2)#  (z*zl>J 
2  3  3  3 

+  Ajq  +  Aj^x  +  Aj2x  +  A^3x  +  Aj4(x-x2)^  (z-Zj)^ 

+  A60  +  A61X  +  A62x  +  A63X  +  A64^*“x2^#  )  ^z"z3^tf 
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A.  WATERLINE  AND  STATION  EQUATIONS 


The  equations  for  waterlines  and  stations  may  be  extracted  from  the 
surface  equation  by  setting  x  (for  stations)  or  z  (for  waterlines) 
constant  and  summing  coef ficie-its .  Using  coefficients  from  equation 
(2),  Fig.  D-2,  the  equation  for  waterline  3  is  found  as  follows  by 
summing  vertically  with  z  ■  3  . 

*1  ■  A00  +  A10<3)  +  A20<3)2  +  A30<3)3  +  A40«-!)# 

B2  -  Aq1  +  An(3)  +  A21(3)2  +  A31(3)3  +  A41(3-2)3 

B3  -  AQ2  +  A12(3)  +  A22(3)2  +  A32(3)3  +  A42(3-2)3 

®3  -  Am  +  Au(3)  -HjjO)2  +  A33<3)3  +  A43(3-2)J 

»4  '  A04  +  A14<»  +  A24(3>2  +  A34(3>3  +  A«(3-2)# 

B5  -  A05  +  A15(3)  +  A25(3>2  +  A35(3)3  +  A45(3-2)3 

The  equation  for  the  waterline  is: 

y  -  B^  +  B2x  +  BjX2  +  B4x3  +  B^x-x^^ 

The  equation  for  a  station  is  found  by  summing  across  in  Equation  2, 
Fig.  D-2.  The  equation  for  Station  2  is  found  as  follows: 

Set  x  *  2  ,  then: 

C1  -  A00+A01<2)  +A02<2>2  +A03(2)3 
C2  -  A10  +  Au(2)  +  A12(2)2  +  A13(2)3 

C3  *  A20  +  A21(2)  +A22(2)2  +A23(2)3 

C4  -  A3q  +  A31(2)  +  A32(2)2  +  A33(2)3 

S  "  A40  +  \l(2>  +  A42(2>2  +  A43(2>3 

The  equation  for  the  station  is: 

Y  -  Cj_  +  C2z  +  C3z2  +  C4z3  +  C5(z-z2)J 
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B.  STANDARD  CUBIC  EQUATIONS  FOR  WATERLINES  AND  STATIONS 


Using  Che  above  two-dimensional  equations  we  can  now  develop  the 
standard  cubic  equations  between  each  of  the  points  where  an  additional 
coefficient  is  added.  On  single-spline  curves  there  is  one  in  every 
interval  between  offsets.  On  double-spline  curves  there  will  be  a 
standard  cubic  between  every  odd  numbered  point  (0,1,3,...). 

The  following  formulas  are  used  to  find  the  standard  cubic  equation  on 
a  waterline  curve  between  the  ith  and  ith+^  points  of  discontinuity. . 
Using  the  waterline  coefficients  developed  above: 

i 

°i  ■  bi-z  .  ■ 

i 

»2  '  B2+3IV*x? 

J-l 

i  ...  :  ' 

D3  B3  '  3X!  BJ+4  Xj 

.  J-l 

‘  i+4  . 

D4  “  B4  +  I  Bj 

J-5 

The  equation  is 

y  -  Dx  +  D2x  +  D3x2  +  D4x3  Xj  4  x  <  xj+1 


Using  the  example  of  Fig.  D-2  and  the  coefficients  of  Waterline  3,  the 
standard  cubic  coefficients  between  points  2  and  4  would  be: 
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Dx  -  B1  -  B5(l)“ 


D2  "  B2  +  3(B5)  ur 


D3  -  B3  -  3(B5)(1) 


D.  -  B.  +  B. 
4  4  5 


The  standard  cubic  coefficients  for  an  interval  on  a  station  are  found 
in  a  similar  manner.  Using  the  station  curve  coefficients  developed 
above,  the  coefficients  between  the  i  and  i^  +  1  points  of 
discontinuity  are: 


El  ‘  °1  '  L.  CJ«  *J 

]-l 

i 

I 

J-l 

i 

I  CJ* 

J-l 


E2  "  C2  +  3  L,  CJ+4  zj 


E3  "  C3  -  3 


The  equation  is : 


y  ■  E^  +  E2z  +  E3z2  +  E^z3 
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C.  BUTTOCKS 


The  following  procedure  is  used  to  obtain  the  heights  of  the  inter¬ 
section  of  buttock  planes  (y  ■  constant)  and  the  hull  surface. 

At  each  position  along  the  length  of  a  surface  where  a  height  is 
to  be  calculated,  the  two-dimensional  equation  for  the  frame  is 
obtained,  *  f(z)  .  The  equation  for  the  buttock  is  yfi  *  C  , 
where  C  is  the  distance  from  the  centerplane.  Subtracting  the  two 
equations  gives : 

6y  ■  yf  -  yb  -  f(z)  ■  c 

Trial  and  error  solution  until  |£y|<  £  ,  where  £,  is  a  small  number 

provides  the  correct  value  of  z  . 


D.  DIAGONAL  FLANKS 

It  is  sometimes  of  interest  to  calculate  the  intersections  of 
diagonal  planes  with  the  surface.  GOBACK  1  has  the  capability  of 
handling  diagonal  planes  which  pass  through  the  intersection  of  the 
centerline  plane  and  a  waterline  and  whose  intersection  with  the 
hull  surface  is  below  that  waterline  (Fig.  D-3) 


rD*  (**,-*) 

T 

Fig.  D-3 
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The  procedure  is  similar  to  that  used  in  finding  heights  of  the 
intersection  of  the  surface  and  buttock.  First,  the  equation  for 
the  frame  is  found,  then  the  equation  for  the  plane  is  subtracted 
from  it,  giving! 

&y  “  ys  -  yd  -  f(z)  -  tan  X  (zm-z) 

The  equation  is  solved  for  z  by  trial  and  error. 


E.  SEGMENTATION  OF  WATERLINE S  AND  FRAMES 

It  is  frequently  necessary,  primarily  for  purposes  of  numerical 
control,  to  be  able  to  approximate  waterline  and  station  curves 
with  straight  lines  such  that  the  deviation  between  curve  and  line 
is  no  greater  than  some  given  tolerance.  The  GOBACK  1  program  has 
the  ability  to  segment  a  station  or  waterline  curve  so  the  greatest 
deviation  will  always  be  exactly  the  specified  tolerance.  A 
development  of  the  method  used  Is  given  below. 

For  the  purposes  of  the  development,  a  waterline  curve  has  been 
selected.  The  method  applies  equally  well  to  a  station  curve. 

Figure  D-4  shows  the  geometry  of  the  problem. 
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Statement  of  the  problem: 


2  3 

Given  curve  y  -  f(x)  ■  A  +  Bx  +  Cx  +  Dx  ,  a  point  on  that 
curve  P^  and  a  tolerance  ;  find  a  point  on  the  curve  T!^  ,  with 
a  greater  x  coordinate,  such  that  a  line  joining  P^  and  P^  will 
have  its  greatest  deviation  from  the  curve  equal  to  i, . 


As  shown  in  Pig.  D-4,  the  greatest  deviation  from  the  chord  (L^) 
occurs  at  the  point  (P^)  where  the  tangent  to  the  curve  (1^)  is 
parallel  to  the  chord.  The  slope  of  the  tangent  at  this  point 
m^  ■  f'(x)  evaluated  at  x^  ,  which  is  equal  to  the  slope  of  the 
chord 


The  equations  for  these  lines  are : 

ll  ,  [»-«,] 

h  •  Y  *  Y3  “  “i  [*"*3] 

or 

L,  ,  Y  ■  Y^  +  m^  x-  iBjX^ 

1*2  ,  Y  ■  Yj  +  m^  x  -  mjXj 


Now  let 


Y1  ’  mlxl 


C2  "  Y3  "  mlX3 


now  the  length  of  the  line  joining  P^  and  P^ 

jL 


C1  -°2 


cos  (tan 


Therefore  Y^  -  Y^  -  ttijX^  +  m^x^ 


1  +  m^ 

% 

r 

equals 
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Substituting  the  equation  for  y3  ,  its  first  derivative  for 
and  rearranging  gives : 

£  1  +  (B  +  Cx3  +  Dx3)  -  (A  -  yx  +  Bx^)  +  (2Cx1)  x3 

J  2  3 

+  (3Dx1  -Ox*-  2Dx3 

x3  is  the  only  unknown  in  this  equation.  It  may  be  solved  by  trial 
and  error  to  yield  x3  . 

Substituting  x3  into  the  first  derivative  of  the  equation  of  the 
curve  yields  the  slope  at  that  point  m^  • 

Now  from  the  equation  for 

Y2  -  Yj.f  ij  (x2  -  xx)  '  .  .  . 

And  from  .the  equation  of  the  curve 

•  Y2.  A  +  Bx2  +  Cx22  +  Dx23 

Subtracting 

f  2  3 

0  -  .  Yx  +  mx  (x2  -  Xj)  A  +  Bx2  +  Cx2  +  DXj. 

This  equation  may  now  be  solved,  for  its  on ly  unknown  ,  x2  ,  by  trial 
and  error . 

Finally  substituting  x2  in  the  equation  for  the  curve  gives: 
y2  -  A  +  Bx2  +  Cx22  +  Dx23 

Note  that  we  have  assumed  in  this  derivation  that  f(x)  is  a  standard 
cubic  equation.  GOBACK  1  solves  for  the  standard  cubic  equation 
coefficients  in  each  interval  before  segmenting  the  curve. 
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F.  INTERSECTIONS  WITH  CURVES 


It  is  frequently  necessary  to  find  intersections  of  the  hull  surface 
equation  y  ■  f(x,z)  with  previously  defined  curves.  For  example, 
to  find  offsets  of  sight  edges  or  longitudinals. 

The  curves  are  usually  defined  by  a  faired  Theilheimer  equation 
z  ■  G(x)  representing  the  trace  of  the  heights  of  the  curve  on 
the  centerplane  of  the  ship.  To  find  the  intersection  of  the  curve 
with  a  given  frame,  GOBACK  first  solves  G(x)  at  the  given  x  of 
the  frame  for  the  height  z  ,  then  plugs  x  and  z  in  f(x,z)  to 
obtain  the  offset  y  .  The  curve  G(x)  must  have  the  same  origins 
and  be  scaled  the  same  as  f(x,z)  , 

Although  GOBACK  expects  a  form  for  z  -  G(x)  ,  such  as; 

■z  ■  A  +  Bx  +  Cx2  +  Dx^  -  (x-x^)^  +  . . .  , 

it  is  clearly  possible  to  use  degenerates  of  this  function,  such 
as  z  ■  Cx  by  making  the  other  coefficients  zero. 
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Section  IV 


SURFACE  EQUATION  WITH  END  CONDITIONS 


When  the  surface  which  has  been  faired  contains  an  end  profile,  the 
equation  requires  a  special  function  which  guarantees  that  the  surface 
will  assume  the  shape  of  the  pronle.  The  profile  function  has  the 
effect  of  squeezing  the  surface  down  to  the  required  profile.  This 
function  usually  is  effective  up  to  about  one  station  space  from  the 
profile. 


The  surface  equation  now  becomes 

y(x,z)  -  F(x,z)  *  T(x,z) 

where  F(x,z)  is  the  original  surface  equation  and  T(x,z)  is  the 
profile  function. 


The  function  equals: 

T(x,z) 


I  1 
L 


3 


P 


0  ±9L*± 

D 


# 


where  G(z)  is  a  two-dimensional  Thielheimer  equation  of  the  profile. 
G(z)  has  the  same  origin  and  is  scaled  the  same  as  the  surface  equation. 
D  is  the  fraction  of  the  station  spacing  over  which  the  function  is 
effective  and  P  determines  if  the  slope  at  the  profile  is  finite  or 
infinite.  If  P  has  a  value  1  the  slope  is  finite  ,  a  value  of  1/3 
would  give  infinite  slope. 


A.  EQUATIONS  FOR  WATERLINE S  AND  STATIONS 

The  equations  for  water  lines  and  stations  are  obtained  in  the  same 
manner  as  in  GOBACK  1  .  Offsets  are  obtained  from  these  equations 
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in  the  same  manner  as  GOBACK  1,  except  a  value  for  the  end  profile 
function  must  be  found  and  the  result  multiplied  by  it.  This  is 
done  by  GOBACK  2  .  First  and  second  derivatives  are  calculated  as 
shown  in  Appendix  B. 

B.  BUTTOCKS 

Because  the  GOBACK  2  program  is  required  to  calculate  data  in  the 
extreme  ends  of  the  ship,  a  somewhat  more  sophisticated  routine 
for  finding  buttocks  was  included  in  it. 

As  shown  in  Fig.  D-5  the  buttock  for  a  ship  can  be  single,  double, 
or  triple  valued  on  any  given  station.  Also,  different  buttocks  on 
the  same  ship  are  different  cases. 

If  the  hull  shape  is  such  that  the  case  shown  in  Fig.  D-5a  can  be 
guaranteed,  the  program  will  calculate  buttocks  as  shown  in 
Section  III. 

If  there  is  a  possibility  that  one  of  the  other  cases  is  present, 
the  following  procedure  is  used. 

(1)  The  equation  of  the  frame  is  calculated  from  the  surface 
equation. 

(2)  The  intersection  of  the  frame  with  the  baseline  is 
found  (y  )  . 

V 

(3)  The  coordinates  of  Point  P  (yp,Zp)  are  found  by  examining 
the  slope  of  the  curve  at  intervals  starting  at  the  base¬ 
line  until  the  slope  changes  from  negative  to  positive. 

The  Point  P  is  the  exact  point  on  the  curve  where  this 
takes  place.  If  no  point  P  is  found,  then  this  is  also 
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Hull 


Fig.  D-5 
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the  case  shown  in  Fig.  5 -a  . 


(4) 


The  values  y^  and  y^  and  compared  with  the  y  value 
of  the  buttock  plane  (y^)  . 


If: 


yq<yb<£  yp  the  case  is  that  of  buttock  1  Fig.  5-b 
and  only  one  intersection  will  be  found 

'  yq<yb  >yp  .the  case  is  that  of  buttock  2  Fig  5-b, 
and  three  intersections  must  be  checked 
for 

y^>yb<y^  the  case  is  that  of  buttock  1  of  Fig.  5-c 
and  no  intersections  are  present 

y^yb^>yp  the  case  is  shown  as  buttock  2  Fig.  5-c, 
and  two  intersections  may  be  present. 


(5)  The  coordinates  of  the  intersections  themselves  are  found 
by  finding  a  trial  and  error  solution  of  the  frame  equation 
minus  the  buttock  equation. 

The  first  frame  cut  for  a  given  buttock  line  is  found  by  finding 
the  intersection  of  the  buttock  with  the  uppermost  waterline  on 
the  surface.  From  this  point,.-  frame,  cuts  are  made  at  a  specified 
interval  along  the  ship,  and  the  intersections  of  each  frame  with 
the  buttock  is  found.' 
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Section  V 


GOBACK  1 


OPERATING  INSTRUCTIONS 

This  version  of  GOBACK  will  accept  coefficients  of  a  surface  equation 
such  as  shown  in  Fig.  D-2  which  does  not  have  end  profile  requirements. 
The  equation  must  be  single  splined  in  the  x  direction  and  may  be 
single  or  double  splined  in  the  z  direction. 

The  following  data  can  be  calculated  with  this  program: 

—  Offsets,  first  and  second  derivatives  along  a  waterline 
at  a  given  interval  or  at  an  interval  determined  by  a 
tolerance 

—  Offsets,  first  and  second  derivatives  along  a  station  at 
a  given  interval  or  at  an  interval  determined  by  a 
tolerance 

—  Standard  cubic  coefficients  for  each  section  of  a  waterline 

—  Standard  cubic  coefficients  for  each  section  of  a  station 

—  Offsets  at  a  given  interval  along  buttocks 

—  Offsets  at  a  given  interval  along  diagonal  planes 


Fortran  Input  Svmbols 

Symbol 

Definition 

XO 

The  actual  full  scale  coordinate  of  the  first 
station  of  the  surface 

zo 

The  actual  full  scale  Z  coordinate  of  the  first 
waterline  of  the  surface 

FINTX  - 

This  value  is  used  as  a  scale  factor  in  the  x 

direction.  If  the  surface  equation  is  scaled  so 

the  stations  are  one  unit  apart,  FINTX  equals  the 

actual  full  scale  station  interval  (x, -X  )  .  If 

X  0 

FINTZ  - 


M 


N 


SEGX  - 


SEGZ 


KS 


XS 


the  stations  of  the  surface  equation  are  less  than 
one  unit  apart,  FINTX  becomes  some  multiple  of  the 
station  spacing.  For  example,  if  the  stations  in 
the  equation  are  1/4  unit  apart,  FINTX  becomes  four 

fimes  the  actual  full-scale  station  spacing 

4<Vxo>] 


Same  definition  as  FINTX, 
(Z1  ‘  Z0) 


except  in  the 


Z  direction 


The  total  number  of  points  of  discontinuity  along 
a  station  (points  where  z  coefficients  are  added) 
including  the  points  at  the  first  and  last  water- 
lines.  For  example,  if  a  surface  containing  seven 
waterlines  were  single  splined,  M  would  equal  seven, 
since  there  are  third  derivative  discontinuities  at 
each  waterline.  If  the  seven -water line  surface  were 
double  splined,  M  would  equal  four,  since  there 
are  discontinuities  only  at  waterlines  0,2,4, 
and  6. 


The  total  number  of  points  of  discontinuity  along  a 
waterline  (points  where  x  coefficients  are  added). 
The  example  used  for  M  is  valid  except  using 
stations  instead  of  waterlines. 


Desired  Interval  between  consecutive  calculated 
offsets  along  water  lines,  buttocks,  and  diagonal 
planes . 


Desired  interval  between  consecutive  calculated 
offsets  along  stations. 


The  value  SEGX  can  be  changed  when  working  along 
a  curve  in  the  x  direction.  KS  is  the  number  of 
changes  to  be  made. 


The  x  distance  from  the  origin  of  the  surface  to 
a  point  where  SEGX  is  to  be  changed.  This  must 
coincide  with  a  station  or  a  specific  point  (see 
definition  of  K  ). 
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SEGXS  - 

A 

Z 

X 

IDENT  - 

K 

ZT 

XT 

NN 

ISP 


Value  to  which  SEGX  Is  to  be  changed  when  XS  is 
reached . 


Coefficients  of  the  surface  equation  produced  by 
linear  program.  There  are  (M+2)(N+2)  coefficients. 


z  coordinates  of  waterlines  where  coefficients  are 
added,  plus  that  of  the  last  waterline.  There  are 
(M-l)  since  the  value  for  the  first  waterline  isn't 
included. 


Same  as  Z  except  for  stations  (N-l) 


Code  which  tells  the  program  what  data  is  required 
from  the  surface  (See  Input  Format  description) 


The  program  makes  a  provision  for  obtaining  offsets 

at  specific  distances  along  a  curve  which  may  not 
coincide  with  the  intervals  specified  by  SEGX  or  SEGZ. 
These  are  called  "specific  points."  K  is  the  number 
of  points  asked  for  along  a  curve. 


The  z  coordinate  of  a  waterline  along  which 
data  is  to  be  calculated. 


The  x  coordinate  of  a  station  or  frame  along  which 
data  is  to  be  calculated. 


A  symbol  consisting  of  two  letters  used  to  identify 
specific  points  in  the  output.  The  same  symbol  must 
be  used  for  all  specific  points  on  any  one  curve 
(see  ISP  below) 


A  companion  symbol  to  NN  providing  two  additional 

characters  for  identification.  For  example,  if  the 
specific  point  on  a  station  corresponded  to  longi¬ 
tudinal  13,  NN  might  become  LG  and  ISP  become  13, 
as  the  specific  point  would  be  identified  LG13. 
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XE 


The  x  coordinate  of  specific  points  on  uaterlines, 
diagonals,  and  buttocks  (see  definition  of  K) . 


ZE  -  Same  as  XE  except  for  stations 

ANGUS  -  The  angle  in  Fig.  D-3)  between  a  diagonal 

plane  and  the  vertical  centerplane  of  the  ship. 

BUTTK  The  y  distance  from. the  centerplane  of  the  ship 

to  a.  buttock- along  which  heights  are  to  be  calculated 

SX  -  .  The  stair  ting  point  on  a  curve  when  segmenting  with 

a  tolerance  in  .the  x  direction 

SZ  -  The  starting 'coordinate  when  segmenting  a  station 

according  to  a  tolerance 

FX  -  The  finish  coordinate  corresponding  to  SX 

FZ  .  -  The  finish  coordinate  corresponding  to  SZ 

TOL  -  The  maximum  deviation  between  a  curve  and  the 

straight  line  segments  approximating  it  when  seg¬ 
menting  according  to  a  tolerance. 

Input  Data  Cards 

There  are  two  distinct  kinds  of  sets  of  data  used  for  input  to 
GOBACK  1  .  The  first  set  consists  entirely  of  data  describing  the 
surface  equation  to  the  program.  No  output  is  produced  from  this 
set.  The  second  set  of  data  consists  of  packets  of  data  cards. 

Each  packet  describes  some  information  required  from  the  surface 
and  gives  the  geometric  data  necessary  to  obtain  the  information 
from  the  surface  equation.  Only  one  of  the  first  sets  of  data  are 
entered  per  problem.  There  may  be  many  small  packets,  each  of 
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calls  for  offsets  along  a  waterline  or  station,  or  perhaps  standard 
cubic  coefficients  etc. 


There 

are  limitations  on  some 

of  the  input  parameters  of  the  program 

These 

limitations  follow. 

Variable 

Minimum 

Maximum 

M 

3 

20 

N 

3 

40 

K  (no.  of  specific 
points /curve) 

0 

20 

In  describing  the. various  data  cards,  the  actual  FORTRAN  format 
field  description  is  used  in  most  cases.  These  fields  come  consecu¬ 
tively  across  the  card  with  no  gaps  or  blank  columns  between.  The 
field. descriptions  are  the  FORTRAN  F  field  which  uses  the  FORTRAN 
fixed  point  decimal  number,  and  the  I  field  that  used  the  FORTRAN 
integer  number  which  is  always  right  justified.  The  card  numbers 
are  not  punched  on -.the. data  cards. 


First  Data  Set 

This  data  set  is  the  one  that  describes  the  surface  equation 


Contents  of  Card 


Card  No. 


The  first  card  is  a. header  card  and  may  contain 
any  alphanumeric  description  of  the  problem  in 
Columns  1-50 


Format 

F15.8 

F15.8 

F15.8 

Variable 

.  xo 

ZO 

FINTX 

Format 

F15.8 

F15.8 

F15.8 

Variable 

SEGX 

SEGZ 

KS 

Format 

F15.8 

F15.8 

Variable 

SEGXS 

XSEG 

Format  F15.8 
Variable  A 


1 

F15.8  15  15  2 

FINTZ  M  N 

3 


next 

KS  cards 
next 

(M+2 ) (N+2 )  cards 
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Format 

F15.8 

next 

Variable 

Z 

(M-l) cards 

Format 

F15.8 

next 

Variable 

X 

(N-l)  cards 

This  completes  the  first  data  set. 

The  coefficients  of  the  surface  equation  (A)  are  presented 
in  the  following  order  (see  Fig.  D-2) : 


‘00  »  A01 

’  A02 

*  A03 

*  A04  * 

»  A10  * 

A11  ’  A12  ’ 

A13  * 

14  ’  **• 

’  A20 

’  A21 

»  A22  *  A23 

’  A24  » 

•”  *  A31  , 

A32  , 

A33  »  A34  »  ‘ *  *  A  40  *  . 

This  is  generally  the  order  in  which  they  will  be  received 
from  the  L.P.  Coefficient  Aqq  is  the  constant  term  for  the 
equation.  It  is  equal  to  the  preliminary  offset  of  the  first 
point  on  the  surface  (Xq  ,  Zq)  .  Any  coefficients  that  weren't 
in  the  final  solution  must  have  a  zero  value  included  in  the 
data  deck. 

Second  Data  Set 

This  data  set  contains  the  packets  of  cards,  each  of  which 
describes  some  information  to  be  extracted  from  the  surface. 

Any  number  of  these  packets  may  be  used  in  a  problem.  They 
may  be  entered  in  any  desired  order  by  simply  stacking  them 
behind  the  first  data  set  in  the  card  reader.  A  description 
of  each  data  packet  follows.  The  first  card  for  each  of  these 
packets  is  an  identification  card  and  contains  a  specific 
value  of  IDENT  in  Column  4  . 
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Packet  for  Offsets  on  a  Waterline 


<<■ 


Contents  of  Card  Card  No. 


The  first 

card 

is  the  identification  card  and 

contains 

a  1 

in  Column  4. 

1 

Format 

14 

F15.8 

2 

Variable 

K 

ZT 

Format 

12 

13  3X  F15.8 

next 

Variable 

NN 

ISP  3  blank  ZE 

K  cards 

Cols . 


(This  set  of  K  cards  contains  the  identification  and 
location  of  any  specific  points.  They  are  to  be 
arranged  so  the  smallest  ZE  is  the  first  card  pro¬ 
gressing  in  order  to  the  largest  ZE  in  the  last  one) 


Packet  for  Offsets  of  a  Station 


Contents  of  Card  Card  No. 


This  card 

has 

a  2 

in  Column  4 

1 

Format 

14 

F15 

.8 

2 

Variable 

K 

XT 

Format 

12 

13 

3X 

F15.8 

next 

Variable 

NN 

ISP 

3  blank 

ZE 

K  cards 

Cols . 


(Begin  with  the  smallest  ZE  and  progress  to  largest) 


Packet  for  Offsets  of  a  Diagonal  Plane 


Contents  of  Card  Card  No. 

This  card  contains  a  3  in  Column  4  1 

Format  F15.8  2 

Variable  ANGLE 

Packet  for  Offsets  of  a  Buttock 

Contents  of  Card  Card  No. 

This  card  contains  a  4  in  Column  4  1 

Format  F15.8  2 

Variable  BUTTK 
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Packet  for  Standard  Cubic  Coefficients  on  a  Waterline 

Contente  of  Card  Card  Ko. 

This  card  contains  a  5  in  Column  4  1 

Format  715.8  2 

Variable  XT 

Packet  for  Standard  Cubic  Coefficients  on  a  Station 

Contents  of  Card  Card  No. 

This  card  contains  a  6  in  Column  4  1 

Format  F15.8  2 

Variable  XT 

Packet  for  Segmentation  of  a  Waterline 

Contents  of  Card  Card  Wo. 

This  card  contains  a  7  in  Column  4  1 

Format  F15.8  F15.8  F15.8  F15.8  2 

Variable  ZT  SX  FX  TOL 

Packet  for  Segmentation  of  a  Station 

Contents  of  Card  Card  No. 

This  card  contains  an  8  in  Column  4  1 

Format  F15.8  F15.8  F15.8  F15.8  2 

Variable  XT  SZ  FZ  TOL 


Output  Data 

The  output  data  for  the  different  types  of  information  requested  is 
given  below.  When  this  output  data  consists  of  offsets  of  curves 
which  may  need  to  be  plotted,  GOBACK  1  includes  cards  which  place 
the  plotter  pen  up  or  down  at  the  proper  times.  The  data  deck  is 
arranged  in  order  and  punched  in  such  a  format  that  the  deck  can 
be  directly  plotted  using  the  plotting  program  of  Appendix  F. 


1.0. 0-2 


D-26 


* 

4- 

Offsets  of  a  Waterline 


The  first  N+2  cards  punched  when  the  offsets  of  a  waterline 
are  called  for  contain  the  coefficients  for  the  two-dimensional 
Thielheimer  equation  of  the  waterline.  Following  this  is  a  header 
card  and  finally  cards  each  containing: 

(1)  A  waterline  identification  number  (WL) 

(2)  The  x  coordinate  of  the  offset 

(3)  The  offset  on  the  waterline 

(4)  The  first  derivative  of  the  waterline 

(5)  The  second  derivative  of  the  waterline 


Offsets  of  a  Station  or  Frame 


Information  corresponding  to  that  of  a  waterline  is  punched.  The 
frame  identification  (FR)  and  the  z  coordinate  of  the  offset 
are  given. 

Offsets  of  a  Diagonal  Plane 

First,  a  header  card  is  punched  and  then  cards  containing  the 
following  data: 

(1)  An  identification  of  the  diagonal  plane  (DP) 

(2)  The  x  coordinate  of  the  offsets 

(3)  The  z  coordinate  of  the  offset 

(4)  The  offset  y^  as  shown  in  Fig.  D-3 

^rfsets  of  a  Buttock 

First,  a  header  card  is  punched  and  then  cards  containing  the 
following : 

(1)  An  identification  of  the  buttock  (BK) 

(2)  The  x  coordinate  of  the  point  is  given 

(3)  The  z  coordinate  of  the  point  is  given 
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Standard  Cubic  Coefficients  for  a  Waterline 


The  first  card  punched  contains  an  identification  of  the  waterline. 
Following  this  are  a  set  of  N  cards,  each  of  which  contains  the 
following  information  for  an  interval  along  the  waterline: 

(1)  The  coefficients  A,B,C,D,  of  the  cubic 

(2)  The  starting  and  ending  x  coordinates  of  the  interval. 

Standard  Cubic  Coefficients  of  a  Frame 

The  same  data  is  given  as  for  a  waterline,  except  that  there  are 
M  intervals  and  the  z  coordinates  of  the  start  and  finish  of 
an  interval  are  given. 

Segmentation  of  a  Waterline  to  a  Tolerance 

First,  a  header  card  is  punched.  Then  cards  containing  the 

following  information  are  punched: 

(1)  The  identification  of  the  waterline 

(2)  The  x  coordinate  of  the  offset 

(3)  The  offset 

(4)  The  tolerance 

Segmentation  of  a  Station  to  a  Tolerance 

The  same  information  is  punched  as  for  a  waterline,  except  the 
z  coordinate  of  the  offset  is  given. 


Sense  Switches 


Sense  Switch  4  ON  -  No  second  derivatives  are  punched  for 

waterlines  or  stations 

Sense  Switch  4  OFF  -  Whenever  an  offset  is  punched  on  a  waterline 

or  station  the  second  derivative  is  punched 

All  other  sense  switches  are  ignored. 
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B.  SAMPLE  PROBLEM 


The  sample  problem  is  a  seven  station  by  eight  waterline  surface 
from  the  DE-1037  class  ships.  The  stations  were  5  through  11  and 
the  water lines  were  the  4'  through  32'.  The  surface  was  single 
splined  in  the  z  direction. 

Offsets  were  calculated  for  one  frame,  one  waterline,  and  one 
buttock.  Standard  cubic  equations  were  found  for  one  waterline 
and  one  frame,  and  the  heights  and  halfbreadths  of  the  intersection 
of  the  01  level  and  the  surface  were  found. 
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B.  SAMPLE  PROBLEM  FOR  GOBACK  1 


*  *  SAMPLE  INPUT  *  * 


DE  1037 


87.5 
4 

20. 

2.5 

1  . 

41.5  _ , 

7.395835 
3.4021024 
.43256810 
-.23899897 
-.05975052 
.30222755 
8.5007951 
-45.380311 
86.461133  l 
-37.816704 
45.392666 
23.518814 
-12.388975 
164.50268 
-298.36348 
129.49442 
-154.92941 
-79.080046 
8.2593167 
-139.79112 
252.23775 
-109.53571  s 
131.40512 
64.742903 
-8.8703614 
176.86673 
-320.34323 
139.37182 
-167.77884 
-77.970667 
1 .4378967 
-46.174359 
85.877451 


4.  35. 

4.  3 

107.5 

150. 

151. 


Coefficients  of  the  surface  equation 


8. 


8 
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< 


-37.57123 
45.564253 
14.238035 
-1 .0426593 
12.468567 
-22.614727 
9.4637234 
-10.732676 
.15734228 
.72400226 
-6.8605798 
9.8515627 
-3.6732793 
4.0526781 
-3.8607566 
-.79698505 
4.5303580 
-6.1544445 
2.3878876 
-3.2185290 
3.8711793 
.(‘6938527 
12.769376 
-20.364931 
8.1774277 
-9.6133634 
1 .0232311 

8. 

12. 

16. 

20. 

24. 

28. 


Waterline  where  new 


coefficients  are  to  be  added  (z) 


32. 

122.5 

157.5 

192.5 


1 


- 1 - 

4 

30. 

r .  i 

107.5 

F  2 

115. 

"0  3 

150. 

F  4 

151 . 

O 

/. 

3 

145. 

•;t  3 

25. 

:r  2 

27. 

Stations  where  new  x  coefficients  are  to  be  added 
_  First  Data  Set 

I  4 


1 

Second  Data  Set 


(x) 
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20.52666 


.)  i  1 


3 

45.'- 

4 

8.  i 

5 

30. 

6 

1 45 . 

7 


30. 

8 

130. 

170. 

.001 

1  45 . 

17. 

26. 

.001 

9 

J-1  LEVEL 

30.56853 

-.1967334 

-.13065 

.06418331 

-.■72799998 
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*  *  SAMPLE  OUTPUT  *  * 


DE  1037  - 


15 .08539200 
2.56392320 
.00166476 
-.29505930 
.53263620 
-.45633232 


1  DENT . 

V 

A 

WL 

30.  CO 

87.50000000 

PEN 

DOWN 

5000.00000000 

FR 

1 

107.49999000 

FR 

2 

114.99999000 

'WL 

30.00 

122.50000000 

WL 

30.00 

127.49999000 

'WL 

30.00 

147.49999000 

FR 

3 

149.99999000 

FR 

4 

150.99999000 

WL 

30.00 

157.50000000 

WL 

30.00 

167.49999000 

WL 

30.00 

187.49999000 

WL 

30.00 

192.50000000 

PEN 

UP 

6000.00000000 

GO 

TO 

87.50000000 

13.07689200 
1 1  .68305200 
-14.38853300 
8.61199300 
-9.49738800 
2.87433700 
-2.4841 1930 
.83143300 
-.23205720 
-.21185000 
IDENT.  Z 

FR  145.00  4.00000000 

PEN  DOWN  5000.00000000 
FR  145.00  8.00000000 


Y 

15.08539200 

16.49597900 

16.95781000 

17.35592000 

17.57887700 

18.19321000 

18.24735500 

18.26805900 

18.39205900 

18.55238000 

18.66692400 

18.63029900 

18.63029900 


Y 

13.07689200 

16.39778400 


FIRST  DER. 
.07325494 


.04805928 
.040  46389 
.01573200 


.01793654 
.03746109 
.1348681 1 
.17204520 


FIRST  DER. 
1 .46038150 

.46918918 
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SECOND  DER. 
.00000271 


-.00144247 

-.00127623 

-.00061129 


-.00027882 

-.00058495 

-.00119721 

-.00135028 


SECOND  DER. 
-.44964165 

-.04595448 


12.00000000 

16.00000000 

20.00000000 

24.00000000 

25.52666000 
27.02666000 
28.00000000 

28.52666000 
32.00000000 

6000.00000000 

4.00000000 

X 

87.50000000 

5000.00000000 

107.49999000 

122.50000000 

142.49999000 

157.50000000 

177.49999000 

192.50000000 
6000.00000000 

87.50000000 

X 

87.50000000 

5000.00000000 

107.49999000 

122.50000000 

142.49999000 

157.50000000 

177.49999000 

192.50000000 
6000.00000000 

87.50000000 


FR  145.00 
FR  145.00 
FR  145.00 
FR  145.00 
ST  3 
ST  2 
FR  145.00 
ST  1 
FR  145.00 
P£fi  UP 
SO  TO 
I  DENT. 

DP  45.00 

PEK  down 

OP  45.00 
OP  45.00 
DP  45.00 
DP  45.00 
DP  45.00 
DP  45.00 
PEN  UP 
GO  TO 
IDENT. 

BK  8.00 
PEN  DOWN 
BK  8.00 
BK  8.00 
BK  8.00 
BK  8.00 
BK  8.00 
BK  8.00 
PEN  UP 
GO  TO 
WL  30.00 
A 

15.08539200 
14.55275600 
18.20341400 
FR  145.00 
A 

13.07689200 
14.26406500 
1 1 .38972800 
19.77363000 
13.12216600 
16.74805900 
22.46800900 


B 

.07325494 
.  1  1 890948 
-.03  75  4  730 

B 

1.46038150 

.57000137 

1.64787770 

-.44809775 

.79905175 

.25516768 

-.45982606 


17.79623100 

18.15465200 

18.29494700 

18.26018300 

18.22452200 

18.18990000 

18.16960000 

18.15974900 

18.10387900 

18.10387900 

y 

17.30200700 

20.43979100 

22.77103300 

25.27905000 

26.59382700 

27.78706900 

28.17376800 

12.07812400 

Z 

4.64062500 

4.00000000 

4.00000000 

4.00000000 

4.00000000 

4.00000000 

4.00000000 

4.00000000 

c 

.00000135 

-.00130305 

.00093204 

c 

-.22482082 

-.00222579 

-.13696034 

.03770428 

-.04024256 

-.01304835 

.01674304 


.64755556 

-.09573533 

-3.87459450 

-13.09499600 


-30.13434200 

-57.38450500 


Z 

19.76562500 

17.54687500 

15.89843700 

14.12500000 

13.19531200 

12.35156200 

12.07812400 


D 

-.00000688 

.00000554 

-.00000510 

D 

.01682029 

-.00172928 

.00388465 

-.00096714 

.00065674 

.00020351 

-.00021025 


-.08745737 

.00577428 

-.01743714 

-.00167504 


.00320919 

-.00183689 


START  FINISH 
87.500  122.500 

122.500  157.500 

157.500  192.500 

START  FINISH 
4.000  8.000 

8.000  12.000 
12.000  16.000 
16.000  20.000 
20.000  24.000 
24.000  28.000 

28.000  32.000 
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l DENT. 

X 

Y 

TOLERANCE 

WL 

30.00 

129.99999000 

17.67813200 

.00100000 

WL 

30.00 

132.63960000 

17.77484800 

.00100000 

WL 

30.00 

135.38715000 

17.86734600 

.00100000 

WL 

30.00 

138.26559000 

17.95605100 

.00100000 

WL 

30.00 

141.30284000 

18.04141700 

.00100000 

WL 

30.00 

144.54118000 

18.12415200 

.00100000 

WL 

30.00 

148.04749000 

18.20537800 

.00100000 

WL 

30.00 

151.93962000 

18.28706100 

.00100000 

WL 

30.00 

156.47632000 

18.37354600 

.00100000 

WL 

30.00 

157.50000000 

18.39205900 

.00100000 

WL 

30.00 

162.26849000 

18.47386700 

.00100000 

WL 

30.00 

166.32337000 

18.53596100 

.00100000 

WL 

30.00 

169.96221000 

18.58406100 

.00100000 

WL 

30.00 

169.99999000 

18.58451700 

.00100000 

IDENT. 

Z 

Y 

TOLERANCE 

FR 

145.00 

17.00000000 

18.19557000 

.00100000 

FR 

145.00 

18.38721400 

18.25104600 

.00100000 

FR 

145.00 

19.25546800 

18.27884000 

.00100000 

FR 

145.00 

19.98029800 

18.29463200 

.00100000 

FR 

145.00 

20.00000000 

18.29494300 

.00100000 

FR 

145.00 

20.70570100 

18.30189200 

.00100000 

FR 

145.00 

21 .48541600 

18.30114000 

.00100000 

FR 

145.00 

22.38815800 

18.29159400 

.00100000 

FR 

145.00 

23.56892300 

18.26968400 

.00100000 

FR 

145.00 

24.00000000 

18.26017000 

.00100000 

FR 

145.00 

24.58768000 

18.24666900 

.00100000 

FR 

145.00 

26.00000000 

18.21334600 

.00100000 

0-1  LEVEL 


X  Z  Y 

87.50000000  30.56853000  15.28253600 

107.49999000  30.42542500  16.59634800 

122.50000000  30.30533000  17.39770100 

127.49999000  30.26864300  17.60595300 

147.49999000  30.14414200  18.18858100 

157.50000000  30.09313100  18.38372500 

167.49999000  30.04800600  18.54872600 

187.49999000  29.97059000  1.8.67598500 

192.5COOOOOO  29.95303000  18.63659300 
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IDNT 


i-2 
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< 


*  52 


K3=0 


[iWij 


-2 


iWITCH  4 


r  PUNCH 
CTjXZjY.RDERT 


415 
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XZx  OLD+  SEG 
OLD»  X2 


1.0. 0-2 


n-AA 


1.0. 0-2 


D-47 


-2 


D-48 


1.0. 0-2 


D-49 


irwr 


D-50 


IK--TK+1 


L 


J 


IND-IND  +  1 
KZ*0 

XZ-ZX(IND) 
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PINT  =  PlNTZ 
SE6Z 
piNTZ 
R0R16  -  XORlG 
TOT  =TOTz: 


i. 0.0-2 
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f- 

'4* 


1.0. 0-2 


503 


D-55 


807 


*2.81 


I. 0.0-2 
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1.0. 0-2 
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l.o.o-'’ 


D-59 
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D.  LISTING  -  GG3ACK  1 


C 

101 

102 

103 

104 

105 
DER 


SURFACE  GO  BACK  (NO  PROFILE) 

FORMAT (4F1 5 »8. 215) 

F0RMAK4F15.8) 

F0RMAT(5HG0  T0,4X,F14.8,3X,F15.8) 

FORMAT  ( 2F15 .8,  15)  _ 

FORMAT (7H  IDENT.,9X,1HZ,17X,1HY,11X,10HFIRST 


DER. ,4X, 1 1HSEC0ND 


DER. ,4X,  1 1HSEC0ND 


C.) 

106  F0RMAT(I4,F15.8) 

107  F0RMAT(2HWL,F7.2,F14.8,3X,3F15.8) 

108  F0RMAT(8HPEN  DOWN, FI 5 . 8 ) 

109  F0RMAT(2HFR,F7.2,F14.8,3X,3F15.8) 

110  F0RMAT(6HPEN  UP,2X,F15.8) 

111  F0RMAT(2HLG, I3,3X,F15.8,3X,F15.8) 

112  FORMAT (7H  I  DENT. ,9X, 1HZ , 1 7X , 1HY , 1 2X .9HT0LERANCE) 

113  FORMAT ( 7H  I  DENT. ,9X, 1HX, 1 7X, 1HY.1 2X,9HT0LERANCE) 

114  FORMAT ( 2HDP ,F6.2, F15 .8, 3X, 2F15 .8) 

115  FORMAT ( 7H  I  DENT. ,9X, 1HX, 1 7X, 1HY , 1 1 X, 10HF I RST 
DER 

1 16° FORMAT ( 1X,6HIDENT.,9X,1HX,17X,1HY,7X.14HFIRST  DER I  VAT.) 
lie  DERIVAT.) 

!12S}BHxKi^:!^i^:(l5c:?lx.,HO,9X.5HSTART,8H  FINISH) 

121  FORMAT ( 7H  I  DENT. . 10X, 1HX, 1 7X, 1  HZ) 

122  FORMAT (4F1 4.8, 2F8. 3)  . 

123  FORMAT (50H  .  > 

124  FORMAT (25H  x  ) 

DIMENSION  SEGXS(5) ,XSEG(5) ,CZ(42) 

100  READ  123 
PUNCH  123 

READ  101 ,XOR I G,ZOR I G, FINTX.FI NTZ,M,N 
READ  104,SEGX,SEGZ,KS 
I F(KS-1 ) 1 70, 1 71 , 1  72 

171  READ  102 , SEGXS( 1 ) , XSEG( 1 ) 

GO  TO  170 

172  DO  173  1-1, KS 

173  READ  102, SEGXSC I ) , XSEG( I ) 

170  NP-N+2 

MP-M+2 

FIVE-5000. 

SIX-6000. 

DO  201  1-1, MP 
DO  202  J-1 , NP 
202  READ  102,  A( I ,J) 
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201  CONTINUE 
M-MP-3 
N-NP-3 

DO  203  1-1, M 
READ  102, Z(l) 

IF(M-I ) 230, 232, 203 

232  TOTZ-Z(I) 

203  Z( I )-(Z( I ) — ZOR I G ) / F l NTZ 
230  DO  204  l-1,N 

READ  102, X(l) 

I F( N- I ) 205 , 233 , 204 

233  TOTX-X(I) 

204  X( I )-(X( I ) — XOR 1 G ) / F I NTX 

205  READ  106,1 DNT 
K2-0 

K3-1 

KSR-KS 

XZ-O. 

KID-1 

IND-O 

OLD-O. 

I F( I DNT— 1 ) 600, 209, 600 

600  I F( IDNT-2)601 ,310,601 

601  I F( I DNT-3 )602 ,231,602 

602  IF( IDNT-4)603, 234,603 

603  I F ( I DNT-5 )604, 801 ,604 

604  I F( IDNT-6)605 ,831 ,605 

605  IF(IDNT-7)606,800,606 

606  I F ( I DNT-8) 700,820 , 700 

C  *  *  A)  OFFSETS  OF  A  WATER  LINE  *  * 

209  READ  106,K,ZT 

I F ( K— 1 ) 200, 206, 207 
207  DO  210  1-1, K 

READ  1 1 1 , I SP( I ) , XE ( I ) 

210  XE( I )-(XE ( I )-XOR I G) /FI NTX 
GO  TO  200 

206  READ  111,  I SP( 1 ) , XE ( 1 ) 

XE ( 1 )- ( XE ( 1 )-XOR 1 G) /F I NTX 
200  FINT-FINTX 

SEG-SEGX/F I NTX 
RORIG-XORIG 
TOT-TOTX 
IK-0 

804  INX-NP 

DO  236  1  =  1,  N 
236  ZX(I)-X(I) 

ZZ=(ZT-ZOR I G ) / F I NTZ 
740  DO  211  1-1, M 


1.0. 0-2 


D-62 


MAZ-l-1 

I F(ZZ-Z( l))212,212,211 

211  CONTINUE 

212  I F( IDNT-9) 721 , 722, 721 
721  DO  213  1-1. INX 

213  C( I )-A(1, I )+A( 2, I )*ZZ+A(3, I )*ZZ*ZZ+A(4, I )*ZZ**3 
I F ( MAZ- 1)221,215,216 

215  DO  217  1-1,  INX 

217  C(l)-C(l)+A(5, l)*(ZZ-Z(1))**3 
221  I F( IDNT-5) 214,503,503 

216  I NZ-MAZ+4 

DO  218  1-1, INX 
DO  219  J-5.INZ 

219  C( l)-C(l)+A(J, l)*(ZZ-Z(J-4))**3 

218  CONTINUE 

I F( IDNT-5) 214,503,503 

214  INZ-INX 
GO  TO  314 

C  *  *  B)  OFFSETS  OF  A  FRAME  OR  STATION  *  * 

310  IK-0 

READ  106,K,ZT 
I F ( K— 1 ) 90 1 ,403.404 

403  READ  111, I S P ( 1 ) ,XE( 1) 

XE( 1 )-( XE( 1 ) — ZOR I G ) / F I NTZ 
GO  TO  901 

404  DO  400  1-1 ,K 

READ  111,1 SP( 1 ) , XE( I ) 

400  XE( I )«(XE(  1 )— ZOR I G ) /F I NTZ 

901  FINT-FINTZ 
SEG-SEGZ/FINTZ 
RORIG-ZORIG 
TOT-TOTZ 

903  DO  401  1-1, M 

401  ZX(I)-Z(I) 

XX«( ZT-XORIG) / F I  NT X 
360  DO  311  1-1, N 
MAX- 1-1 

IF(XX~X(  0)312,312,311 

311  CONTINUE 

312  INZ-MP 

DO  313  1-1, INZ 

313  C( l)-A( I .1)+A( l,2)*XX+A( I ,3)*XX*XX+A( I ,4)*XX**3 
I  F  ( MA  X- 1 )  90  2 , 3 1 5 , 3 1 6 

315  DO  317  1-1, INZ 

317  C( I )-C( I )+A( I ,5)*( XX-X( 1 ) )**3 

902  I F( ID NT-5) 31 4, 906,906 
906  I F( IDNT-9) 503 , 702, 503 

316  INX-MAX+4 

DO  318  1-1, INZ 
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DO  319  J=5» I  NX 

319  C(l)-C(l)+A( l,J)~(XX-X(J-4))**3 
318  CONTINUE 

I F ( I DNT-5) 31 4,906,906 
314  I F( I D NT— 3 ) 905 *321  ,321 
90S  DO  320  1-1 , INZ 

320  PUNCH  102, C( I ) 

I F (SENSE  SWITCH  4)322,323 

322  IF(!DNT-2)181,182,182 

181  PUNCH  116 
GO  TO  405 

182  PUNCH  118 
GO  TO  405 

323  I F( IDNT-2) 1 83, 1 84, 184 

183  PUNCH  115 
GO  TO  405 

184  PUNCH  105 

405  Y-C( 1 )+C(2)*XZ+C(3)*XZ*XZ+C(4)*XZ**3 
TDERI«C(2)+C(3)*2.*XZ+C(4)*3.*XZ*XZ 
SDER l-C(3)*2.+C(4)*6.*XZ 

I F( I ND-1 ) 406,408,408 

408  DO  409  1-1, 1  NO 
DIF-XZ-ZX(I) 

I F(D I F-. 000001 ) 406,406,491 
491  Y-Y+C ( l+4)*(D I F)**3 

TDER l-TDERI+3.*C( l+4)*(D I F)**2 
SDER!=SDERI+6.*C(I+4)*DIF 

409  CONTINUE 

406  XZ-XZ*F I NT+RORIG 
RDERI-TDERI/FINT 
SDERI-SDERI/(FINT*FINT) 

I F(K2)438,439,438 

439  I F ( SENS E  SWITCH  4)410,411 

410  I F ( I D NT- 2) 437»436 ,436 

437  PUNCH  107,ZT,XZ,Y,RDERI 
GO  TO  412 

436  PUNCH  109,ZT,XZ,Y,RDERI 
GO  TO  412 

411  I F( IDNT-2) 441, 440, 440 

441  PUNCH  107,ZT, XZ,Y,RDERI , SDER I 
GO  TO  412 

440  PUNCH  109»ZT, XZ,Y,RDERI , SDER I 
GO  TO  412 

438  PUNCH  111. I S P ( IK) ,XZ,Y 

412  IF(XZ-TOT)497,51 ,51 
51  PUNCH  110, SIX 

PUNCH  103,R0RIG,Y 
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GO  TO  205 

497  IF(K3)52.415,52 
52  PUNCH  108, F I VE 
K3-0 

415  i F(KSR-1 ) 1 74,1 75, 1 76 

175  I F( XSEG(1 )-XZ) 1 74, 1 78, 1 74 

178  SEG-SEGXS(1)/FINTX 
KSR=0 

GO  TO  174 

176  DO  177  I-1.KSR 

! F( XSEG( I )-XZ) 1 77, 1 79, 1 77 

179  SEG-SEGXS(  D/FINTX 

I F ( KSR— I ) 1 77, 1 80, 1 77 

180  KSR=0 

177  CONTINUE 

174  I F( IDNT-3) 186,185,185 

186  !F(OLD+SEG-ZX( IND+1 ) )416 ,41 7,41 8 

418  IF(K)419,419,420 

420  I F(K-K I D)419, 421,421 

421  IF(ZX( IND+1 )-XE (KID) >419,422, 423 

423  XZ=XE (K I D ) 

K2=  1 

I F(K-K 10)405,424,424 

424  K I D=K I D+1 
IK-IK+1 
GO  TO  405 

422  I ND= I ND+1 
GO  TO  423 

419  I ND= I ND+1 
K2=0 

XZ=ZX( i NO) 

GO  TO  709 

417  IF(K)430,430,428 

428  IF(K-KID)430,429,429 

429  I F(ZX( I ND+1 )-XE (KID) >430,431,423 

431  OLD=OLD+SEG 
GO  TO  422 

430  OLD=OLD+SEG 
GO  TO  419 

416  I  F ( K ) 43  2 , 43  2 , 43  3 

433  1 F ( K-K 10)432,434,434 

434  I F(OLD+SEG-XE (KID) >432,435, 423 

435  OLD=OLD+SEG 
GO  TO  423 

432  XZ-OLD+SEG 
Ol.D=OLD+SEG 
K2=0 

709  I F( I DHT-9) 405 , 710,405 

710  XX-XZ 
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GO  TO  360 

C  *  *  E)  STANDARD  CUBIC  COEFFICIENTS  FOR  A  WATER  LINE  *  * 
801  READ  102, ZT 
PUNCH  107, ZT 
PUNCH  120 
MA*»i\ 

ROR I G-XOft  I G 
F I  NT=F  i  N'TX 
GO  TO  804 

503  IF( I DNT-7) 81 1 ,810,810 

810  DO  805  1-1, NA 
I  N  X<=  I  —1 

IF(S-ZX(  0  )  506,805,805 

805  CONTINUE 
GO  TO  506 

811  I  NX-0 

506  XE(1)-C(1) 

X E(2)-C(2) 

:<e(3)=c(3) 

XE (4)=C (4) 

i F( I  NX-1 ) 508, 50 7, 50 7 

507  DO  504  I-1.IHX 

XE  ( 1  )*=XE  ( 1  )-C  ( !+4)*ZX( I )**3 
XE(2)-XE(2)+3.*C( l+4)*ZX( l)*ZX( I ) 

XE(3)»XE(3)-3.*C( l+4)*ZX( 1 ) 

XE(4)=XE(4)+C( 1+4) 

IF( l-l NX) 504,508, 508 

504  CONTINUE 

508  ;f( I DNT-7) 807, 806, 806 

806  XS=S 

I F ( XS-F ) 23 , 22, 22 

22  XS=F 

23  V=XE ( 1 ) +XE (2)*XS+XE (3 )*XS*XS+XE (4)*XS*XS*XS 
!F( IDNT-8)840,841 ,841 

840  XT«XS*FI MTX+XORIG 
PUNCH  107, ZT, XT, Y , TOL 
GO  TO  842 

841  XT=XS*F I NTZ+ZOR ! G 
PUNCH  109 , ZT, XT, Y , TOL 

842  I F( XS-F) 34, 205, 205 
34  AA=3.*XE(4)*S-XE(3) 

AC=XE(3)*S*S+XE(4)*S*S*S 

AAA-2. *XE(4)*XS*XS*XS-AA*XS*XS-2.*XE(3)*S*XS+AC 
DER«(XE(2)+2.*XE(3)*XS+3.*XE(4)*XS*XS)/FINTX 
Y=AAA*AAA-TOL*TOL*( 1 ,+DER*DER) 

I F(Y)4, 12,4 

4  H**2  .*TOL 

5  CA-Y 
XS-XS+H 
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I F ( XS— ZX ( INX+1 )) 7.7,44 
44  S-ZX( INX+1) 

GC  TO  810 

7  AAA-2. *XE(4)*XS*XS*XS-AA*XS*XS-2.*XE(3)*S*XS+AC 
DER«(XE(2)+2.*XE(3)*XS+3.*XE(4)*XS*XS)/FINTX 
Y«AAA*AAA-T0L*T0L*(1 .+DER*DER) 

I F(CA*Y)9» 1 2,5 

9  IF(ABS(H)-. 0001)10, 11, 11 

11  H--H/2. 

GO  TO  5 

10  XS-XS-H/2. 

12  RAD«8.*XE(3)*XE(4)*XS+12.*XE(4)**2*XS**2-2.*XE(3)*XE(4)*S 
RAD=SQRT(RAD-3 .*XE (4)*XE (4)*S*S+XE (3)*XE (3) ) 

13  XX— .5*(S+(XE(3)+RAD)/XE(4)) 

I F( XX-S) 14,15,15 

14  RAD— RAD 
GO  TO  13 

15  I F( XX-ZX( I NX+1 ) ) 1 9 .44,44 
19  I F( XX-F) 16,17,17 

16  S=XX 

GO  TO  806 

17  S-F 

GO  TO  806 

807  IF( INX)280, 280,281 

280  S-RORIG 

F=ZX( 1 )  *F  I NT+ROR I  G 
GO  TO  282 

281  S=ZX( I  )*F I NT+ROR I G 
F=ZX( I +1 )*F I NT+ROR I G 

282  DO  283  1=2,4 

283  XE( I )«XE( I )/FINT**( 1-1) 

PUNCH  122,XE(1),XE(2),XE(3),XE(4),S,F 
!F( INX-NA+2)509,509,205 
509  I NX= I NX+1 
GO  TO  506 

*  *  H)  SEGMENTATION  OF  A  FRAME*  * 

820  READ  102,ZT,S,F,T0L 

S=(S-ZOR I G)/F I NTZ 
F=( F-ZOR I G) /F I NTZ 
NA-M 

PUNCH  112 
GO  TO  903 

*  *  F)  STANDARD  CUBIC  COEFFICIENTS  FOR  A  FRAME  ** 

831  READ  102, ZT 

NA-M 

PUNCH  109, ZT 
PUNCH  120 
RORIG-ZORIG 
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402 

C 

800 


c 

234 


C 

231 


1360 


321 

452 


454 

455 

456 
465 

ill 

480 


478 
453 

479 

481 

482 
471 


FINT-FINTZ 
I NZ-MP 
GO  TO  903 

*  *  G)  SEGMENTATION  OF  A  WATER  LINE  *  * 
READ  102,ZT,S,F,T0L 
S-(S-XORIG)/FINTX 

F-  ( F— XOR I G ) / F I  NT X 
NA-N 

PUNCH  113 
GO  TO  804 

*  *  D)  OFFSETS  OF  A  BUTTOCK  *  * 

READ  102.BUTTK 

PUNCH  121 
GO  TO  1360 

*  *  C)  OFFSETS  OF  A  DIAGONAL  PLANE  *  * 
READ  102,  ANGLE 

TAN-ANGLE*3. 14159/180. 

SI-SIN(TAN) 

CO-COS (TAN) 

TAN— S  I  /CO 
PUNCH  117 
SEG-SEGX/F I  NTX 
SEGTZ-SEGZ/FINTZ 
XX-O. 

I NZ-MP 
GO  TO  360 
K-1 
KID-1 

ZXZ-XZ*F I NTZ+ZOR I G 

Y-C ( 1 ) +C ( 2 ) *XZ+C ( 3 ) *XZ*XZ+C ( 4 ) *XZ**3 

IF(  IND-D465.454.455 

Y-Y+C(5 )*( XZ-Z( 1 ) )**3 

GO  TO  465 

DO  456  1-1, 1  NO 

Y-Y+C(I+4)*(XZ-Z(I))**3 

IF( lDNT-4)478,477,477 

I F(BUTTK) 480,485 , 480 

ZXZ-ZOR I G 

GO  TO  457 

Y-Y-BUTTK 

GO  TO  453 

y-tan*(totz-zxz)+y 
IF(ABS(Y)-. 005)457, 457, 479 
GO  TO  (481, 482), KID 
KID-2 

IF(Y)473,457,457 
GO  TO  (471, 472), K 
I F(Y)473 ,457,474 
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473  IF(ZXZ+SEGZ-T0TZ)483,483,484 
484  ::x~=totz 

,-u  TO  457 
483  OlD-XZ 
x;>xz+segtz 

GO  TO  458 

474  '<-2 

GO  TO  475 

472  IF(Y)476, 457,475 

475  XT=XZ 

XZ*XZ-(ABS( XZ-OLD) )/2. 

OLD=XT 
GO  TO  458 

476  XT=XZ 

XZ=XZ+(ABS( XZ-OLD) )/2. 

OLD=XT 

458  DO  459  I-1.M 

I  F ( XZ— Z ( I )) 452, 45 2, 459 

459  I  NO- 1 

GO  TO  452 

457  XZ=XX*F I NTX+XOR I G 

I F ( I DNT-4) 467,468,468 

467  Y=( TOTZ-ZXZ) /CO 

PUNCH  1 14, ANGLE, XZ,Y,ZXZ 
GO  TO  466 

468  PUNCH  1 19.BUTTK.XZ.ZXZ 
466  I F( XZ-TOTX) 497. 64,64 
64  PUNCH  1 10.SIX 

66  PUNCH  103.X0RIG.ZXZ 
GO  TO  205 

185  I F( XX+SEG-X(K2+1 )) 463 ,462,461 

461  K2=K2+1 
XX=X(K2) 

324  XZ=0. 

I  ND=0 
GO  TO  360 

462  K2=K2+1 

463  XX=XX+SEG 
GO  TO  324 

C  *  *  I)  OFFSETS  OF  LONGITUDINAL  ELEMENTS  *  * 

700  READ  124 
PUNCH  124 
PUNCH  125 

DO  701  1=1 , NP 
READ  102,CZ( I ) 

701  CONTINUE 

DO  707  1=1, N 
707  ZX( 1 )=X( I ) 
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K*=0 

TOT-TOTX 

SEG-SEGX/FIMTX 

XX-O. 

GO  TO  360 

702  ZPCH»CZ(1)+CZ(2)*XX+CZ(3)*XX*XX+CZ(4)*XX**3 
I F(MAX-1 ) 1000, 704,704 

704  DO  705  1-1, MAX 
ZPCH»ZPCH+CZ( l+4)*( XX— X ( I ) )**3 
I F ( MAX- 1)1000,1 000 , 705 

705  CONTINUE 

1000  ZZ* ( ZPCH-ZOR I G) /F I NTZ 
GO  TO  740 

722  Y-C(1 )+C(2)*ZZ+C(3)*ZZ*ZZ+C(4)*ZZ**3 
I F(MAZ-1 ) 703 ,724, 724 

724  DO  725  ! » 1 , MAZ 

Y*Y+C( l+4)*(ZZ-Z( I )  )**3 
I F ( MAZ- I ) 703 , 703 , 725 

725  CONTINUE 

703  XPCH=(XX*F I NTX)+XOR I G 
PUNCH  101 ,XPCH,ZPCH,Y 

xz-xx 

I  F(XPCH-TOT) 186,205,205 
END 
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Section  VI 


GOBACK  2 
3-D  PROFIIE 


A.  OPERATING  INSTRUCTIONS 

This  version  of  GOBACK  will  accept  coefficients  of  a  surface  equation 
such  as  given  in  Fig.  D-2,  modified  as  shown  in  Section  IV,  for  end 
profile  requirements.  The  equation  must  be  single  splined  in  the  x 
direction  and  may  be  single  or  double  splined  in  the  z  direction. 


The  following  data  can  be  calculated  with  this  program: 

—  Offsets,  first  and  second  derivatives  along  a  waterline 
at  a  given  interval 

—  Offsets,  first  and  second  derivatives  along  a  station 
at  a  given  interval 

—  Offsets  at  a  given  interval  along  buttocks  as  given  in 
Section  IV-B.of  this  appendix. 


Fortran  Input  Symbols 
Symbol 


Definition 


XO 

ZO 

FINTX 


The  actual  full  scale  coordinate  of  the  first 
station  of  the  surface 

The  actual  full  scale  z  coordinate  of  the  first 
waterline 


This  value  is  used  as  a  scale  factor  in  the  x 
direction.  If  the  surface  equation  is  scaled  so 
the  stations  are  one  unit  apart,  FINTX  equals  the 
actual  full  scale  station  interval  (x.-x  ).  If 
the  stations  of  the  surface  equation  are  ?ess  than 
one  unit  apart  FINTX  becomes  some  multiple  of  the 
station  spacing.  For  example,  if  the  stations  in 
the  equation  are  1/4  unit  apart  FINTX  becomes  4  times 
the  actual  full  scale  station  spacing 


>] 
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FINTZ 

M 


N 

K 

NAB 

SEGX 

SEGZ 

POINT 

PWR1 


PWR2 


Same  definition  as  FINTX,  except  in  the  z  direction 


The  total  number  of  points  of  discontinuity  along  a  station 
(points  where  z  coefficients  are  added)  including  the 
points  at  the  first  and  last  water  lines.  For  example, 
if  a  surface  containing  seven  waterlines  were  single 
splined,  M  would  equal  seven  since  there  are  third 
derivative  discontinuities  ac  each  waterline.  If  the 
seven-waterline  surface  were  double  splined,  M  would 
equal  four,  since  there  are  discontinuities  only  at  water- 
lines  0,2,4,  and  6  . 

The  total  number  of  points  of  discontinuity  along  a 
waterline  (points  where  x  coefficients  are  added). 

The  example  used  for  M  is  valid  except  using  stations 
instead  of  waterlines. 


The  number  of  coefficients  in  the  profile  equation 


This  value  tells  the  program  whether  the  surface  has 
been  single  splined  in  the  z  direction  (NAB  ■  +1)  or 
or  double  splined  in  the  z  direction  (NAB  ■  -1) 

Desired  interval  between  consecutive  calculated  offsets 
along  water lines,  buttocks,  and  diagonal  planes. 

Desired  interval  between  consecutive  calculated  offsets 
along  stations. 

The  actual  full  scale  interval  over  which  the  end 
condition  is  to  be  effective.  Usually  less  than  or 
equal  to  the  first  station  interval  (D  in  Section  IV) 


In  the  end  condition  function  of  Section  IV 


T(x,z)  = 


l  .  ( P+gl2)-* 


PWR1  is  the  value  of  n  .It  is  nearly  always  equal 
to  3  . 


Equals  the  power  (P)  to  which  the  function  T(x,z)  is 
to  be  raised. 


CZ  -  The  coefficients  of  the  profile  equations.  These 

coefficients  must  be  scaled  the  same  and  have  the  same 
origin  as  the  surface  equation. 

A  -  The  coefficients  of  the  surface  equation  produced  by 

the  linear  program.  There  are  (M+2)  (N+2)  coefficients. 
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z 


X 

ZP 


IDENT 


z  coordinates  of  waterlines  where  coefficients  are 
added,  plus  that  of  the  last  waterline.  There  are 
(M-l)  since  the  value  for  the  first  waterline  isn't 
included. 

Same  as  Z  except  for  stations  (N-l) 

If  the  surface  equation  has  been  double  splined  in 
the  z  direction  (NAB  ■  -1),  the  z  coordinates  of 
the  waterlines  where  new  coefficients  are  added  must 
be  read  in.  These  are  ZP  and  there  are  (K-3) .  If 
the  equation  is  single  splined  in  the  z  direction 
(NAB  *  +1)  and  these  are  not  necessary. 

Code  which  tells  the  program  what  data  is  required 
from  the  surface 


ZT  -  z  coordinate  of  a  waterline  along  which  data  is  to 

be  calculated 


NBOW  -  An  indicator  for  telling  the  program  if  the  buttocks 

will  be  single  valued  or  multi-valued.  If  NBOW  *  1 
the  case  must  be  that  shown  in  Fig.  D-5a,  Section  IV-B 
If  NBOH  »  2  the  program  will  expect  to  find  the  more 
complicated  cases . 

BUTTK  -  The  y  distance  from  the  centerplane  of  the  ship 

to  a  buttock  along  which  heights  are  to  be 
calculated 


Input  Data  Cards 

There  are  two  distinct  kinds  of  sets  of  data  used  for  input  to 
GOBACK  2  .  The  first  set  consists  entirely  of  data  describing  the 
surface  equation  to  the  program.  No  output  is  produced  from  this 
set.  The  second  set  of  data  consists  of  packets  of  data  cards. 
Each  packet  describes  some  information  required  from  the  surface 
and  gives  the  geometric  data  necessary  to  obtain  the  information 
from  the  surface  equation.  Only  one  of  the  first  sets  of  data  are 
entered  per  problem.  There  may  be  many  small  packets,  each  of 
which  calls  for  offsets  along  a  waterline  or  station  or  perhaps 
standard  cubic  coefficients,  etc. 

There  are  limitations  on  some  of  the  input  parameters  of  the 
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program.  These  limitations  follow: 


Variable 

Minimum 

Maximum 

M  (no.  of  water lines) 

3 

11 

N  (no.  of  stations) 

3 

15 

K  (no.  of  coefficients  in 
profile  equation) 

1 

16 

In  describing  the  various  data  cards,  the  actual  FORTRAN  format 
field  description  is  used  in  most  cases.  These  fields  come  con¬ 
secutively  across  the  card  with  no  gaps  or  blank  columns  between, 
except  where  indicated.  The  field  descriptions  are  the  FORTRAN  F 
field,  which  uses  the  FORTRAN  fixed  point  decimal  number,  and  the 
I  field  that  used  the  FORTRAN  integer  number,  which  is  always 
right  justified.  The  card  numbers  are  not  punched  on  the  data 
cards . 

First  Data  Set 

This  data  set  is  the  one  that  describes  the  surface 
Contents  of  Card 

The  first  card  is  a  header  card  and  may  contain 
any  alphanumeric  description  of  the  problem  in 
columns  1-50  . 


Format 

Variable 

F15.8 

XO 

F15.8 

ZO 

F15.8 

FINTX 

F15.8 

FINTZ 

15 

M 

15  15 

N  K 

15  2 

NAB 

Format 

Variable 

F15.8 

SEGX 

F15.8 

SEGZ 

F15.8 

POINT 

F15.8 

PWR1 

F15.8 

FWR2 

15 

NBOW 

3 

Format 

Variable 

F15.8 

CZ 

next 

K  cards 

Format 

Variable 

F15.8 

A 

next 

(M+2) (N+2) 
cards 

Format 

Variable 

F15.8 

Z 

next 

(M-l)  cards 

Format 

Variable 

F15.8 

X 

next 

(N-l)  cards 

equation 
Card  No. 
1 
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Format  F15.8 

Variable  ZP 


next 

(K-3) 


(This  set  of  K-3  cards  need  be  entered  only  if 
NAB  ■  -1  ,  see  Input  Symbol  Definition) 


This  completes  the  first  data  set. 


The  coefficients  of  the  surface  equation  (A)  are  presented 
in  the  following  order  (See  Fig.  D-2) : 


o 

o 

< 

a 

’  Aoi 

CM 

O 

< 

’  A03 

'  A04  * 

’  A10  ’  A11 

'  A12 

’  A13  * 

> 

1— * 

■p' 

,  . .  • 

’  A20 

»  A21 

‘  A22  *  A23 

»  A24  * 

’  A31 

*  A32  • 

A33 

’  A34 

)  *  *  * 

*  A40 

>  •  •  •  >  ‘  ' 

This  is  generally  the  order  in  which  they  will  be  received 
from  the  L.P.  Any  coefficients  that  weren't  in  the  final 
solution  must  have  a  zero  value  included  in  the  data  deck. 

Second  Data  Set 

This  data  set  contains  the  packets  of  cards,  each  of  which 
describes  some  information  to  be  extracted  from  the  surface. 
Any  number  of  these  packets  may  be  used  in  a  problem.  They 
may  be  entered  in  any  desired  order  by  simply  stacking  them 
behind  the  first  data  set  in  the  card  reader.  A  description 
of  each  data  packet  follows.  The  first  card  for  each  of 
these  packets  is  an  identification  card  and  contains  a 
specific  value  of  IDENT  in  Column  4. 

Packet  for  Offsets  on  a  Waterline 

Contents  of  Card  Card  No. 

The  first  card  is  the  identification  card  and 
contains  a  1  in  Column  4  1 

Format  F15.8  2 

Variable  ZT 
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Packet  for  Offsets  of  a  Station 


Contents  of  Card 
This  card  has  a  2  in  Column  4 

Format  14  F15.8 

Variable  K  XT 

Packet  for  Offsets  of  a  Buttock 
Contents  of  Card 

This  card  contains  a  4  in  Column  4 

Format  F15.8 

Variable  BUTTK 


Card  No. 
1 

2 


Card  No. 
1 
2 


Output  Data 

The  output  data  for  the  different  types  of  information  requested  is 
given  below.  Since  the  output  data  consists  of  offsets  of  curves 
which  may  need  to  be  plotted,  GOBACK  2  includes  cards  which  place 
the  plotter  pen  up  or  down  at  the  proper  times .  The  data  deck  is 
arranged  in  order  and  punched  in  such  a  format  that  the  deck  can 
be  directly  plotted  using  the  plotting  program  of  Appendix  F. 

Offsets  of  a  Waterline 

The  first  N+2  cards  punched  when  the  offsets  of  a  waterline  are 
called  for  contain  the  coefficients  for  the  two-dimensional 
Thielheimer  equation  of  the  waterline.  Following  this  is  a  header 
card  and  finally  cards  each  containing: 

(1)  A  waterline  identification  (WL) 

(2)  The  x  coordinate  of  the  offset 

(3)  The  offset  on  the  waterline 

(4)  The  first  derivative  of  the  waterline 

(5)  The  second  derivative  of  the  waterline 

Offsets  of  a  Station  or  Frame 

Information  corresponding  to  that  of  a  waterline  is  punched.  The 
frame  identification  (FR)  and  the  z  coordinate  of  the  offset 
are  given. 
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Offsets  of  a  Buttock 

First,  a  header  card  is  punched  and  then  cards  each  containing 
the  following: 

(1)  An  identification  of  the  buttock 

(2)  The  x  coordinate  of  the  offset 

(3)  The  z  height  of  the  offset 

(4)  The  point  number  on  the  frame 

If  the  buttock  is  multi-valued  on  a  frame,  the  points  are  numbered 
in  Columns  48  and  49  of  the  card.  The  points  are  numbered  con¬ 
secutively  beginning  with  the  lowest  (least  z  )  value  on  each 
frame.  If  a  buttock  is  multivalued  during  only  part  of  its  length, 
the  series  of  numbers  describing  it  may  change  as  shown  in  the 


All  sense  switches  are  ignored. 
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B .  SAMPLE  PROBLEM 


The  sample  problem  was  taken  from  a  surface  describing  the  bow  of 
the  DLG-26  class  frigate.  The  surface  has  eleven  water lines  and 
five  stations.  It  has  been  double  splined  in  both  the  x  and  z 
directions.  The  offsets  of  one  waterline,  one  station,  and  one 
buttock  have  been  solved  for. 

In  this  case  the  buttock  is  double  valued  overpart  of  its  length. 
The  offsets  for  the  two  lines  must  be  separated  before  plotting. 
This  can  be  done  by  examining  their  z  coordinates. 
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*  * 


SAMPLE  INPUT 


*  * 


1 


*  COMPLETE  BOW  DLG 

0.0  0.0 

13  -1 

1.0  2. 

2  _ 

.8562 

-.272763473 

u.O 

-.6288 
.1531 
-.0704 
-.0041 

0.0 
0.0 
0.0 
0.0 
o.c 
0.0 

7.16765334 
15.10834002 
3.90044334 

0.0 

1 .37436000 
5.17799440 
-  21.63735908 

6.09769242 

0.0 

1.46282033 
2.89710319 
10.92294755 
3.04203734 

0.0 


5.10 

1.275 


8. 


3. 


V  Coefficients  of  the  profile  equation 


^Coefficients  of  the  surface  equation 


.01033209 

2.33844915 

1.72304778 

.46882580 

.26489650 

3.81770188 


6  3 

.3333333 
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0.0 


35464391 

12961316 


.57317627 

1.54632108 

2.35139965 

.65726751 

.0 

.35806726 

.34638152 

.16817131 

.00715955 

.00013420 

.04213764 

.44640127 

.59797708 

.23047406 

0.0 

.13067829 

8.  — i 

16. 

24.  Water  lines  where  new  z  coefficients  are 


32.  r  to  be  added  (z) 

40.  _ J 


» 


*  *  SAMPLE  OUTPUT  *  * 


*  COMPLETE  BOW  DLG  <  - 

.65988975 
1 .68165847 
.08525161 
-.00001677 
-.01417503 


1  DENT. 

X 

Y 

FIRST  DER. 

SECOND  DER. 

WL 

28.00- 

2.29857000 

0.00000000 

0.00000000 

0.00000000 

PEN 

DOWN 

5000.00000000 

WL 

28.00- 

1 .29857000 

.23643444 

.32099980 

.00823214 

WL 

28.00- 

.29857000 

.56173237 

.32777970 

.00655552 

WL 

28.00 

.70142999 

.89278971 

.33433484 

.00655476 

WL 

28.00 

1.70142999 

1.23040181 

.34088923 

.00655400 

WL 

28.00 

2.70142999 

1.57456792 

.34744285 

.00655324 

WL 

28.00 

3.70142999 

1.92528728 

.35399572 

.00655248 

WL 

28.00 

4.70142999 

2.28255912 

.36054783 

.00655172 

WL 

28.00 

5.70142999 

2.64638268 

.36709917 

.00655096 

WL 

28.00 

6.70142999 

3.01675722 

.37364976 

.00655021 

WL 

28.00 

7.70142999 

3.39368197 

.38019959 

.00654945 

WL 

28.00 

8.70142999 

3.77715617 

.38674867 

.00654869 

WL 

28.00 

9.70142999 

4.16717906 

.39329698 

.00654793 

WL 

28.00 

10.20000000 

4.36407894 

.39656149 

.00654755 

WL 

28.00 

11.19999999 

4.76380722 

.40278809 

.00590563 

WL 

28.00 

12.19999999 

5.16944115 

.40837277 

.00526372 

WL 

28.00 

13.19999999 

5.58033880 

.41331553 

.00462180 

WL 

28.00 

14.19999999 

5.99585825 

.41761638 

.00397988 

WL 

28.00 

15.19999999 

6.41535759 

.42127531 

.00333797 

WL 

28.00 

16.19999999 

6.83819491 

.42429233 

.00269605 

WL 

28.00 

17.19999999 

7.26372829 

.42666742 

.00205414 

WL 

28.00 

18.19999999 

7.69131580 

.42840061 

.00141222 

WL 

28.00 

19.19999999 

8.12031554 

.42949187 

.00077030 

WL 

28.00 

20.19999999 

8.55008558 

.42994122 

.00012839 

WL 

28.00 

20.40000000 

8.63607554 

.42995406 

0.00000000 

PEN 

UP 

6000.00000000 

GO 

TO 

2.29857000 

8.63607554 

6.61527891 

-13.45460022 

12.38146087 

-3.85844631 

3.57647199 
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.49019330 

.11412742 

-.82734926 

IDENT.  2 

Y 

FR 

7,50 

0.00000000 

6.61527891 

PEN 

DOWN 

5000.00000000 

FR 

7.50 

2.00000000 

3.96518193 

FR 

7.50 

4.00000000 

2.50103822 

FR 

7.50 

6.00000000 

1.86111844 

FR 

7.50 

8.00000000 

1.68369324 

FR 

7.50 

10.00000000 

1  .66291566 

FR 

7.50 

12.00000000 

1 .71646822 

FR 

7.50 

14.00000000 

1.81791584 

FR 

7.50 

16.00000000 

1.94082342 

FR 

7.50 

18.00000000 

2.06641513 

FR 

7.50 

20.00000000 

2.20655225 

FR 

7.50 

22.00000000 

2.38075530 

FR 

7.50 

24.00000000 

2.60854480 

FR 

7.50 

26.00000000 

2.91122454 

FR 

7.50 

28.00000000 

3.31723123 

FR 

7.50 

30.00000000 

3.85678487 

FR 

7.50 

32.00000000 

4.56010542 

FR 

7.50 

34.00000000 

5.44448553 

FR 

7.50 

36.00000000 

6.47550851 

FR 

7.50 

38.00000000 

7.60583035 

FR 

7.50 

40.00000000 

8.78810703 

PEN 

UP 

6000.00000000 

GO  ' 

TO 

0.00000000 

8.78810703 

IDENT. 

X 

Z 

BK 

4.00 

.36766906 

40.00000000 

PEN 

DOWN 

5000.00000000 

BK 

4.00 

1.36766906 

38.35546875 

BK 

4.00 

2.36766906 

36.87158203 

BK 

4.00 

3.36766906 

35.52490234 

BK 

4.00 

4.36766906 

34.27441406 

BK 

4.00 

5.36766906 

.31738281 

BK 

4.00 

5.36766906 

33.07373046 

BK 

4.00 

6.36766906 

1 .21630859 

BK 

4.00 

6.36766906 

31 .87304687 

BK 

4.00 

7.36766906 

1 .89062500 

BK 

4.00 

7.36766906 

30.62207031 

BK 

4.00 

8.36766906 

2.37829589 

BK 

4.00 

8.36766906 

29.29235839 

BK 

4.00 

9.36766906 

2.69628906 

BK 

4.00 

9.36766906 

27.86230468 

BK 

4.00 

10.20000000 

2.82910156 

BK 

4.00 

10.20000000 

26.57812499 

BK 

4.00 

11.19999999 

2.81616210 

BK 

4.00 

11.19999999 

24.89819335 

FIRST  DER. 

-1.68182502 

-.99841605 

-.49587176 

-.17419213 

-.03337717 

.01039666 

.04095296 

.05829172 

.06241293 

.06480549 

.07695833 

.09887142 

.13054478 

.17465327 

.23387175 

.30820021 

.39763866 

.48279610 

.54428153 

.58209496 

.59623637 


1 

1 

1 

1 

1 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


SECOND  DER. 
.38692065 

.29648831 

.20605598 

.11562364 

.02519131 

.01858253 

.01197376 

.00536499 

.00124378 

.00363634 

.00851648 

.01339661 

.01827674 

.02583174 

.03338673 

.04094172 

.04849672 

.03666071 

.02482471 

.01298870 

.00115270 


1.0. 0-2 


D-82 


t 


BK 

4.00 

12.19999999 

2.61108398 

1 

BK 

4.00 

12.19999999 

23.02709960 

2 

BK 

4.00 

13.19999999 

2.20458984 

1 

BK 

4.00 

13.19999999 

20.85302734 

2 

BK 

4.00 

14.19999999 

1.57788085 

1 

BK 

4.00 

14.19999999 

18.12866210 

2 

BK 

4.00 

15.19999999 

.70214843 

1 

BK 

4.00 

15.19999999 

14.48730468 

2 

BK 

4.00 

16.19999999 

11.06445312 

1 

BK 

4.00 

17.19999999 

9.12500000 

1 

BK 

4.00 

18.19999999 

7.85693359 

1 

BK 

4.00 

19.19999999 

6.72802734 

1 

BK 

4.00 

20.19999999 

5.40087890 

1 

BK 

4.00 

20.40000000 

5.08984375 

1 

PEN 

UP 

6000.00000000 

GO 

TO 

.36766906 

5.08984375 

Section  VI  -  GOBACK  2 
C.  FLOW  DIAGRAM 


REAP 

,K,  NAB 


COKlVf-NiyiQNJ 


E-X.IT  E>oiNT 


REAP X 
S EG-X,  SEG2/  \ 

Point,  pwr.  i,  ) 

\^WR2J^B0W/ 


EMTI5T  rtoiKlj 


♦ 


i. 0.0-2 


D-"4 


1.0. 0-2 


D-3G 


1.0. 0-2 


D-38 


t 


1.0. 0-2 


D-?0 


f 


1.0. 0-2 


D-91 


©~* 


xma*  =  function 

K.5-3 

of  XZ  I N 

FRAME 

EQuatio  n 

552. 


1.0. 0-2 
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*  ACTUAL  FORTRAN  NUMBER 


*  ACTUAL  FORTRAN  NUMBER 


4 


1.0. 0-2 


D-99 


GOBACK  2  -  WITH  PROFILE 


D.  LISTING 


*12C5 

C  CiO  HACK  PROFILE 

0!Ke:.'.!0;  ::{ i 5 ) , a( i 3 , i 7> , z( i  i  )  .c ( i 8)  ,cz( i 6)  f zx( i 5 )  ,zp ( i  1 ) 

ICO  READ  123 

pu::ch  123 

READ  101,  XOR !  S ,  7. OR  I G ,  F  I  HTX,  F I  NTZ ,M,N,K,  KAO 
READ  102,  5E3X,SESZ, POINT, PWR1 ,PWR2,K30W 
PC  I  ;  T* '  PO I  ‘.T-XOR  •  5 ) / F ! KT  X 

MP--1+2 

DO  2C6  ! - 1 ,  ■  \ 

re \d  io2,cz(;) 

20b  COOT  1 1 "JE 
FiVE-ICOO, 

;  i  :=fcooc . 

30fi  :  =  1 ,  mp 

JO  2C2  0-1, i  P 

2C2  r:a:  10  2,  i  ,o) 

3 OR  nO;"|J.!c 

M-MP-3 

N-NP-3 

I)U  203  !-1,M 
READ  102,2(1) 

. F ( M- I ' 230 , 232 , 203 

232  7GTZ=Z( i  * 

203  /. ( I  ' - u  •:  !  •  -ZuR  i  G )  / F :  NTZ 
230  DO  204  ,-1,JJ 

READ  102,:<(i) 

I  F( K-1 ) 340 , 233 , 204 

233  TOTX-X( 1 ' 

204  X( I )=(X; ; )-XOR!G)/FlKTX 

340  I F ( N  >  3 41 ,341 , 3 43 

341  KS-K-3 

DO  342  i=1,RS 
READ  102,/.P(  . ) 

342  ;:p(  i  )-rzp( ;  )-zurig)/fiutz 

DO  TO  203 

343  KS=M 

DO  344  1  =  1 ,11 

344  7.P (  !  )-Z(  1 1 

205  READ  106,1 DKT 
K1-1 


1.0. 0-2 


D-103 


NUMBR-1 

ISTRT-3 

K2-0 

K3-1 

K4-1 

K5-1 

KSR-KS 

KID-1 

XZ-O. 

IND-1 

OLD-O. 

NCOMP-1 

I F( I DNT-1 >600,209,600 

600  I F ( I DNT-2)601 ,310,601 

601  IF(IDNT-3)205,205,234 

C  POINT  ON  THE  PROFILE  CORRESPONDING  TO  W.  L.  ZZ 

690  CALL  FUNCT1  (CZ( 1 ) , CZ (2) , CZ(3) ,CZ (4) , ZZ, FZ) 
FUPZ-CZ(2)+2.*ZZ*CZ(3)+3.*CZ(4)*ZZ*ZZ 
FUSZ-2.*CZ(3)+6.*CZ(4)*ZZ 

I F ( ZZ-ZP ( 1 ) ) 1 79  *  1 79 , 69 1 

691  DO  700  1-1, KS 
DIF-ZZ-ZP(I) 

I F(D I F-. 0000001 ) 1 79, 1 79»  720 
720  FUPZ-FUPZ+3.*CZ( l+4)*D I F*D I F 
FUSZ-FUSZ+6.*CZ( I +4)*D I F 
700  CALL  FUNCT2  ( FZ,CZ( I +4) , 0 1 F, FZ) 

179  GO  TO  (697,405,265,265), I DNT 
697  XX-FZ 
XZ-FZ 
GO  TO  77 

C  EQUATION  OF  W.L.  ZZ 

212  DO  213  1-1, NP 

CALL  FUNCT1  ( A( 1 , I ) , A(2, I ) , A(3 , I ) , A(4, I ) , ZZ, C ( I ) ) 

213  CONTINUE 

I F(ZZ-Z(1 ) )21 6,216,21 5 

215  DO  218  1-1, NP 
DO  219  J-5.HP 
DIF-ZZ-Z(J-4) 

I F(D I F-. 00000001 ) 218, 21 8,82 
82  CALL  FUNCT2  (C ( I ) , A( J , I ) ,D I F,C ( I ) ) 

219  CONTINUE 
218  CONTINUE 

216  I F( I DNT-4) 75,690,690 

75  DO  76  1-1, NP 

76  PUNCH  101,C(I) 

PUNCH  115 

GO  TO  690 

C  POINT  ON  THE  PROFILE  CORRESPONDING  TO  FRAME  XX 
268  GO  TO  (5 71. 5 72, 572) ,NBOW 
572  SEGVA-ZP( 2)-ZP( 1 ) 


1. p.0-2 


D-104 


1 


GO  TO  321 

571  SEGVA=Z(2)-Z(1 ) 

321  I F(XZ-Z(M) ) 581 ,581,580 

580  PAUSE  1 

581  CALL  FUNCT1  (CZ(1 ) ,CZ(2) ,CZ(3) ,CZ(4) ,XZ,Y) 
I F(XZ-ZP( 1 ) ) 3 22, 3 2 2, 300 

300  DO  323  1-1, KS 
DIF-XZ-ZP(I) 

I F(D IF-. 0000001)322, 322, 323 
323  CALL  FUNCT2  (Y,CZ( 1+4) ,DI F,Y) 

322  OIF-XX-Y 

I F(ABSF(D IF)-. 00001)430, 430, 326 

326  GO  TO  (642, 642, 643), K1 

643  I F ( D I F ) 641 ,430,778 
642  I F(D I F) 327, 430, 3 28 

327  K1*»2 
778  OLDZ-XZ 

XZ-XZ+SEGVA 
GO  TO  321 

328  GO  TO  (640,641 ),K1 

640  K1-3 

GO  TO  778 

641  XZ-XZ-SEGVA/2. 

329  CALL  FUNCT1  (CZ(1),CZ(2),CZ(3) ,CZ(4),XZ,Y) 
I F(XZ-ZP( 1 )) 330, 330,331 

331  DO  332  1  =  1  ,KS 
DIF-XZ-ZP(I) 

I F(DIF-. 0000001)330, 330, 332 

332  CALL  FUNCT2  (Y,CZ( 1+4) ,D I F,Y) 

330  D I F=XX-Y 

I F(ABSF(D IF)-. 00001)430, 430, 333 

333  GO  TO  (644, 644, 645), K1 
645  I F ( D 1 F)335 ,430,334 

644  IF(DIF)334, 430,335 

334  SEGVA=XZ-OLDZ 
OLDZ=XZ 
XZ=XZ+5EGVA/2. 

GO  TO  329 

335  SEGVA=XZ-OLDZ 
GO  TO  641 

430  IFODNT-3)  85,292,292 
85  DO  84  1-1, M 
K2-I-1 


84 

798 


292 


IF  (XZ-Z(  0)798,84,84 

CONTINUE 

Y-0. 

RDERI-O. 

SDERI-O. 

GO  TO  270 

GO  TO  (566,566,452) , I STRT 


1.0. 0-2 


D-105 


C  EQUATION  OF  FRAME  XX 
360  DO  313  1-1, MP 

313  CALL  FUNCT1  ( A ( I ,1 ),A( I ,2) ,A( I ,3),A(I ,4),XX,C( I )) 
IF(XX-X(1))902, 902,315 
315  DO  318  1-1, MP 
DO  319  J-5 , NP 
D I F-XX-X( J-4) 

I F(D I F-.00000001 )318,318,83 
83  CALL  FUNCT2  (C( I ) , A( I , J) ,D I F,C ( I ) ) 

319  CONTINUE 
318  CONTINUE 
902  K5-1 

I F( I DNT-4) 78,79,79 

78  DO  80  1-1, MP 
80  PUNCH  101 ,C( I ) 

PUNCH  105 

79  I F(XX-CZ( 1 ) )268, 796, 797 

796  I F ( I DNT-4) 798,566,566 

797  IF(IDNT-4)552,780,780 
780  I  NO-2 

GO  TO  566 

C  *  *  D)  OFFSETS  OF  A  BUTTOCK  *  * 

234  SEG-SEGX/FI NTX 
SGTZ-SEGZ/FINTZ 
ZXZ-TOTZ 
ZZ-Z(M) 

READ  102.BUTTK 
PUNCH  121 

I F(BUTTK) 263,566, 263 
263  DO  293  1-1, N 
293  ZX(I)-X(I) 

TOT-X(N) 

NM-N 

FUPZ-0 

FUSZ-0 

GO  TO  (621,597,597) ,NBOW 
621  XGREA-CZ(l) 

ISTRT-1 
GO  TO  212 

C  INTERSECTION  OF  BTK  WITH  A  WL  OR  FR 

480  Y-Y-BUTTK 

453  IF(ABSF(Y)-. 0001)457, 457, 479 
479  GO  TO  (481 ,482) , K I D 

481  KID-2 
IF(Y)473,457,585 

585  K1-3 

GO  TO  772 

482  GO  TO  (471, 472), K4 
471  GO  TO  (770, 770, 771 ) , K 1 
770  I F(Y)473 ,457,474 


1.0. 0-2 


D-106 


771  IF(Y)474,457,772 

473  K1-2 

772  OLO-XZ 
XZ-XZ+SEGTZ 
GO  TO  275 

474  K4-2 

GO  TO  475 

472  GO  TO  (583, 583, 584), K1 

583  IF(Y)476,457,475 

584  IF(Y)475,457,476 

475  XT-XZ 

XZ-XZ-(ABSF(XZ-0LD))/2. 

OLD-XT 
GO  TO  275 

476  XT-XZ 

XZ-XZ+(ABSF(XZ-0LD))/2. 

OLD-XT 
GO  TO  275 

457  GO  TO  (277,452), I STRT 
452  ZXZ-XZ*F I NTZ+ZOR I G 
273  XZ-XX*F I NTX+XOR i G 
468  CALL  PUNCH3  (BUTTK,XZ,ZXZ,NUMBR) 
NUMBR-NUM8R+1 

466  GO  TO  (650,582,582), I STRT 
650  ISTRT-2 
GO  TO  573 

582  IF(ZXZ-ZMAX)574,573,573 

574  GO  TO  (573,575,575) ,NBOW 

575  I F(BUTTK) 752, 753, 752 

752  IF(BUTTK-XMAX)573,620,576 

753  IF(XX-XMAX)576,620,573 

576  XZ»(ZXZ-ZORIG)/FI NTZ 
SEGVA-ZP(2)-ZP(1 ) 

SEGTZ-SEGVA 

IF(BUTTK)754,755,754 

754  K4-1 

GO  TO  (585, 585, 473), K1 

755  GO  TO  (640, 640, 327), K1 
C  CHECK  FOR  LAST  VALUE 

573  IF(XZ-T0TX)60,64,64 
64  PUNCH  110, SIX 
65  CALL  PUNCH2  (XPLOT.ZXZ) 

GO  TO  205 

60  IF(K3)61 ,460,61 

61  PUNCH  108, FIVE 
K3-0 

460  IF(XX+SEG-X(K2+1))463,461,461 

461  K2-K2+1 
XX«X(K2) 

324  XZ-O. 


1.0. 0-2 


D-107 


NUMBR-1 

IND-1 

C  CHECK  FOR  FR  OUTSIDE  OF  AREA  WHERE  FUNCT.  T  IS  DEFINED 

I F( XX-PO I NT-XGREA)674,675 . 675 
675  NCOMP-2 
GO  TO  673 
674  NCOMP-1 
673  K4-1 
KID-1 
K1-1 

GO  TO  (360,360,266), I STRT 
463  XX-XX+SEG 
GO  TO  324 

C  FUNCT.  T  AND  ITS  1ST  AND  2ND  DERI V 

364  CALL  PR0FL2  (XX, FZ.PWR1 , FUPZ, FUSZ , PWR2.P0 I  NT, T, TP, TS , NCOMP) 
I F(K5-3) 272.590. 750 
272  IF(XZ-TOT)586,586,587 
587  XZ-O. 

GO  TO  452 

586  CALL  FUNCT1 (C(1),C(2),C(3),C(4),XZ,Y) 

I F ( XZ-ZX ( 1 ) ) 368 , 368 , 369 
369  DO  346  1-1, NM 
DIF-XZ-ZX(I) 

I F(D IF-. 000001)368, 368, 346 
346  CALL  FUNCT2  (Y,C( 1+4) ,D I F,Y) 

368  GO  TO  (721, 722), K5 

721  D I FF1-T*D I FF1+Y*TP 
GO  TO  (724, 726), KID 

722  Y=Y*T 

GO  TO  480 

265  GO  TO  (670,364,670), I STRT 
670  XX-FZ 

301  GO  TO  (271,364,269), I STRT 
269  XPLOT-XX*F I NTX+XORI G 
GO  TO  273 

266  I F ( XX-XLARG) 268, 267,267 

267  GO  TO  (275,275,291), ISTRT 
291  ZXZ-ZORIG 

GO  TO  273 
271  XZ-XX 

SEGTZ-SEG 

275  GO  TO  (550,272) ,NCOMP 
550  GO  TO  (553,552,364), ISTRT 
552  ZZ-XZ 

GO  TO  690 
277  DO  294  1-1, M 
294  ZX(I)-Z(I) 

TOT-Z(M) 


1.0. 0-2 


D-108 


♦ 

i) 


NM— M 
XX- XZ 

SEGTZ-ZP(2)-ZP( 1 ) 

XPLOT-XX*F I NTX+XORI G 
77  DO  370  1-1, N 
K2-I-1 

IF(XX-X(l))  81,370,370 
370  CONTINUE 

81  i F ( I D NT-4) 901 ,273,273 

553  XX-XZ 

GO  TO  364 

C  DETERMINE  THE  PT  WHERE  1ST  DERIV.  CHANGES  FROM  -  TO  + 
C  THE  PROFILE  OR  FRAME  EQUATIONS 
566  GO  TO  (596,597,597) ,NBOW 

596  GO  TO  (589,589,598) , I STRT 
598  XMAX-CZ(I) 

ZMAX-ZOR I G 
XLARG-CZ(I) 

GO  TO  567 

597  SEGVA- (ZP(2)-ZP( 1 ) ) / 2 • 

557  GO  TO  (568,568,671), ISTRT 
671  I F(XZ-Z(M) ) 710,710,596 

710  DIFF1-CZ(2)+2.*CZ(3)*XZ+3.*CZ(4)*XZ*XZ 
I F( XZ-ZP( 1 )) 562, 562, 551 
551  DO  563  1-1, KS 
DIF-XZ-ZP(I) 

I F(D I F- .000000 1 ) 562, 562, 563 

563  D I FF1-D I FF1+3.*CZ( l+4)*D I F*D  I  F 
562  GO  TO  (564, 565), KID 

564  GO  TO  (725,725,724), ISTRT 

724  IF(DIFF1)556,555,555 

555  OLD-XZ 
XZ-XZ+SEGVA 
GO  TO  557 

556  K I D=2 

GO  TO  555 

565  GO  TO  (725,725,726), ISTRT 

725  ZZ-XZ 

GO  TO  (779,690), I ND 
779  FZ-XX 
IND-2 
GO  TO  364 

726  IF(ABSF(DIFF1/FINTZ)-.00001)560,560,782 

782  IF(ABSF(OLD-XZ)-.001/FINTZ)560,560,783 

783  I F(D I FF1  ) 55 8 , 560 ,561 

558  GO  TO  (555,559) ,K4 
568  I F(XZ-Z(M) ) 588,588, 589 

588  Dl FF1-C(2)+2,*C(3)*XZ+3.*C(4)*XZ*XZ 
I F(XZ-Z( 1 )) 56 2, 562, 776 
776  DO  777  1-1, M 


i. 0.0-2 
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DIF«XZ-Z( I) 

I F(D I F-.0000001 )562,562, 777 
777  DIFF1«DIFF1+3.*C(I+4)*DIF*DIF 
GO  TO  562 

589  K5-3 
XMAX-C ( 1 ) 

GO  TO  751 

774  K5-3 

CALL  FUNCT1(C(1),C(2),C(3),C(4),XZ,XMAX) 
IF(XZ-Z(1))552,552,591 

591  DO  592  1-1, M 
DIF-XZ-Z(I) 

I F(D IF-. OOOOOOI >552,552,592 

592  CALL  FUNCT2(XMAX,C( 1+4) ,D I F,XHAX) 

GO  TO  552 

590  XMAX-XMAX*T 
K5-4 

751  FZ-CZ(I) 

GO  TO  364 
750  COEF-C( 1 )*T 

I F(XMAX-COEF)594f593,593 

593  XLARG-COEF 
GO  TO  595 

594  XLARG-XMAX 

595  K4-1 
K5“2 
KID-1 
NCOMP-1 

xz-o. 

GO  TO  275 

559  SEGVA-0LD1-XZ 
OLD-XZ 

XZ-XZ+SEGVA/2. 

GO  TO  557 
561  K4-2 

SEGVA-XZ-OLD 

0LD1-XZ 

XZ-XZ-SEGVA/2. 

GO  TO  557 

560  ZMAX-XZ*F I NTZ+ZOR I G 

GO  TO  (774,774,775), ISTRT 

775  CALL  FUNCT1 ( C Z ( 1 ) ,CZ(2),CZ(3),CZ(4),XZ, XMAX) 
„  I  F(:<Z-ZP(  1 ))  577,5  77, 569 

569  DO  570  1-1, KS 
DIF-XZ-ZP(I) 

I F(D IF-. 0000001)577, 577, 570 

570  CALL  FUNCT2  ( XMAX,CZ( I +4) ,D I F, XMAX) 

577  I F ( XMAX— CZ ( 1 ) ) 5 78 , 5 78 , 5 79 

578  XLARG-CZ(I) 

GO  TO  567 


1.0. 0-2 


D-110 


579  XLARG-XMAX 
567  K4-1 
K5-2 
KID-1 

XGREA-XLARG 

IF( I DMT-4) 795 • 794, 795 
795  XZ-O. 

GO  TO  360 

794  IF(BUTTK)622,690,622 

622  GO  TO  (212,690,623) , ISTRT 

623  ISTRT-1 
GO  TO  212 

620  ZXZ-ZMAX 
GO  TO  273 

C  *  *  A)  OFFSETS  OF  A  WATER  LINE  *  * 

209  READ  102, ZT 
FINT-FINTX 
SEG-SEGX/FI NTX 
RORIG-XORIG 
TOT-TOTX 
DO  236  1-1 ,N 
236  ZX(I)-X(I) 

ZZ-( ZT-ZOR I G ) / F I NTZ 
GO  TO  212 

405  GO  TO  (901,70), NCOMP 

C  FUNCT.  T  AND  ITS  1ST  AND  2ND  DERI V 

901  CALL  PR0FL2  (XX, FZ, PWR1 , FUPZ, FUSZ.PWR2, PO I  NT, T, TP , TS , NCOMP) 
70  Y-C(1)+C(2)*XZ+C(3)*XZ*XZ+C(4)*XZ**3 

TDERI-C(2)+C(3)*2.*XZ+C(4)*3.*XZ*XZ 
SDER I =C (3)*2.+C (4)*6.*XZ 
I F(K2-1 )406 ,408 ,408 

408  DO  409  1-1, K2 
D I F»XZ-ZX( I ) 

I F(ABSF(D IF)-. 00000001) 409, 409, 252 
252  PRODCT-C ( l+4)*D I F 
Y*=Y+PRQDCT*D  I  F*D  I  F 
TDERI-TDERI+3.*PR0DCT*DI F 
SDERI-SDERI +6.*PR0DCT 
I F( I -K2) 409, 406, 406 

409  CONTINUE 

406  RDER l«(TDER I *T+Y*TP) /FI  NT 
SDERI«(SDERI*T+2.*TDERI*TP+Y*TS)/(FINT*FINT) 

Y«Y*T 

270  ZXZ-XZ*F I NT+ROR I G 

CALL  PUNCH1  (ZT,ZXZ,Y,RDERI , SDER I , IDNT) 

412  IF(ZXZ-TOT)497,51,51 
51  PUNCH  110, SIX 

CALL  PUNCH2  (XPLOT,Y) 

GO  TO  205 


1.0. 0-2 
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497  I F(K3)52.415 ,52 
52  PUNCH  108, FIVE 
XPLOT-ZXZ 
K3-0 

415  I F(XZ+SEG-ZX(K2+1 )) 416,501 ,501 
501  K2-K2+1 

XZ«ZX(K2) 

GO  TO  432 

416  XZ-XZ  +SEG 

432  GO  TO  (903,799), I DNT 

799  IF(XX-P0INT-XGREA)552,405,4O5 

903  XX- XZ 

GO  TO  405 

C  *  *  B)  OFFSETS  OF  A  FRAME  OR  STATION  *  * 

310  READ  102, ZT 
FINT-FINTZ 
SEG-SEGZ/FINTZ 
RORIG-ZORIG 
TOT-TOTZ 
DO  401  1-1, M 
401  ZX(I)-Z(I) 

XX-(ZT-XORIG)/FI NTX 
GO  TO  566 

101  F0RMAT(4F15 .8,41 5 ) 

102  FORMAT(5F15.8, 15) 

104  FORMAT (2F 15. 8, 15) 

105  FORMAT (7H  IDENT. ,9X, 1HZ, 1 7X, 1HY, 1 1 X, 10HFI RST  DER. ,4X, 1 1HSEC0ND 
DER 

C.) 

106  F0RMAT(I4,F15.8) 

107  F0RMAT(2HWL,F7.2,F14.8.3X,3F15.8) 

108  F0RMAT(8HPEN  D0WN,F15.8) 

109  F0RMAT(2HFR,F7.2,F14.8,3X,3F15.8) 

110  F0RMAT(6HPEN  UP,2X,F15.8) 

111  F0RMAT(2HLG, I3,3X,F15.8,3X,F15.8) 

112  F0RMAT(7H  I  DENT. ,9X, 1HZ , 1 7X, 1HY, 1 2X,9HT0LERANCE) 

113  F0RMAT(7H  I  DENT. ,9X, 1HX, 1 7X, 1HY. 1 2X,9HT0LERANCE) 

114  FORMAT (2HDP, F6.2, FI 5.8,3X, 2F15 .8) 

115  FORMAT (7H  IDENT. ,9X,1HX,17X,1HY,11X,10HFIRST  DER. ,4X, 1 1HSEC0ND 
DER 

C.) 

116  FORMAT (1 X,6HIDENT.,9X,1HX,17X,1HY,7X,1 4HF I RST  DERIVAT.) 

117  FORMAT ( 1  X,6H I  DENT. , 10X, 1HX, 1 7X, 1HY, l4X, 1HZ) 

118  F0RMAT(1X,6HIDENT.,9X,1HZ,17X,1HY,7X,i4hFIRST  DERIVAT.) 

120  F0RMAT(7X, 1HA, 13X, 1HB, 13X, 1HC, 13X, 1HD,9X,5HSTART,8H  FINISH) 

121  FORMAT (7H  I  DENT. , 10X, 1HX, 1 7X, 1  HZ) 

122  FORMAT (4F1 4.8, 2Fo. 3) 

123  FORMAT (50H  ) 

124  FORMAT(25H  ) 
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125 

*1205 

) 


2 


4 


7 

6 

1 

5 

3 

*1205 


1 

2 

5 

7 

8 

6 

9 

10 
3 


FORMAT ( 8X, 1HX, 14X, 1HZ , 14X, 1HY) 
END 


SUBROUTINE  PR0FL2(XD,FD,EXPD1 ,TZPD,TZSD,EXPD2,X1 ,TD,TPD,TSD,NCD 

TD-( X1+FD-XD) /XI 
IF(TD-. 00001)1,1 ,2 
DIFD-EXPD1-1 . 

PROD-TD**DIFD 

TPD1-EXPD1*TZPD*PR0D 

TSD1-EXPD1*(D I FD*TD**(D I FD-1 . )*TZPD**2+TZSD*PR0D) 
TD-1.-TD**EXPD1 
I F(TD-. 00001 ) 7,7,4 
DIFD-EXP02-1 . 

TDEXP«TD**DIFD 
TPD— EXPD2*TPD1*TDEXP 

TSD-EXPD2*(DIFD*TD**(EXPD2-2.)*TPD1*TPD1-TSD1*TDEXP) 

TD-TD**EXPD2 

IF(TD-.00001)7,7,6 

TD-O. 

GO  TO  3 
NCD-1 
RETURN 
NCD-2 
TD-1  . 

TPD-O. 

TSD-O. 

RETURN 
NCD-1 
GO  TO  5 
END 


SUBROUTINE  PUNCH1 

XD-XID 

YD-Y  I  D 

I  F(  XD ) 1 ,2,2 

XD— XD 


( ZD ,  X I  D ,  Y  I  D ,  RD  f  SD ,  I 


I F ( YD ) 3 » 4, 4 
IF(YD)5,6,6 
YD— YD 


GO  TO  (7,8), ID 

PUNCH  20, ZD, XD,YD,RD,SD 

RETURN 

PUNCH  21 , ZD,XD,YD,RD,SD 
RETURN 

GO  TO  (9, 10), ID 

PUNCH  2 2 , ZD , XD , YD , RD , SD 

RETURN 

PUNCH  23 , ZD , XD , YD , RD , SD 
RETURN 


D) 


YD— YD 
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GO  TO  (11, 12),  ID 

11  PUNCH  24 , ZD , XD , YD , RD , SO 
RETURN 

12  PUNCH  25»ZD,XD,YD,RD,SD 
RETURN 

4  GO  TO  (13, 14), ID 

13  PUNCH  26 , ZD, XD, YD ,RD, SD 
RETURN 

14  PUNCH  2 7 , ZD , XD , YD , RD , SD 
RETURN 

20  FORMAT(2HWL,F7.2,F14.8,4X,1H-,F13.8,2F15.8) 

21  F0RMAT(2HFR,F7.2,F14.8,4X,1H-,F13.8,2F15.8) 

22  FORMAT(2HWL,F7.2,F14.8,3X,3F15.8) 

23  F0RMAT(2HFR,F7.2,F14.8,3X.3F15.8) 

24  FORMAT(2HWL,F7.2,1H-,F13.8,4X,1H-,F13.8,2F15.8) 

25  FORMAT(2HFR,F7.2,1H-,F13.8,4X,1H-,F13.8,2F15.8) 

26  FORMAT(2HWL,F7.2,1H-,F13.8,3X,3F15.8) 

27  FORMAT(2HFR,F7.2,1H-,F13.8,3X,3F15.8) 

END 


*1205 

SUBROUTINE  PUNCH2  (XID.YID) 

XD-XID 
YD-YID 
I F( XD) 1 , 2 , 2 

1  XD— XD 

I F(YD)3,4,4 

2  I F(YD)5 ,6,6 

3  YD— YD 

PUNCH  30,XD,YD 
RETURN 

4  PUNCH  31 , XD.YD 
RETURN 

5  YD— YD 

PUNCH  32, XD.YD 
RETURN 

6  PUNCH  33, XD.YD 
RETURN 

30  FORMAT  (5HG0  T0,4X,1H-,F13.8,4X,1H-,F13.6) 

31  FORMAT  (5HG0  T0,4X,1H-,F13.8,3X,F15.8) 

32  FORMAT  (5HG0  TO ,4X,F1 4.8.4X, 1 H-.F13.8) 

33  FORMAT  (5HG0  T0,4X,F14.8,3X,F15.8) 

END 


*1205 

SUBROUTINE  PUNCH3  (BUTTK, XZ.ZXZ.NUMBR) 
IF  (XZ)  1,2,2 

1  XZ4— XZ 

PUNCH  126, BUTTK, XZ4.ZXZ.NUMBR 
RETURN 

2  PUNCH  119, BUTTK, XZ.ZXZ.NUMBR 
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RETURN 

119  F0RMAT(2H8K,F6.2,F15.8,3X,F15.8,6X,I2) 

126  FORMAT ( 2H8K , F6 , 2 , 2H  -,F13 .8,3X, F15.8.6X, 1 2) 

END 

*1205 

SUBROUTINE  FUNCT1 (Cl ,C2,C3,C4,XD,YD) 

YD-C1+C2*XD+C3*XD*XD+C4*XD*XD*XD 

RETURN 

END 

*1205 

SUBROUTINE  FUNCT2  (YD1 ,C5 ,DI FD,YD2) 

YD2-YD1 +C5*D I FD*D I FD*D I FD 

RETURN 

END 

*1205 

SUBROUTINE  DERI VI  (C2D,C3D,C4D,FXD,FZPD,FZSD) 

FZPD«C2D+2.*C3D*FXD+3.*C4D*FXD*FXD 

FZSD-2.*C3D+6.*C4D*FXD 

RETURN 

END 

*1205 

SUBROUTINE  DERIV2  (C5D,DI FD,FPD,FSD,FZPD,FZSD) 

PROD-3. *C5D*D I FD 

FZPD-FPD+PROD*D I FD 

FZSD-FSD+2.*PR0D 

RETURN 

END 
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GOBACK  3 


A.  OPERATING  INSTRUCTIONS 


This  version  of  GOBACK  will  accept  coefficients  of  two-dimensional 
curves,  either  single  or  double  splined,  as  shown  in  Fig.  D-l. 

From  these  coefficients  it  will  do  the  following: 

(1)  Calculate  offsets 

(2)  Calculate  first  and  second  derivatives 

(3)  Plot  a  picture  of  the  curve 


The  program  is  written  in  FORTRAN  for  the  IBM-1620  computer.  In 
order  to  plot  the  results,  one  of  the  subroutines  of  Appendix  F 
must  be  included  in  the  FORTRAN  compiler  and  the  hardware  requirements 
of  Appendix  F  must  be  met. 


Fortran  Symbol  Definitions 

Symbol  definition 

RORIG  -  The  x  coordinate  of  the  first  point  on  the  full 

size  curve 

SEG  -  Interval  between  calculated  offsets  on  the  full  size 

curve 

CA  -  Value  of  the  offset  of  the  first  point  on  the 

full  size  curve 

N  -  Number  of  coefficients  in  the  equation  without 

the  first  coefficient 

L  -  Number  of  offsets  given  in  the  data,  except  the 

first  offset 

K  -  Number  of  specific  points  where  offsets  must  be 

calculated  (  0  if  none  required) 
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cx 

CY 

FINT 

A 

XA 

XE 


Scale  factor  for  plotting  in  the  x  direction 
on  the  curve  (lengthwise  on  the  plotter),  x  plotted 
equals  x  calculated  times  CX 

Scale  factor  for  plotting  in  the  y  direction  on 
the  curve  (crosswise  on  the  plotter),  y  plotted 
equals  y  calculated  times  CY 

The  interval  between  the  first  two  stations  on 
the  original  curve  (x^-Xq  ) 

The  coefficients  of  the  equation 

The  x  coordinates  of  the  original  data  offsets 
(stations)  from  Station  2  to  the  last  station 

List  of  x  coordinates  of  specific  points  where 
offsets  are  to  be  calculated  in  order  from  least  to 
greatest  x  value 


Input  Data  Cards 

The  description  of  the  cards  given  below  is  in  terms  of  the  actual  FORTRAN 
format.  The  standard  F  field  description  -  for  fixed  point  decimal 
numbers,  and  the  I  field  -  for  integer  numbers  which  must  always  be  right 
justified,  are  used.  The  fields  given  come  consecutively  across  the  card. 

Card  Formats 

Contents  Card  No. 


Format 

F15.8 

F15.8 

F15.8 

Variable 

RORIG 

SEG 

CA 

1 

Format 

15 

15 

15 

Variable 

N 

L 

K 

2 

Format 

F15.8 

F15.8 

F15.8 

Variable 

CX 

CY 

FINT 

3 

Format 

F15.8 

Variable 

A 

Next  N  cards 

Format 

F15.8 

Variable 

XA 

Next  L  cards 

Format 

F15.8 

Variable 

XE 

Next  K  cards 
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Output  Data 


The  first  output  will  be  the  input  data  punched  in  the  same  format  as 
it  was  read  in,  except  cards  3  and  2  are  punched  on  one  card.  Follow¬ 
ing  that,  a  header  card  will  be  punched  and  then  the  x  coordinate, 
y  coordinate,  and  derivatives  of  the  curve  as  required.  The  curve  may 
be  plotted  simultaneously  with  this  output,  or  the  output  may  be  omitted 
and  only  the  plot  obtained. 


Sense  Switch  Settings 


Switch 

1 

OFF  - 

Both  the  first  and  second  derivatives  of  the 
curve  will  be  punched 

Switch 

1 

ON 

Only  the  first  derivative  will  be  punched 

Switch 

2 

OFF  - 

The  curve  will  not  be  plotted 

Switch 

2 

ON 

The  curve  will  be  plotted 

Switch 

3 

OFF  - 

All  output  will  be  punched 

Switch 

3 

ON  - 

No  output  will  be  punched 

All  the  switch  positions  may  be  changed  any  time  during  the  execution 
of  the  program  to  obtain  partial  data. 
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B.  SAMPLE  PROBLEM  FOR  GOBACK  3 


Input  Data 


0.0 

11  9  0 

0.05 

0.0 

0.0 

0.0000000 
-0.2068768 
-0.8883194 
1 .7682645 
-1 .9088810 
2.0039195 
-1 .6010168 
2.7202405 
-0.25504729 

25.5 

38.25 
51.0 
63.75 

76.5 

89.25 
102.0 
114.75 
127.5 


5.0  19.50 


x  coordinates  of 
original  data  offsets 


Output 


.00000000  5.00000000  19.50000000 

.05000000  .12500000  25.50000000  11  9 

1  .00000000 

2  .00000000 

3  .00000000 

4  -.20687680 

5  -.88831940 

6  1  .76826450 

7  -1 .90888100 

8  2.00391950 

9  -1.60101680 

10  2.72024050 

11  -.25504720 

1  25.50000000 

2  38.25000000 

3  51.00000000 

4  63.75000000 

5  76.50000000 

(Continued  on  next  page) 
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6 

7 

8 
9 


89 . 25000000 

102.00000000 

114.75000000 

127.50000000 

X 

.00000000 

5.00000000 

9.99999990 

15.00000000 

20.00000000 

25.00000000 

25.50000000 

29.99999800 
35.00000000 
38.25000000 

39.99999900 

44.99999800 
50.00000000 
51.00000000 

54.99999900 

59.99999800 

63.75000000 
65 .00000000 

69.99999900 

74.99999800 

76.50000000 
80.00000000 

84.99999900 

89.25000000 

89.99999800 

94.99999800 

99.99999900 
102.00000000 
105.00000000 
110.00000000 

114.75000000 
115.00000000 
1 20 .00000000 
125.00000000 

127.50000000 


Y 

19.50000000 
19.50000000 
19.50000000 
19.50000000 
19.50000000 
19.50000000 
19.50000000 
19.49886400 
19.48930300 
19.47414100 
19.46167700 
19.39101300 
19.22961200 
19.18208500 
18.93476200 
18.51419200 
18.13450500 
17.99801600 
17.38924400 
16.63805100 
16.37656200 
15.69371800 
14.55902700 
13.47049500 
13.26826900 
11 .83795200 
10.23781400 
9.54177140 
8.43461480 
6.46138250 
4.50586050 
4.40217670 
2.34013220 
.35073487 
-.59364910 


I  RST  DERIVAT. 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
-.00075794 

-.00337800 
-.00608461 
-.00836172 
-.02155531 
-.04465637 
-.05046548 
-.07254610 
-.09466698 
-.10726189 
-.11123872 
-.13413327 
-.16820734 
- . 1 8060960 
-.20901960 
-.24369906 
-.26771338 
-.27159293 
-.30178948 
-.33952080 
-.35672312 
-.38035745 
-.40608937 
-.41472047 
-.41475082 
-.40760490 
-.38569304 
-.36919978 
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Section  VII  -  GOBACK  3 


C.  FLOW  DIAGRAM 


99 
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GOBACK  3 
D.  LISTING 


C  *  *  2D  GO  BACK  -  *  * 

DIMENSION  XA( 100) ,A( 100) ,B( 100) , XZ( 100) , XE( 100) ,XEZ( 100) 
99  READ  2,R0RIG,SEG,CA 
READ  7,N,L,K 
READ  2,CX,CY,FI NT 
PRINT  42 
PAUSE 

IF  (SENSE  SWITCH  3)  360,361 
361  PUNCH  2,  RORIG,SEG,CA 
PUNCH24,CX,CY,FINT,N,L,K 
360  X-0. 

KID-1 

I F(SENSE  SWITCH  2)701,702 
702  DUMMY«PLOT( 11111.) 

701  XX-O. 

IND-4 

OLDX-ROR I G 

Y-0. 

DO  22  1-1, N 
READ  2 , A ( I ) 

IF  (SENSE  SWITCH  3)22,36 
36  PUNCH  26, I ,A( I ) 

22  CONTINUE 
DO  33  1-1, L 
READ  2 , XA( I ) 

IF  (SENSE  SWITCH  3)  410,39 
39  PUNCH  26, I , XA( I ) 

410  XZ ( I ) — ( XA ( I )-RORIG)/FINT 
33  CONTINUE 
L-L+3 

I F ( K ) 109* 109 , 1 10 
110  DO  104  1-1, K 
READ  2, XE( I ) 

XEZ ( I ) — ( XE ( I )-ROR I G) /F I  NT 
I F(SENSE  SWITCH  3)104,112 
112  PUNCH  26, I , XE ( I ) 

104  CONTINUE 

109  I F(SENSE  SWITCH  3)  55,400 
400  I F( SENSE  SWITCH  1)183,184 

183  PUNCH  44 
GO  TO  55 

184  PUNCH  43 

55  Y Y-C A+A ( 1 ) *  X  X+A ( 2 ) * ( X  X**  2 ) + A ( 3 ) * ( XX**3 ) 

TDER l-A( 1 )+A(2)*2,*XX+A(3 )*3 .*XX*XX 
SDER l-A( 2)*2.+A(3)*6.*XX 
I F( I ND-5 ) 100 ,  66,77 
66  YY-YY+((XX-XZ(1))**3)*A(4) 

TDER I -TDER 1+3. *A( 4)* ( (XX-XZ( 1 ) )**2) 
SDERI-SDERI+6,*A(4)*(XX-XZ(1 ) ) 
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GO  TO  100 

77  I  N*»  I  ND-1 

503  1)0  78  1-4,  IN 

TDERI=T0ERi+(3.*A( I )*((XX-XZ( l-3))**2)) 
SDERI»SDERI+(6.*A(I )*(XX-XZ( 1-3))) 

78  YY-YY+((XX-XZ( l-3))**3)*A(  I ) 

100  XX«XX*F I NT+ROR I G 

YR-YY 

XR«XX 

RDERI»TDERI / F I  NT 
SOERI=SDERI /(FI NT*F I  NT) 

203  iFf SENSE  SWITCH  3)81,700 
700  iF( SENSE  SWITCH  1)185,186 

185  PUNCH  40,  XR,YR,RDERI 
GO  TO  81 

186  PUNCH  40,  XR,YR,RDERI , SD ER I 

81  ! F( SENSE  SWITCH  2)201,103 

201  I F( SENSE  SWITCH  3)202,1 

202  PRINT  42 
PAUSE 

GO  TO  203 
103  yR=Y";*CY 
XR-XR*CX 
DUMMY-PLOT(YR) 

1  l f(XX-XA(L-3))88,99,99 
88  IF((0LDX+SEG)-XA(IND-3))87,97J97 

197  5 F C K ) 30 1 , 30 1  ,311 

311  I  F<K-K 10)301 ,402,402 

40  2  I F ( XA ( I ND-3 ) -XE ( K I D ) ) 30 1  , 309 , 303 

303  XX=XEZ(KIDl 

i  F ( K-K 10)55,308, 308 

308  K I D=K 1 0+1 
GO  TO  55 

309  i  F '  L- 1 ND ) 303 , 303 , 30  2 
302  I ND= I ND+1 

GO  TO  303 
301  XX*XZ(lM0-3.' 

I  F(L-i  ?;D)55 ,55,310 

310  IND-I.YD+1 
GO  TO  55 

97  !F(K)304.304,312 

312  .F(K-Ki 0)304,403,403 

403  » F ( :<A ( !  '•  J-3 ) -XE  (K 1 0) ) 304, 305 , 303 
305  OLDX-OLDX+SEG 
GO  TO  309 

304  OLDX=OLDX+SE G 
GO  TO  301 

87  IF(K)307,307,313 

313  I F(K-K i 0)307,401 ,401 
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40 1  I F(Ol.OX+SEG-XE(K I D) ) 307,306 , 303 

306  0LDX=ULDX+5EG 
30  TO  303 

307  XX« (OLDX+SEG-ROR I  G)  /F I  NT 
QLDX-OLDX+SEG 

50  TO  55 

2  F0AMA7(4F15 .8) 

7  FORMAT ( 15, 15, 15) 

24  R>  iK.vr(3F15.8, 15, 15, 15) 

26  FORMAT  (1 6, Ft 5. 8) 

40  F0„r,AT(8;<,F15.8,3X,3F15.8) 

42  f'0'lFAT(29HSET  SWITCHES,  THEN  PUSH  START) 

43  FORMAT '16X,  1 HX,  1  7'X,  1HY,7X,1 4HF IRST  DER I  VAT. ,  1  X,  15HSEC0.J 

44  FORMAT ( 16X,1HX,17X,1HY,7X,14HFIRST  DER I  VAT.) 

END 


DERI  VAT 
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PROGRAM  FOR  CONVERSION  OF  STANDARD 
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PROGRAM  FOR  CONVERSION  OF  STANDARD 
CUBIC  COEFFICIENTS  FOR  NUMERICAL  CONTROL 

OPERATING  INSTRUCTIONS 

This  program  accepts  coefficients  of  cubic  equations  of  the  standard  form: 

Y  «  A  +  BX  +  CX2  +  DX3 

Xx  <x  <  x2 

with  the  following  conditions  : 

(1)  Y  Units  are  feet 

(2)  X  units  are  feet 

(3)  X^  and  X2  units  are  feet 

(4)  X,  X^  ,  X2  may  be  scaled 

(5)  The  coefficients  are  in  the  normal  fixed  point  format  of  FORTRAN, 
and  are  on  cards  as  punched  by  GO  BACK  (See  Fig.  E-l) 

These  coefficients  are  converted  to  meet  the  following  conditions 

(1)  Y  units  in  inches 

(2)  X  units  in  inches 

(3)  X^  and  X2  in  inches 

(4)  X  ,  X^  and  X2  rescaled  -  usually  to  full  size 

(5)  The  coefficients  are  presented  in  normalized,  excess  50 
floating  point  format  similar  to  that  used  in  SPS-I  for  the 
1620  computer . 

(6)  The  origin  of  the  equation  and  its  limits  may  have  been 
translated  along  the  X  axis. 

(7)  The  coefficients  and  limits  are  punched  on  a  card  along  with 
the  appropriate  AUTOMAP  statement  as  shown  in  Fig.  E-l,  and  in 
an  order  such  that  the  equation  now  reads : 

Y  -  DX3  +  DX2  +  BX  +  A 

As  may  be  noted  from  the  flow  chart,  this  program  contains  a  fairly  concise 
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TYPICAL  INPUT  CARDS 


4. 54152120  1.0953 1033 

I  I 


>.453  32  740 

.41337164 

1 

1 

:  0.000 

1.333 

1 

1 

.19681363  -. 02331354  1.933  3.333 

I  I  II 


0  0  0  0  0  0  0  oooo oo|oooooo|oooo|o|ooooaooooooooooo  0  0  0  0|  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

I  I  1  4  !  I  7  I  I  II II  77  II 14  II II  I?  «  II BI1  77  IS  74  75  B  77  74  71  B  >1 11  B  II B  B II  »  B 4141  U« 41 4S  41 47  ■  41 II 51  17  11  S4  B  SI 57 B  SI  « II  U U 44 B  «I7 II  B  7|  71  71  71 74  IS  71  77  7|  71  B 

iiiiim|ii|iiiiiijiim|iiiimii|iit|tiiiiiiiitiii|inm|i!iimimiiiiinii 

2222222222|2|2222222222222222222222222222222Z2222|22222222222222222222222222Z22? 

0  3  3  3  )|  3  3  3  3  3  3  3  3  3  3  3  3  3|  3  3  3  3  3  3| 3  3  3  3  3  3|  3  3  3  3  3|3|  3  3  3  1  l|  0  0  3  3  3  Ill  3  3  3  3|  3  3  3  3  3  3||  1  3  3  3  )  3  3  3  3  3  3 
4444|44|44444444 4444 4 44444444444444 4444444444 4 444444444|44 4444444 4444444444 444 44 
9bb6bb|bb|bb8bbbbbbbbb|bb5bbbbbbbbbbbbbbb5bbbbbbbbbbbb|bbbbbbb5bbbbbbbbbbbbbbbbb 

0010066006066(6000(600600166806010101 00 01060000000006060606600600010608086600000  S 

J  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  J  7  7  7  1  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  1  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  1  1  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 

(l888|B8M88(((((l(|80  6|6B6((((0(|((l|mS81S8l|((||l|(((S((|((((((l|((8(((((0(B 
99S999999999999999999|9999999999939|99999999S9S99999999999999|||3S999|||99999999 

t  2  3  4  5  I  3  •  I  II  11  12  13  U  IS  15  17  t«  1|  ft  }1  2?  73  24  25  »  27  21  29  *  )1  U  33  M  »  ft  J7  ft  ft  «  <1  42  43  4«  45  41  (7  #«  ft  51  52  59  M  ft  M  57  ft  ft  ft  II  12  ft  14  ft  ft  17  ft  ft  TO  71  72  7J  74  75  7|  77  7«  71 10 


TYPICAL  OUTPUT  CARDS 


M 64030323  4333-93*33 

u _ 

-4720012373  4916401135  5110353103  5037346010 


23.  3s  47. 3i 

I  I 


■ 

oooooo ooooo ooo oo||o oo ooo oo oqo|qooooooo|oooo|ooo|ooooo|o|ooooooo|oo  0  0  0  0  0|  0  0  0  0  0  0  0  0 

I  1  1  I  S  •  I  ■  ■  II  li  II  II 74  IS  II  17  II  IS  a  II  11  71  H  13  B  77  B  B  S  II  B  B  M  8  B  17  B  B  •  II  II  IS  II  4!  II  41 II  «  SI  SI  SI  SJ  51  B  »  51  B  H  10 II  H  I]  14  B  M  17  n  M  70  It  71  II  II  »  71 77  II  71  B 

1  1  11  1  1  1  1  1  1  1  1  1  1  11  I  1 1  1  1  11  1  1  1  |  1  I  1||  11  1  1  ||  1  1  1  t |  1  1  1  1  I  1  1  1  1  1  1 | 1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 
222222222222222|222|222222222222222222222222222222222222222|22222222222222222|2? 

0  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3|  3  3  3  3  3  3  3  3  333|3I3|(33  3  3  3  3  3  3  3  3||  3  1  3  3  3  3  3  1  3|  3  3  3  3  3  3  3  3 
444444444444 4|4444444444| 44 4| 444444444444444444444 4|444444444444444|444444444444 
5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  b  5|  5  5  5  5  5  5  5  5  5  5|  5|  5  5  5  5|  5  5  5  b  5|  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5|  5  5  5 
60006606(6606006006666600001006000601006066086006000106680600(600000608606068060  * 
77777777777777|T77777|777777777777T77777777777777|777777777777777777|77777777777 
I  (((((( 6  (0(1(1110 (li| (i 10 8S((Bi(ll(lll(((|l(((((((|((((((((((|(| 0(0  I (|(| I (0(01 0  0 
9999909999999999000999999I 9990999909 80 S| 999009990000099909990 9| 0(0099 9| 909999990 

I  1  1  I  S  I  I  I  I  l|  II  II  I]  14  IS  II  II  II  l|  B  H  11 11 »  IS  B  17  71  B  B  II  S  S  14  B  B  17  8  B  «  II  II  41 14  B  ««  41  •  SI  SI  S  B  S4  It  »  17  B  B  M II  Cl  IS  II  «  a  II  8  M  II II  71  71  II  IS  71  77  71  II  8 


Fig.  E-l  Typical  Data  Cards 


l.O. 0-2 


routine  for  converting  from  fixed  point  to  excess  50  normalized  floating 
point  in  FORTRAN,  as  well  as  overcoming  FORTRAN  format  difficulties  in 
punching  the  resulting  numbers. 

A  number  of  sets  of  coefficients,  which  together  describe  a  given  curve, 
may  be  processed  at  the  same  time.  The  AUTOMAP  program  requires  that  the 
first  card  of  a  particular  set  be  punched  with  a  different  format  than 
that  of  the  others.  The  program  does  this  automatically  and  also  punches 
a  curve  identification  number  in  each  card  of  each  set. 

The  program  is  written  in  FORTRAN  II  for  the  IBM-1620  computer.  The 
variable  word  length  feature  of  this  version  of  FORTRAN  has  been  used. 

Integer  numbers  have  been  increased  in  length  to  ten  digits .  The  length 
of  floating  point  numbers  remains  at  eight  digits. 


A.  Input  Data 


Symbol 

Format  Card  Columns 

Definition 

1st 

Card 

N 

110 

1-10 

No.  of  segments 

T 

F10.1 

11-20 

Translation  Factor 

SF 

F10.5 

21-30 

Scale  Factor* 

2nd 

and  Remaining 

Cards 

rk 

1-14 

A,B,C,D 

4F14.8  ^ 

B 

lC 

15-28 

29-42 

Cubic  Coefficients 

i 

'd 

43-56 

XB 

F8.1 

57-64 

Beginning  Point 

XE 

F8.1 

65-72 

Ending  Point 

The  X  values  of  the  limits  and  assumed  in  the  equation  are  multiplied 
by  the  scale  factor.  Thus,  the  scale  factor  for  doubling  the  curve 
length  would  be  2 
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b.  Output  Data 


Each  Card 


Symbol  Format  Card  Columns 


I 

'  D 

13-23 

D.B.C.A 

1 

4111  H 

lc 

24-34 

)B 

35-45 

1 

'a 

46-56 

XB 

F7.1 

57-63 

XE 

F7.1 

65-71 

ISTA 

15 

73-77 

Definition 

Converted  Coefficients 

Beginning  Point 
Ending  Point 
Curve  Identification 
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SAMPLE  PROBLEM  -  CONVERSION  OF  STANDARD  CUBIC  COEFFICIENTS 
FOR  NUMERICAL  CONTROL 


Input  Data 

N  T  SF 


r 

:> 

0.0  1.0C00C 

,\ 

i*» 

o 

r 

D 

5TA.IT 

Fit:  0 

r> 

: .00030000 

6.40809230 

-2.45932740 

.41387163 

0.000 

1.99 

9 

Ci 

4.^41 521 20 

1 .09581030 

.19681363 

-.02881854 

1.999 

3.99 

9 

n 

2.: 70 12650 

2.94933630 

-.26657288 

.00979699 

3.999 

7  ■  99 

9 

6.39770760 

1.32631330 

-.06371750 

.00134469 

7.999 

12.00 

6.69046240 

1.25332460 

-.05761844 

.00117527 

12.000 

15.999 

Output 

f- 

1  j 

D  C 

3 

A 

X3 

XE 

r 

r  U  i — 

'  -CU3/ 

4828741 085-5020494395 

5164080923 

4983333333 

0.0, 

23.9 

52 

-4720012875  4916401135 

5110958103 

5037846010 

23.9, 

47. 

9, 

52 

4668034652-4922214406 

5129493563 

5017251054 

47.9, 

95. 

9, 

52 

4593381250-4853097916 

5113265133 

5053314230 

95.9, 

144.0 

i 

52 

4581615972-4848015366 

5112533246 

5055753853 

144.0, 

191 . 

9, 

52 
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PROGRAM  FLOW  CHART  -  CONVERSION  OF  STANDARD  CUBIC  COEFFICIENTS 

FOR  NUMERICAL  CONTROL 
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SUBPROGRAM  CU8ICCA,&;C,D,T,Sr  j 


A.  A / 12 

a.  s/sp 

c  -  c/  (_ii  ■»  SP*) 

D  m  O  j  (1AM*  SF*] 


T2.  T 
T3.T* 


As  A  -  B*T  +  C*T2  _  B*TJ 
S»B_2*C*T  +  1ABPT1 
C,C  _3*D*T 


RETURN  TO 
CALLING  PROGRAM 


CONVERT  COEFFICIENTS 
FOR  X  ANt>  Y  \N  INCHES 
AN6  FOR  NEW  SCALE 


CONVERT  THE  COEFFICIENTS 
FOR  TRANSLATION  ALONG 
THE  X-AXIS 
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PROGRAM  LISTING  -  CONVERSION  OF  STANDARD  CUBIC  COEFFICIENTS 
FOR  NUMERICAL  CONTROL 


*0810 

1  READ  1 50 , N, T, SF 
T=T*1 2 . 

ACCEPT  149,1  STA 
DO  21=1 , N 

HEAD  1J1 ,A,3,C,D,XB,XE 

XB*(XH*12.*SF>+T 

HE* ( XE* 1 2 . *5F } +T 

CALL  CUB  I C ( A, 3 , C , D , T , S F) 

J1-0 

J2-0 

v,3=0 

J4-0 

CALL  FHT(A,J4) 

CALL  FMT (3  ,  J3 ) 

CALL  FMT ( C , J  2 ) 

CALL  FMT (D , J 1 ) 

!F(I-1)3,3,4 

3  PUNCH  1 52 , J 1 , J2,J3,J4,XB,XE, ISTA 
GO  TO  2 

4  PUNCH  153, J1,J2,J3,J4,XB,XE, ISTA 

2  CONTINUE 
PAUSE 

GO  TO  1 

1 49  FORMAT ( 15) 

153  FORMAT ( 1 2 X , 4 1  1 1 ,F7.1 , 1H, f F7. 1 , 1H, , 1 X, I  5) 

152  FORMAT ( 1 2HHULL  =CUB/ ,4 1 1 1 , F7. 1 , 1H, , F7. 1 , 1 H , , 1 X , I  5 ) 

151  FORMAT (4F1 4.8, 2F8. 1 ) 

150  FORMAT ( I 10.F10.1 ,F10.5) 

END 

*0810 

SUBROUT .  r.C  CUB  I  C  ( A,B ,  C  ,D ,  T ,  SF) 

A-A/12. 

B=B/SF 

C=C/(12.*SF*SF) 

D=D/ ( 1 44.*SF*SF*SF) 

T2=T*T 

T3=T2*T 

A=A-B*T+C*T2-D*T3 

B-B-2.*C*T+3.*D*T2 

C=C-3.*D*T 

RETURN 

END 

*0810 

SUBROUTINE  FMT(CF.J) 

AF=CF 

CF=ABSF(CF) 

J=50 

1  IF(CF-10.)3,5,2 

2  J=J+1 
CF=CF/ 10 . 
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GO  TO  1 

3  I F ( C F— 1 .)4,6,6 

4  J-J-1 
CF«CF*10. 

GO  TO  3 

5  J-J+2 

CF-CF* lOOOOOO. 

GO  TO  7 

6  J-J+1 

CF-CF*10000000. 

7  NF-CF 
IF(AF)8,9,9 

8  J-(J*100000000+NF)*(-1) 
GO  TO  10 

9  J«J*100000000+NF 
10  RETURN 

END 
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Appendix  P 


PLOTTING  BOUT INKS 


CMKBHXf 

lection  Pege 

jjfHOBDCTiOK . r-i 

i  run  smuornat . r-3 

n  msxKDcnoMS  to  us*  or  the  pLormc  no gun  .  .  r-7 

IHPUI  BAZA . r-7 

srasx  much  skttiugs.  .  . . r-9 

PHOCKBUXI . r-9 

SAMnX  PROBLEM . P-10 

ni  un  or  ihk  subkouiikk  a  a  romux  pkockam  .  .  .  r-11 

BAXPU  PX0BIEM . P-12 

IT  rwi  BIAGXAM  FOR  7URUC  PKOCKAM . P-U 

▼  nit  mo  m  run  pkogxam . r-21 

n  Lisnna  rat  rasur  1  maomn . r-33 

tii  ujhhb  rat  pokoum  iz  uuounm . r-47 
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PLOTTING  ROUTINES 


INTRODUCTION 

When  working  with  mathematical  equations  for  curves,  particularly  if 
there  are  certain  requirements  on  the  characteristics  of  the  curves,  it 
is  desirable  to  have  some  type  of  automatic  plotting  device.  This  will 
provide  an  accurate  picture  of  the  contour  of  the  curve  and  permit  a 
subjective  analysis  of  the  curve's  characteristics. 

Tne  plotting  device  used  in  developing  this  system  consists  of  a  small 
digital  incremental  plotter  (CalComp  Model  560-R)  which  is  connected 
directly  to  the  IBM-1620  computer.  By  means  of  a  switch,  output  signals 
normally  routed  to  the  paper  tape  punch  of  the  computer  can  be  channeled 
directly  to  the  plotter  instead.  By  feeding  the  plotter  a  series  of 
digits  (0,  ...  ,9)  which  it  recognizes  as  instructions  (See  Section  I 

of  this  Appendix),  line  increments  can  be  drawn. 

The  purpose  of  the  plotting  program  is  to  feed  the  plotter  the  series 
of  digits  which  will  cause  the  plotter  to  draw  a  line  between  any  two 
given  points.  If  enough  points  to  describe  the  periphery  of  an  object 
are  made  available  to  the  plotting  program  in  order,  it  will  cause  the 
plotter  to  draw  a  picture  of  that  object. 

Plotting  routines  which  accept  data  points  from  two  different  sources 
are  included  with  this  report:  first,  the  plotting  program,  which  reads 
data  directly  from  cards  processed  through  the  card  reader  (see  Section 
II  of  this  Appendix  for  operating  instructions);  and  second,  the  plotting 
subroutine  which  is  used  with  a  master  FORTRAN  program. 

The  master  FORTRAN  program  collects  or  calculates  the  desired  data  points, 
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and  places  them  in  certain  computer  memory  locations  It  then  calls  on 
the  plotting  subroutine  which  plots  the  calculated  points. 

There  are  two  versions  of  this  subroutine,  one  is  meant  to  be  placed  in 
a  FORTRAN  I  subroutine  deck  (see  the  program  listing  in  Section  V  of  this 
Appendix);  the  other  is  meant  for  the  FORTRAN  II  subroutine  deck  (see 
program  listing  in  Section  VI)  The  use  and  operating  instructions  of 
each  of  these  programs  are  the  same  (Section  III),  only  the  program 
instructions  themselves  change 

All  of  the  plotting  programs  are  written  in  SPS-1  for  the  IBM-1620  computer 
The  subroutine  programs  are  intended  for  a  machine  with  60-K  digits  of 
memory  The  plotting  program  can  be  used  with  any  size  memory. 
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PLOT  SUBROUTINE 


METHOD 

The  plotter  traces  straight  line  increments  along  one  of  the  eight 
directions  (four  axes)  shown  in  Fig.  1.  The  increments  are  .01"  in  the 
X  and  Y  directions,  and  .01414"  for  the  inclined  directions. 

A  punch  tape  instruction,  1,  2,  3,  ...,  8,  will  move  the  plotter  in 
one  of  eight  specific  directions.  Instruction  0  will  raise  the  pen, 
and  9  will  lower  the  pen.  ,  Y 


Fig.  1  Direction  of  Movement  of  Plotter 


Given  a  particular  quadrant  (Fig.  2),  and  a  specific  point,  the  line 
connecting  this  point  to  the  origin  must  fall  in  one  of  the  five  cases 
listed  below: 


Case  1  Point  P  :  AY  ■  0 

Case  2  Point  Q  :  AX>  Ay 

Case  3  Point  R  :  AX  -  AY 

Case  4  Point  S  :  Ay  >  A  X 

Case  5  Point  T  :  AX  =  0 

Case  5 


X  Case  1 


Fig.  2 

After  determining  the  case,  the  problem  is  reduced  to  positioning  the 
point  to  one  side  of  the  bisector  of  the  quadrant. 
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For  Point  P  having  a  null  or  less  than  .005  increment  in  the  Y  direction, 
the  procedure  is  as  follows: 

1.  Round  off  the  increment  at  the  1/100  position  by  adding  or 
subtracting  .005  to  its  value  if  positive  or  negative  respectively. 

2.  Add  this  rounded  value  to  the  corresponding  coordinate  of  Point  0 
and  store  it  as  the  origin  to  determine  the  increment  for  the  next 
point  to  be  joined. 

3.  Check  for  digit  in  the  units  position  of  the  absolute  value  of  the 
increment  for  fast  one-inch  moves. 

4.  Check  for  digit  in  the  1/10  position  for  fast  l/10th-inch  moves. 

5.  Test  for  digit  in  the  1/100  position  for  single  moves. 

By  subtracting  one  from  the  corresponding  position  after  each  move,  the 
increment  will  be  zero  at  the  completion  of  the  total  move. 

For  Point  R  the  procedure  is  the  same  as  mentioned  above  for  Point  P 
with  the  difference  that  both  increments  X  and  Y  ...,  (Fig.  2)  will 
be  rounded  off  and  added  to  the  coordinates  of  Point  0  . 


The  method  for  plotting  those  Points  Q  having  unequal  increments  in 
X  and  Y  is  as  follows: 

1.  Get  the  ratio  of  the  smaller  to  the  larger  increment  in  absolute 
value . 

2.  Round  off  the  increments  X  and  Y  and  add  them  to  the  coordinates 
of  the  Point  0  to  use  these  new  values  as  the  origin  to  calculate 
the  increments  for  the  next  point  to  be  joined. 

3.  Approximate  the  line  0Q  (Fig.  3)  with  segments  0A,  AB,  etc.,  in 
the  following  manner. 

4.  Set  an  XY  counter  equal  to  .005 

5.  Add  MN  (Ratio  x  .01")  to  the  counter 

6.  Compare  the  counter  with  AM  (=  .01)  (check  for  digit  in  the 
1/100  position) 

a.  If  the  counter  is  greater,  it  means  that  we  are  closer 

to  the  actual  line  going  through  the  diagonal  rather  than 
the  horizontal  line  __  __ 

MN  +  . 005 >  AM 

AM  -  MN  <.005 

(Note :  The  XY  counter  has  to  be  reduced  by  MA  after  the 
diagonal  move  and  the  cycle  goes  back  to  Step  5.) 
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b  If  the  counter  is  less  than  MA  Che  closest  path  is  the 
horizontal  line 

MN  +  .005  <  AM 
AM  -  MN  >  .005 

(Note;  The  cycle  goes  back  directly  to  Step  5  after 
the  horizontal  move.) 


An  analysis  of  the  moves  from  0  to  B  of  Fig  3  will  help  to 
clarify  the  method 


1. 

XY 

counter  =  MN  4 

005  >  AM  _ 

2. 

Move  0A 

3. 

XY 

Counter  =  MN  + 

.005 

-  AM 

4, 

XY 

Counter  =  MN  + 

.005 

-  AM  +  MN 

=  2MN  + 

.005 

-  Si 

5. 

XY 

Counter  =  LK  + 

005 

-  AM  <  Si 

/.LK  + 

. 005  <  2AM 

005  <  HK  -  LK 


6 .  Move  AB 


K-r, 
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INSTRUCTIONS  FOR  USE  OF  THE  PLOTTING  PROGRAM 


INPUT  DATA 

There  are  three  types  of  input  cards  for  the  plotting  program: 

(1)  Scaling  data 

(2)  Pen  raising  and  lowering  instructions 

(3)  Coordinates  of  data  points 

(1)  Scaling  Card 

If  it  is  desired  to  plot  the  data  at  some  scale  other  than  1"  ■  1  unit 
dimension,  a  scale  factor  must  be  entered  with  a  value  for  scaling 
X  and  a  value  for  scaling  Y.  The  plotted  values  are  determined  by: 

X  *  X  (*)  scale  factor 
Y  *  Y  (*)  scale  factor 
The  format  for  this  card  is  given  below: 


Card  columns 

Contents 

1  -  9 

Blank 

10 

Minus  (-)  for  negative  X  scale  factor 
Blank  for  positive 

11 

-  19 

X  scale  factor 

15 

Decimal  point  for  X  scale  factor 

20 

-  27 

Blank 

28 

Minus  (-)  for  negative  Y  scale  factor 

29 

-  37 

Y  scale  factor 

33 

Decimal  point  for  Y  scale  factor 

38 

-  80 

Any  alphanumeric  information 

If  scaling  is  required,  this  card  must  be  the  first  valid  card  of 
the  data  deck. 

Only  cards  with  a  decimal  point  in  Column  15  are  recognized  by  this  program. 
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(2)  Pen  Raising  or  Lowering  Instruction  Cards 


Cards  must  be  placed  In  the  deck  to  Inform  the  program  when  the 
plotter  pen  must  be  raised  or  lowered.  These  cards  have  the  follow¬ 
ing  format.  Any  card  which  does  not  have  a  decimal  point  in 
Column  15  is  ignored. 

Pen  Down  Card 


Card  Columns 

1-10 

U  , 

12  ,  13  , 

14  ,  15 

16  -  38 

39  -  80 

Contents 

Blank 

5 

0  0 

0 

Blank 

Anything 

Pen  Ud  Card 

Card  Columns 

1  -  10 

H  , 

12  ,  13  , 

14  ,  15 

16  -  38 

39  -  80 

Contents 

Blank 

6 

0  0 

0 

Blank 

Anything 

(3)  Data  Coordinate  Cards 


Coordinates  of  points  to  be  plotted  are  presented  one  set  to  a  card 
in  the  following  format.  The  origin  of  the  picture  will  always  be 
the  first  point  in  the  data  set. 


Card  Columns 
1  -  9 


10 

11 

-  19 

15 

20 

-  27 

28 

29 

-  37 

33 

38 

-  80 

Contents 

Blank  or  any  alphanumeric  information 
(Ignored  by  the  plotting  program) 

Minus  (-)  if  the  X  coordinate  is 
negative,  blank  if  plus 

X  coordinate  of  point  to  be  plotted 

Decimal  point  of  X  coordinate 

Same  as  Cols.  1-9 

Minus  (-)  if  the  Y  coordinate  is 
negative,  blank  if  plus 

Y  coordinate  of  point 

Decimal  point  of  Y  coordinate 

Same  as  Cols .1-9 
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SENSE  SWITCH  SETTINGS 

Switch  1  ON  -  All  values  of  X  and  Y  will  be  scaled  by 

multiplying  them  by  the  scale  factors  given  on 
the  first  card 

Switch  1  OFT  -  Values  will  be  plotted  "full  scale,"  that  is,  each 
inch  on  the  plot  will  be  equal  to  1  unit  of  the 
coordinate  dimensions.  No  scale  card  used. 

PROCEDURE 

The  following  sequence  of  events  takes  place  when  using  this  program 

(1)  Load  program  deck  and  data  deck  In  card  reader,  including  a 
scaling  card  if  desired 

(2)  Push  load  button  on  card  reader,  causing  the  program  to  load 

(3)  Program  loading  is  complete  when  the  typewriter  types  the 
message,  "Switch  1  On  to  Scale" 

(4)  Set  the  switch  as  desired,  pushing  "start"  will  execute  the 
program 
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SAMPLE  PROBLEM  FOR  PLOTTING  PROGRAM 


SAMPLE  DATA  FOR  PLOTTING  PROGRAM 


.5  .5 

2.  2. 

6.  2. 

6.  6. 

2.  6. 

2.  2. 

3.  3. 

6000.  0. 

5.  5. 

5000.  0. 

6.  6. 

6000.  0. 

6.  2. 

5000.  0. 

5.  3. 

6000.  0. 

2.  6. 

5000.  0. 

3.  5. 

5.  5. 

5.  3. 

3.  3. 

3.  5. 


FIGURE  PLOTTED  FROM  SAMPLE  DATA 


1.0. 0-2 


F-10 


USE  OF  THE  SUBROUTINE  IN  A  FORTRAN  PROGRAM* 


The  calling  sequence  for  the  plot  subroutine  is: 


X  =  Any  expression 

Dummy  =  PLOTF  (Y) 

Where : 

X  =  Coordinate  X  of  the  point  to  be  plotted.  This  value 
must  be  on  the  right-hand  side  of  the  instruction 
preceding  the  plot  instruction.  This  is  to  place  this 
value  in  FAC  (Floating  Point  Accumulator). 

Y  =  Coordinate  Y  of  the  same  point 
Dummy  =  Any  floating  point  dummy  variable  not  used  in  computation 


a.  Dummy  «=  PLOTF  (11111.)  is  an  instruction  that  must  be  placed  at 
the  beginning  of  the  program.  This  will  set  the  subroutine  to 
start,  put  the  pen  down,  and  ready  the  subroutine  for  a  set  of 
data. 

b.  Dummy  =  PLOTF  (50000.)  will  only  put  the  pen  down 

c.  Dummy  =  PLOTF  (60000.)  will  put  the  pen  up 


Instructions  b.,  and  c.,  are  complete  by  themselves  and  do  not  require 
any  prior  instruction. 

Attention  is  called  to  the  fact  that  the  origin  of  the  plotted  figure  is 
automatically  set  as  the  first  point. 

This  subroutine  can  only  handle  X  and  Y  values  having  up  to  three 
integers  (at  the  moment  they  are  stored  in  FAC  and  argument,  respectively, 
ready  to  be  plotted). 

*The  User  should  refer  to  the  FORTRAN  Manual  for  the  procedure  to  incor¬ 
porate  this  subroutine  in  the  library  subroutines  deck. 


*  *  *  SAMPLE  FORTRAN  PROGRAM  *  *  * 

c 

SAMPLE  PROBLEM  FOR  PLOT 
JIJL 1 0=0 

SUBROUTINE 

n 

1C 

SET  THE  PLOT  SUBROUTINE 
rUMMYePLOTF (11111.) 
f  E  AO  100 .SCALX, SCALY 

TO  INITIATE 

r 

READ  IN  X  Ai'KJ  Y  COORDINATES  .  F  THE  POINT 

11 

RFAD  100.X. Y 

j-N 

CHECK  FOR-  PEN  UP  OR  DOWN 

IF(X-r>0C0.)12,13,l4 

r- 

13 

PUT  THE  PEN  DOWN  • 
D1JKMY-PL0TF ( £0000 , ) 

GO  TC  11 

* 

14 

PUT  THE  PEN  UP 
DUMMYwPLOTF ( 60000 . ) 

GC  TO  11 

12 

Y-Y*SC V  v 

X=X*$C.’.LX 

DIIMMYnPLOTr  (Y) 

GO  TO  11 

HO 

FO=tMA7(8X,F15.8,3X,F15. 

FI!D 

8) 

*  *  *  SAMPLE  DATA 

*  *  * 

.5 

.5 

2. 

2„ 

6. 

2. 

6. 

6. 

2. 

6, 

2. 

3. 

3. 

600 C , 

C, 

r>. 

5. 

50C0. 

r> 

A 

6. 

6. 

ficco . 

t . 

’*•'  • 

• 

:xco. 

5. 

3. 

vC'C’O , 

t 

9 

2. 

i:* 

5000 , 

u. 

3. 

5. 

5. 

5. 

5. 

3. 

3. 

3. 

3. 

5. 

l.o. n- 

F-12 

Section  IV  FLOW  DIAGRAM  FOR  PLOTTING  PROGRAM 
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FLOW  CHART  OF  PLOT  SUBROUTINE 
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Section  I 


INTRODUCTION 


This  report  gives  e  detailed  description  of  the  input/output,  operating 
instructions,  and  computing  procedures  (flow  diagram)  of  a  shell 
development  program. 

The  program  was  written  to  provide  a  more  accurate  and  economical 
method  of  developing  shell  plate  than  the  graphical  method  now  employed 
by  the  shipbuilding  industry. 

The  development  is  made  from  a  table  of  final  lofted  offsets  which 
locate  the  plate  on  the  hull  surface.  With  these  offsets,  the  program 
computes  a  set  of  two-dimensional  coordinates  which  represent  points  on 
the  perimeter  of  the  plate  at  each  frame.  These  points,  when  plotted 
and  faired,  can  be  used  as  a  pattern  for  cutting  the  shell  plate. 

The  program  is  written  in  FORTRAN  for  an  IBM- 1620  computer. 


Section  II 


THE  PROBLEM 


The  process  of  developing  shell  plete  seperetes  Into  two  distinct 
portions:  — - 

(1)  Calculating  the  true  lengths  of  line  segments  on  the 
hull  surface 

(2)  Assembling  these  calculated  line  segments  into  the  shape 
of  the  plate  to  be  cut. 

The  information  necessary  to  solve  the  problem  is  a  complete  table  of 
offsets,  including  those  at  the  plate  sight  edges. 

FINDING  THE  TRUE  LENGTH  OF  THE  LINE  SEGMENTS 

Since  the  equation  of  the  hull  surface  is  unknown,  the  equation  for 
each  line  segment  is  also  unknown.  The  curve  must,  therefore  be 
approximated  by  passing  a  known  function  through  the  data  points.  The 
length  of  the  known  curve  in  that  interval  may  then  be  found. 

The  method  which  is  presented  is  independent  of  this  approximating 
function.  Because  the  equation  for  passing  a  circle  through  the  given 
points  is  simple  and  convenient  to  calculate,  circles  were  used  for  this 
purpose.  Other  functions,  such  as  parabolas,  could  be  easily  substituted. 

Any  three  points  in  three-dimensional  space  are  co-planner.  For  this 
reason  it  is  most  convenient  to  fit  the  approximating  function  to  groups 
of  three  offsets  at  a  time. 

This  equation  is  most  accurate  if  the  plane  passing  through  the  three 
points  is  very  nearly  perpendicular  to  the  center  plane  of  the  ship, 
i.e.,  the  projection  on  the  centerline  plane  of  a  line  passing  through 
the  three  points  is  nearly  straight.  This  condition  is  met  by  points 
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along  frame  lines  and  along  plate  sight  edges. 

In  finding  the  length  of  lines  such  as  Line  7  of  Fig.  1  (Page  G-12  of 
this  Appendix),  this  condition  is  only  rarely  met.  Usually  the  coordi¬ 
nates  of  the  line  of  Points  ...  B .  are  much  closer  to  one  edge  of 
tlie  plate  than  to  the  other  edge.  In  this  situation  the  plane  passing 
through  the  points  necessary  to  find  the  arc  length  (Points  Cq  ,  , 

A^  for  Line  7)  intersects  the  centerplane  at  some  angle  very  different 
than  90°.  A  circle  passing  through  these  three  points  would  project  as 
an  arc  on  the  centerplane.  The  length  of  the  line  obtained  would  be  in 
error  by  the  amount  the  projected  arc  exceeds  the  length  of  projected 
straight  lines  between  the  points. 

For  this  reason,  the  shell  plate  for  finding  the  length  of  these  lines 
is  assumed  to  be  a  cylinder,  and  the  arc  length  found  using  the  method 
of  Section  IV  of  this  Appendix. 

Establishing  the  Flat  Plate  Contour 

Once  the  true  lengths  of  the  curves  have  been  found,  they  must  be 
assembled  to  accurately  determine  the  shape  of  the  flat  plate.  This 
procedure  is  accomplished  by  triangulation  (See  Fig.  1,  Page  G-12). 

Since  the  lengths  of  the  three  sides  of  each  of  the  triangles  are  known, 
the  triangles  are  completely  determined.  The  procedure  then  becomes 
that  of  assembling  these  triangles  and  solving  for  the  coordinates  of 
points  on  the  periphery  of  the  plate  which  determine  its  shape. 

An  outline  of  this  procedure  is  given  in  Section  V. 
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Section  III 


MX3ATIOM8  FOR  PASSIMS  A  CIRCLE  THBOOGB  3  POIMTS 

Length  of  the  circular  ere  constructed  through  three  erbltrery  points 


Objective:  Determine  the  length  of  the  ere  A  given  only  the  coordinates 
of  the  three  points  ,  P^  »  *3  • 

S  ,  Sj  ,  can  be  calculated  from  the  well-known  formula  for  the 
distance  between  two  points,  l.e., 

sx  -  +  <Y3'Y1>2  +  ^3’Zl>2 

S2  -  I/O^V  +  (I2-y3)2  +  <Z2‘Z2)2 

s3  -  +  (r2“ij)2  +  (Zj-zp2 

For  the  present,  assume  R  (Radius  of  Circle  through  Pj  ,  ?2  ,  P^) 

is  known. 

L.  0.0-2 


6-5 


The  problem  now  becomes :  Determine  A  as  a  function  of  S  end  R 
only.  The  length  of  en  ere  on  e  circle  is  equel  to  the  product  of 
the  redius  and  the  central  angle  determined  by  the  end  points  of 
this  arc.  In  the  immediate  application  this  implies: 

(1)  A  «  RT 

T  must  now  be  expressed  in  terms  of  S3  and  R. 

Employing  the  law  of  cosines  with  respect  to  triangle 
P^  Q  Pj  ,  it  can  readily  be  seen  that: 

S3  -  2R2  -  2R2  cos  (T) 

or 

S3 

— I  -  1  -  cos  (T) 

2RZ 

Using  the  identity  Sin2(T/2)  -  1/2  -  1/2  cos  (T) 

above  expression  becomes 

Sin2  T/2  -  S2/4R2 

or 

Sin  T/2  -  S3/2R 
Equivalently 

(2)  T  -  2  Sin"1  S3/2R 

Combining  (i)  and  (2)  the  expression  for  A  in  terms 
of  S3  and  R  is  obtained: 

(3)  A  -  2R  Sin'1  ~ 
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Using  a  trigonometric  series: 

.3 


r' 


(U)  sin-1(x)  -  X  +  £  +  LU?  +  +  + 

w  sin  w  a -r  6  ^2.4.5  t2-4-6-7  2-4*6-8‘9  2*4*6-8*10*ll 

The  argument  X  must  now  be  expressed  in  terms  of  and  R. 
Utilizing  law  of  sines  with  respect  to  U  P^  P ^  P3  • 

S, 


11 


SinU 


■  Diameter  of  circumscribed  circle 


77-r  *  ZR 
SinU 


(5) 


1  -  2 SinU 

R  S, 


Using  law  of  cosines  with  respect  to  U  P3 

S3  *  S2  +  S1  ‘  2S1  S3  C°8  U 

2  2  2 
Cos  U  -  Sj  +  S*  -  S* 

2S1  S2 


(6) 


also 


Sin2U  -  l-Cos3U 


Sin' 


>  (S2+S1“S3) 

■u  -  1  -?  1  „  J  y 


2  2 
4S2!  S22 


Combining  (5)  and  (6) 


1  Ah8!5; 

1 


*»i:  j  -  (s2 4 1 


e2  2 
S2  S3 
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Combining  (3),  (4),  and  (7),  tha  final  expreaalon  for  arc 
length  becoaea : 


Replacing  by  3^  In  (8)  and  (9)  would  give  the  arc  length 
between  and  ?3  ,  etc. 
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SECTION  IV 


DIAGONAL  ARC  LENGTH  APPROXIMATION 


The  following  describee  the  method  used  to  find  the  arc  lengths  of 
certain  lines  (See  Fig.  I,  lines  3,  7 ,  10,  11  of  Section  V)< 


Inown:  Pi  *  P2  *  P3  *  P5 

Xo  *  X1  -  X3 


£  *  j  FJicV-xf+terZ.t ;  H  =  \y3-Y,U(1,'Z^ 


0)  6  =  V( Yr-Yot+^-g'f  (Z)  DzkYi-Yit+CZ,-^ 

Knowing  lengths  of  K,F,  and  H  ,  arc  lengths  Pj  P^  and  P^Pj  can  be 
computed  using  (1)  and  (2)  with  the  method  described  in  Section  III. 


Then  using  the  Pythagorean  theorem: 


rp>  T^TW^f 
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Section  V 

'  *  ‘  ,‘r'- 

OUTLINE  OF  THE  ANALYTIC  PZOC1DURE  FOR  ESTABLISHING 
THE  FLAT  PLATE  CONTOUR 


In  71{ i  X f  Line  >  . » .  .  C ^ .  end  Line  i  f  .  • .  .  A^ 

ere  bn  the  upper  end  lower  eight  edges  of  the  plate;  Line  Bq, 

■j .  is  either  e  waterline  or  buttock  line  which  passes 

reasonably  close  to  the  mean  line  between  the  sight  edges. 

The  following  assumes  the  coordinates  of  Points  A  ,  ...  ,  A.; 

o  i 

B  ,  ...  ,  B. ;  C  .....  C.  are  known  on  the  ship's  hull 
o  l  o  1 

(heights  and  half-breadths) ,  and  also  assunes  that  the  arc 
lengths  between  the  points  have  been  calculated  using  Appendix 
I -A,  or  some  similar  method. 

Bq  Is  the  plate  origin  (See  Fig.  1)  and  Its  coordinates  are 
therefore  known.  Line  BqB^  Is  assumed  to  be  parallel  to  the 
Z  axis 


B  B,  ,  1  X~  ,  i.k  are  known 
0  1  O  O  1  o 


So  "  2  <Vi +  +  *7? 


ro-  2 


(So  -  J7h)  <8o  -  TO  <So  •  B 


??]' 


cos  I  •  (Mi)  +  (BA)  -  (O  ) 
I  0  '  o  1' _ ^  o  0 _ '  1  o' 


*  <*?? 


o y 


■  2r  s 


ol 


Coordinates  of  Point  A 


ox 


-  o 


O  0 


cos 


Po 
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BoCb  •  *1^0  to0,m 
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“  1  - - - - "S - • - 
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Vi  *  BiBi  known 


S2“s2  <¥o  +  X?1  +  W 

*2  "  [?2  <S2  -  ¥o>  <S2  "  Vl>  <S2  “  Vl> 


'  tan_1 


v 

A  Y 
o 


°<2  -  90°  -  to  "  y2  “  /Sq  »  where  AoX  U  n***tlv* 
°^2  ■  90°  -  £  +6  -  )T,  ,  where  A  X  ie  positive 

0  '  O  L  O 

or  soro 


i 

t 
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tnli-*  Mi-* 


V  ■  V  +  *b*i  *in  *s 


A.x  ■  A  x  +  A  A,  cos 
1  0  0  1 

C  C,  >  C,  B  Known 
oi  11 


•S 


Coordinates  of  A^ 


O^o  +  CoCl  *  C1B1> 


[<s3  ■  <s3  *  W  (S3  ’  W] ' 


-  tan 


-1 


“  Arctan 


2 r 


S3  ’  *1*1 


1  -  r. 


S3  -  BjCj 


•S-  rh 

♦  *3 

+ 

90°  ,  where  c  X 
o 

is  negative 

•s-  f! 

•  ^ 

+ 

y3  -  90°  ,  where  CqX 

Is  positive 

or  zero 

V  •  v 

+  ToCl 

sin 

V| 

C.x  ■  C  x 
1  0 

+  ^*1 

cos 

*3 

/ 

)  Coordinates  of 

C1 

A1B2  Known 


*1*2  *  ’ 


LOOP 


Alcl - 

[<C1X  " 

X4  * 

,  where 

H  “ 

„  -1  fClY 

tan  \cxx 

V 

180°  -  tan 

:iy  •  v>2  ] 


where  A^  < 


■1  N  ’  XA 

Vaxx  -  cpj  . 


where  k^X  >■  CjX 


S4  ”  2  <A1C1  +  A1B2  +  B2C1) 

*4  •  |4  <S4  '  W  <S4  -  *?2>  ‘S4  ’  W  ] 

P ^  »  ArcCan 


*u-  (f-M  -  pu  -x.  - 


'5  "  ’2  4'  r4 

B2y  "  A^y  +  AjBj  «in  Y5 


B2*  ■  A^x  +  A^B2  pos  V5 


4  "  P4 


Coordinate*  of  B. 
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Coordinates  of  C. 
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-  CjX  +  CjC2  cos 


REPEAT  FROM  LOOP  FOR  NEXT  TRIANGLE,  ETC. 


,,,  |-JMf 
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Section  VI 


SHIFT  OF  SIGHT  EDGE 


In  manual  methods  of  shell  plate  development,  the  neutral  plane  is 
graphically  projected  to  define  the  equivalent  flat  plate  pattern. 

The  neutral  plane  is  determined,  conventionally,  by  a  point  1/2  t  from 
and  normal  to  the  molded  line  of  the  plate,  where  t  ■  plate  thickness. 
This  is  necessary  since  the  finished  mold  loft  offsets  lie  on  the  mold 
line,  which  is  normally  the  interior  surface  of  the  plate.  The  point 
which  is  the  intersection  of  the  plate  edge  and  the  mold  line  is  called 
the  sight  edge  offset,  as  shown  in  Fig.  G-2. 

To  exactly  duplicate  the  foregoing  in  a  computer  program,  calculation 
of  the  slope  of  the  curve  would  be  necessary  to  establish  the  normal 
to  the  plate  through  the  sight  edge  offset.  This  complication  can  be 
avoided  entirely  by  permitting  the  final  plate  sight  edge  to  deviate 
slightly  from  that  provided  by  the  loft.  Of  importance  is  that  the 
deviation  in  sight  edge  on  one  plate  is  entirely  compensated  for  on  the 
adjoining  plate.  The  plate  itself  remains  in  exactly  the  same  position 
on  the  hull.  Fig.  G-2  shows  the  method  employed  in  the  computer 
program  wherein  the  neutral  plane  is  approximated  by  finding  the  point 
which  is  from  the  sight  edge  offset  along  the  plane  of  measurement. 
Whether  the  projection  is  to  a  vertical  plane  or  a  horizontal  plane 
depends  upon  whether  the  sight  edge  offset  is  defined  by  the  mold  loft 
from  waterline  half  breadths  or  buttock  heights,  respectively. 

It  is  obvious  that  no  shift  of  sight  edge  will  occur  where  the  slope 
of  the  mold  line  is  either  zero  or  infinite.  For  all  other  slopes,  the 
deviation  varies  as  the  cosine  of  the  angle  of  tangency,  with  maximum 
deviation  occurring  when  that  angle  equals  forty-five  degrees. 
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In  a  one-half  Inch  shell  plate,  the  maximum  deviation  would  be  slightly 
less  than  one-fourth  inch.  Since  such  a  negligible  shift  in  sight  edge 
permits  a  substantial  simplification  in  the  computer  program,  its  adop¬ 
tion  is  considered  justifiable. 

The  foregoing  is  presented  in  detail,  since  the  computer  method  produces 
a  result  which  deviates  from  conventional  lofting  methods.  This  deviation 
must  be  known  by  the  loftsman,  otherwise  it  will  appear  that  the  computer 
solution  produces  a  measurable  error  in  some  cases,  when  checked  in  the 
loft  by  conventional  methods.  Further,  if  both  computer  development  and 
manual  plate  development  are  done  on  the  same  ship,  the  manual  method 
must  conform  to  the  computer  method  in  order  to  achieve  compatibility. 
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Section  VZI 


INPUT  AND  OUTPUT 


When  shell  plating  Is  placed  on  the  ship  it  is  butted  or  joined  between 
frames.  Because  the  offsets  are  provided  on  the  frames  only,  the 
computing  methods  used  in  the  program  require  that  the  offsets  be  given 
one  frame  ahead  of  the  butt  at  the  start  of  the  plate,  and  for  two  frames 
beyond  the  butt  at  the  other  end.  This,  then,  means  that  for  a  plate 
covering  twenty  frames,  offsets  for  twenty-three  frames  are  needed. 

The  offsets  provided  by  the  loft  are  given  for  the  upper  and  lower  sight 
edges,  as  well  as  for  a  waterline  or  butt  line  on  the  plate  which  passes 
between  the  sight  edges.  The  program  provides  for  a  maximum  plate  length 
of  fifty  frames  and  a  minimum  of  two  frames. 

The  offsets  are  given  in  feet,  inches,  eighths,  and  sixteenths,  and  are 
in  the  form  X  ,  Y  ,  and  Z  ,  where  X  is  frame  spacing,  Y  is  half¬ 
breadths,  and  Z  is  height.  When  the  plate  is  on  the  bottom  of  the  hull, 
the  7-  (detailed  above)  value  must  be  subtracted  from  the  Z  terms. 

*  t 

However,  the  program  was  written  so  that  ^  is  algebraically  added  to 
the  Y  terms;  consequently,  the  X  and  Z  values  are  reversed  on  the 
input  data  cards,  and  the  j  input  is  made  negative. 

The  punched  output  containing  the  coordinates  is  in  a  format  compatible 
with  the  input  format  requirements  of  an  existing  program  which  will 
plot  these  points  for  fairing.  However,  the  coordinate  information  can 
readily  be  used  to  manually  plot  the  traces  of  the  plates. 
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Input  Format 


The  card  format  for  X.  ,  Y  ,  and  Z  ,  la  at  follow* : 


S  I  2  i 

«  5  5  5 

00  fH 


13  12  II  12  13  12  II  12  13  12  II  12 


A  complete  Input  data  deck  la  conatructed  in  the  following  manner: 

13  13 

The  first  card  Is  the  plate  Identification,  e.g.,  D-bb7,  B-blA, 
13 

D-101,  ...  .  The  second  card  contains  the  number  of  frames  of 
the  plate  (13  format),  and  the  plate  thickness  (F7.4  format). 
The  third  card  will  be  the  start  of  the  offsets,  as  previously 
described. 


The  offsets  for  the  upper  edge  of  the  plate  must  be  first  in 
order,  followed  by  the  middle  line  offsets,  and  finally  the 
offsets  for  the  lower  edge. 


Output  Format 

The  first  card  of  punched  output  is  the  plate  identification,  which 
is  identical  to  the  first  card  of  input  data.  The  remaining  cards 
are  the  X-Y  coordinates  of  the  plate.  The  coordinates  of  the  upper 
edge  on  every  frame  are  punched  first,  and  then  the  coorespondlng 
coordinates  for  the  lower  edge  are  punched.  The  computed  coordinates 
are  in  inches  at  one -tenth  scale.  The  cards  containing  the  upper 
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edge  points  have  a  format  as  follows : 


13  F9.3  13  F9.3 

XU(n)-+XXXX.XXXbYU(n)+XXXX.XXX 

XU  and  YU  being  X  upper  and  Y  upper 

The  cards  for  the  lower  edge  have  the  same  format,  except  XU  and  YU 
are  replaced  by  XL  and  YL  (X  lower  and  Y  lower). 


Operating  Instructions  for  IBM-1620 


1.  All  switches  to  off  and  program 

2.  Clear  memory  and  reset 

3.  Sense  Switch  3  ON,  if  printout  of  lengths  is  desired 

4.  Load  Object  Deck  and  push  "Load"  Switch  on  1622 

5.  Load  Data  Deck(s)  on  command,  "load  data"  message  printed 
by  typewriter,  and  push  "reader  start"  switch  on  1622 

6.  Load  Punch  Hopper  and  push  "punch  start"  on  1622 

7.  Data  processing  is  complete  when  (a)  all  data  cards  have  been 
read  properly,  (b)  output  has  been  punched,  and  (c)  a  "reader- 
no-feed"  light  on  1620  console 

8.  A  non-process  runout  of  the  1622  card  punch  must  be  made  to 
obtain  the  last-punched  output  card 


Timing 

1.  Sense  Switch  3  ON  -  Load  program,  load  data,  and  execute 
^»1.4  min/ frame 

2.  Sense  Switch  3  OFF  -  Load  program,  load  data,  and  execute 

1.04  min/ frame 
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Sample  Problem  -  The  Set  of  Input  Data  for  the  Shell  Plate  to  be  Developed 


Section  VIII 


SAMPLE  PROBLEM 


Figure  6-3  shows  a  set  of  input  data  for  a  shell  plate  to  be  developed. 

The  form  on  which  this  data  is  recorded  was  developed  for  use  by  the  manual 
loft.  Below  is  the  listing  of  the  coordinates  of  the  pattern  developed 
which  was  typed  as  the  program  was  executed. 


Output  (From  computer  typewriter) 

Listing  of  the  coordinates  of  the  pattern  developed 


A-  10A 


xu( 

1)» 

.0018YU 

xu( 

2)- 

2.4018YU 

xu( 

3)« 

4.8018YU 

xu( 

4)** 

7 .20 19 YU 

xu( 

5)- 

9.6019YU( 

xu( 

6)-  1 

I2.0019YU( 

XU( 

7)« 

14.4020YU 

xu( 

8)- 

16.8021YU 

xu( 

9  - 

19. 20 22 YU 

XL( 

1)~ 

.0004YL 

XL( 

2)« 

2.3995YL( 

XL 

3)- 

4.7995YL 

XL( 

4  - 

7.1 995 YL 

XL  ( 

5)- 

9.5995YL 

XL 

6)- 

1 1 .9995YL ( 

XL 

7)« 

14.3995YL 

XL( 

8)« 

16. 7995 YL 

XL  ( 

9)- 

19.1 995 YL 

1) -  5.2069  FR  116  A 

2) «  5.2006  FR  116  B 

3) -  5.2014  FR  116  C 

4) -  5.2011  FR  116  D 

5) -  5.2019  FR  116  E 

6)-  5.1958  FR  116  F 

7) -  5.1972  FR  116 

8) -  5.1924  FR  117 

9) -  5.1936  FR  118 

1)—  3.6794  FR  116  A 

2)—  3.6797  FR  116  B 

3) —  3.6799  FR  116  C 

4) —  3.6792  FR  116  D 

5) —  3.6795  FR  116  E 

6) —  3.6792  FR  116  F 

7) —  3.6793  FR  116 

8) —  3.6795  FR  117 

9) —  3.6852  FR  118 
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Section  IX 


READ  OFFSETTS  IN  FORM  OF  FEET 

ADD  1 

\  .Map\  t 

INCHES  .EIGHTS  A  SIXTEENTHS  k 

'/t  PLATE  l__ 

L  *  O 

CONVERT  TO  FEET  k  DECIMAL 

THICKNESS! 

J  >  1 

f  '  '  *  11  * 

FRACTIONS  THEREOF.  RESULTING 
NOTAT.ON  >VU.Z13 _ 

TO  Ytj  ' 
OFFSETt j 

NE  * 

,^.L  1 

sci  -(X1+a,j-xi,j)*  + 

Yv)*+  (zi+a,j  -  zx,j )* 

»C2.(X1  +  Ii,-  XIfJ)S 

yl*E,J-vl,3)k  ♦  (*Ul,J  -xx,j)* 

DOtl-  l,NE 
L  -L+l 


A I  -  SCI,  M  ■  |  k 

At  ■  let, AW  ■  »£> 
At  t  >61 


$CU)  •  A  4 


Rev  1, 

8  24  6?  0*29 

1*0. 0-2 


E 


COMPUTe  TRUE  LENGTH  BETWEEN  1*?  U 
FRAME*  ON  MIDDLE  LINE  OF  PLATE 


,xi+a,j  -  *i,j)  +  (Yi  +  a,j  -Yx,j)  ♦  (zi+  a, 3  -  zx,j) 

(XX+  1,3  “  *X,j)*  +  (Yl4l,J  -  Yl,j)a+  (ZX  ♦  l,J  -  z  I  ,j)* 
(Xx+a,7-  *i  +  i,j)a  +  (Yn.a,j- Yx+  i,j^4+  (zi+a,j  -Zi  +  i.j)* 


IN  -  2 
A  I  *  SB  I 
A2  a  SB2 
A3  a  SB3 
A&  -  S8  2 


DO  IS  I  .  I,  N2 
l.  L  +  l 


SB  I  a  AS 


COMPUTE  TRUE  LENGTHS  BETWEEN  FRAMES 
ON  LOWER  EDGE  OF  PLATE 


SAI  -(X1  +  2iJ-XI  j)  +(yi+2ij  -yi,j)  +(^I  +  2,J  “zl,j)* 
sa  a  a(Xi»2ij~xl4l|J)*  4  (Yi42,j  -  Yx  +  i,j)*+fzi*a.j  -  zi  +  i,j)r 


IN  >  S 
A  I  s  SAI 
A  2  a  SA2 
A3  .BAS 
At.  »  SAI 


DO  10  Xa  I,  N 
L  a  L+l 


SAU).AS 


COMPOTE  TRUE  LENGTHS  BETWEEN  UPPER  k  LOWER  EDGES 


5CA  a(Xi(j42 +(YI,3t2  “  YI,'3)t  ♦(zX,j'.*2-zI,j)1 

SC  0  a(Xx,j+|  -Xjtj)2  +(Yi)J4|  -  yI,j)*  +  — Zj  ,J )*" 

|&BA  ■  (*i,j,t.;>-.Xitj+1V+(Yi,3»g  -  YX.J»l)*  +(.ZI,J»2  -zl.3  +  if 


•«#  v- 


L- 

INt  A 

Al.  SCA 

) 

^4  t  SC  o 

A3  a  A4 

A3  -SBA 
A6»BC» 

Rev.  1. 
8-24-62 
1.0. 0-2 
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Rev.  1. 
8-24-62 
1.0. 0-2 


Rev .  1 . 
8-24-62 
1.0. 0-2 


XX  a  Z  A2  Al  +  2  AZ  AS  +  Z  A»  AU  At*-  AS-  A»* 
XX  «(XXA4>)/ j>>3  AZ  At) _ 


A3  *  I  +  xx/i4  +  3XX*/640+  5XX^/ 7168  T- . 

....+  35  XXA/t949\Z  +  63  XXS/2»e3Bft4 
A3  .  A2>Vm> 


COMPUTE 6 
OO  TO  IN  • 

2 ,13, 15, n, 16.21 


COMPUTED 
GO  TO  IP  . 


Section  X 


(  TOtlKAM  PROGRAM  LI ST IMG 

DIMENSION  X(  53,3) ,Y(  53,3), Z(  53,3),SC(  51),SA(  51),ACB(  53) 
DIMENSION  DSCB(  51 5 ,DSABC  51) 

DIMENSION  BX(  53),BY(  53),ABA(  53) ,SB(2) , AC (  53) 

DIMENSION  XL(53),XU(53),YU(53),YL(53) 

90  READ  108, ID 
READ  100, N,T 
DO  9  J-1,3 
DO  3  1-1, N 

1  REA0101 , IXF, IXI , I XE, IXS, IZF, IZI , IZE, IZS, IYF, IYI , I  YE, I YS 
XF-IXF 

XI-IXI 
XE-IXE 
XS-IXS 
YF-IYF 
YI-IYI 
YE- 1  YE 
YS-IYS 
ZF-IZF 
ZI-IZI 
ZE-IZE 
ZS-IZS 

X(l,J)-XF+(XI/12.)+(XE/96.)+(XS/192.) 
Y(l,J)«YF+(YI/12.)+(YE/96.)+  YS/192.)+(T/24.) 

(  3  Z( I ,J)-ZF+(ZI/12.)+(ZE/96.)+(ZS/192.) 

9  CONTINUE 
1-1 
J-1 
L-0 

N2-N-2 
DO  2  1-1, N2 
8  L-L+1 

SC1-(X(I+2,J)-X(I,J))*(X(I+2,J)-X(I,J)) 
SC1-SC1+(Y(I+2.J)-Y(I,J))*(Y(I+2.J)-Y(I ,J)) 

SCI-  (SC1+(Z(I+2,J)-Z(I,J))*(Z(I+2.J)-Z(I,J))) 

SC2-(X(I+1,J)-X(I ,J))*(X(I+1,J)-X(I,J)) 
SC2-SC2+(Y(I+1,J)-Y(I,J))*(Y(I+1,J)-Y(I,J)) 

SC 2-  SC2+(Z(I+1,J)-Z(I .J))*(Z(I+1,J)-Z(I,J)) 

SC3-(X(  l+2.J)-X(  1+1  ,J))*(XU+2,J)-X(  1+1  ,J)) 
SC3-SC3+(Y(l+2lJ)-Y(l+1,J))*(Y(l+2lJ)-Y(l+1,J)) 

SC3-  (SC3+(z( 1+2, J)-Z( 1+1 , J) )*(Z( 1+2, J)-Z( 1+1 , J) ) ) 

IN-1 

A1-SC1 

A2-SC2 

A3-SC3 

A6-SC2 

GO  TO  300 

2  SC(L)-A3 
L-1 


( 
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•1 


It  1-1 

12  SfM-(X(l+2.J)-X(l,J))*(X(l+2.J)-X(l,J)) 

SB1-SB1+(Y(  l+2,J)-Y(l,J))*(YU+2.J)-Y(l,J))  %%% 

SB1-  (SB1+(Z(I+2.J)-Z(I ,J))*(Z(I+2,J)-Z(I,J))) 

SB2-(X( 1+1 ,  J)-X( l,J))*(X(l+1,J)-X(l,J)j 
SB2-SB2+(Y(I+1.J)-Y(I,J))*(Y(I+1.J)-Y(I,J)) 

SB2-  (SB2+(Z(I+1,J)-Z(I,J))*(Z(I+1,J)-Z(I,J))) 

SB3-(X(I+2.J)-X(I+1.J))*(X(I+2.J)-X(I+1.J)) 

SB3-SB3+(Y( 1+2. J)-Y( 1+1 , J))*(Y( l+2.J)-Y( 1+1 ,  J)) 

SB3-  (SB3+(Z( l+2,J)-Z(l+1,J))*vZ(l+2,J)-Z(l+1,J))) 

A1-SB1 

A2-SB2 

A3-SB3 

IN-2 

A6-SB2 

GO  TO  300 

13  SB(L)-A3 
J-3 

1-1 

L-0 

0015  1-1, N2 

14  L-L+1 

SA1-(X(I+2.J)-X(I,J))*(X(I+2.J)-X(I,J))  % 

SA1-SA1+(Y{ l+2,J)-Y( I , J))*(Y{ 1+2. J)-Y( I ,  J) ) 

SA1-  (SA1+(Z(l+2fJ)-Z(l,J)j*(Z(l+2.J)-Z(l,J))) 

SA2-(X( 1+1 ,J)-X( ltJ))*(X(l+1.J)-X(l,J)) 

SA2-SA2+(Y( 1+1 . J)-Y( I , J))*(Y( 1+1 . J)-Y( I , J)) 

SA2-  (SA2+(Z(I+1,J)-Z(I,J))*(Z(I+1,J)-Z(I,J))) 

SA3-(X( 1+2, J)-X( l+1,J))*(X(l+2,J)-X(l+1,J)) 

SA3-SA3+(Y( 1+2. J)-Y( 1+1 , J) )*<Y( l+2.J)-Y( 1+1 I, J)) 

SA3-  (SA3+(Z( 1+2, J)-Z( 1+1 , J) )*(Z( 1+2, J)-Z( 1+1 , J) ) ) 

A1-SA1 

A2-SA2 

A3-SA3 

A6-SA2 

IN-3 

GO  TO  300 

15  SA(L)-A3 
1-1 

J-1 

L-0 

DO  18  t — 1 , N 
L-L+1 

SCA- ( X ( I . J+2)-X( I ,  J))*(X( I . J+2)-X( I , J)) 

SCA-SCA+(Y( I , J+2)-Y( I , J))*(Y( I . J+2)-Y( I , J)) 

SCA-  (SCA+(Z( I , J+2)-Z( I , J))*(Z( I , J+2)-Z( I ,J))) 


Rev.  1. 

HftSh 


( 


( 


< 


t7 


906 

905 

18 


19 


SCB-(X( I .J+1 )-X( I , J))*(X( I .J+1 )-X( I ,  J) ) 
SCB-SCB+(Y(I,J+1)-Y(I,J))*(Y(I . J+1 )-Y( I , J) ) 

SCB-  (SCB+(Z( I ,J+1)-Z( I ,J))*(Z( I ,J+1)-Z(I ,J))) 
SBA-(X( I . J+2)-X( I , J+1 ) )*(X( I . J+2)-X( I ,  J+1 ) ) 

SBA-SBA+(Y( I , J+2)-Y( I , J+1 ) )*(Y( I . J+2)-Y( I , J+1 ) ) 

SBA-  (SBA+(Z( I , J+2)-Z( I , J+1 ) }*(Z( I , J+2)-Z( I , J+1 ) ) ) 

A1-SCA 

A2-SCB 

A3-SBA 

A6-SCB 

IN-4 

GO  TO  300 
ACB(L)-A3 

XX-2 .  ■ *A2*A 1 +2 . *A2*A4+2 . *A4*A1 -A2*A2-A4*A4-A 1 *A 1 
XX-XX/ (A2*A1 ) 

I F(XX-1. >905,905,906 
TYPE  115 

A3-1.+XX/24.+3.*XX*XX/640.+5.*XX*XX*XX/7l68. 

A3-A3+35.*XX**4/294912.+63.*XX**5/2883584. 

A3-SQRTF ( A4)*A3 
ABA(L)-A3 
1-1 
J-1 
L-0 


DO  21  1-1 ,N2 
L-L+1 

CA-(Y( I ,J+2)-Y( I ,J))*(Y( I ,J+2)-Y( I ,J>) 

CA-CA+(Z( I , J+2)-Z( I , J))*(Z( I , J+2)-Z( I , J)) 

CB-(Y( l,J+1)-Y(l,J))*(Y(l ,J+1)-Y( I ,J>) 

CB— CB+ ( Z ( I , J+1 )-Z( I ,J))*(Z( I , J+1 )-Z( I ,J>) 

BA-(Y( I , J+2)-Y( I , J+1 ) )*(Y( I , J+2)-Y( I , J+1 )) 

BA-BA+(Z( I , J+2)-Z( l,J+1))*(Z(l , J+2)-Z( I , J+1 ) ) 

CB1-(Y(I+1 , J+1 )-Y(l,J))*(Y(l+1 .J+1 )~Y( I , J) ) 

CB1-CB1+(Z( l+1.J+1)-Z(I.J))*(Z(l+1 , J+1 )-Z( I , J) ) 

BA1-(Y( 1+1 . J+1 )-Y( I , J+2) )*(Y( I +1 -J+1 )— Y ( I , J+2) ; 
BA1-BA1+(Z( 1+1 ,J+1 )-Z( I , J+2))*(Z( 1+1 , J+1 )-Z( I , J+2)) 

A1-CA 

A2-CB 

A3-BA 

A6-CB1 

IN-5 

GOTO  300 

DSCB(L)-SQRTF(A3*A3+(X( 1+1 , J+1 )-X( I , J+1 ) )*(X( 1+1 , J+1 )-X( 


A1-CA 

A2-CB 

A3-BA 

A6-BA1 
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0-41 


I, J+1)) 


I 


) 


IN-6 

GO  TO  300 

21  DSAB(L)-SQRTF(A3*A3+(X( 1+1 ,J+1)-X( I ,J+1))*(X( 1+1 ,J+1 )-X( I ,J+1)) 

I F( SENSE  SWITCH  3)400,405 

400  N8-N-2 
L-0 

DO  401  1-1, N8 
L-L+1 

401  TYPE  105  ,L,SC( I ) ,L,SA( I ) 

L-0 

N1-N-1 

DO  402  1-1, N1 
L-L+1 

402  TYPE  106  ,L,ACB( I ) ,L,ABA( I ) 

L-0 

DO  403  1-1, N8 
L-L+1 

403  TYPE  107  ,L,DSCB( I ) ,L,DSAB( I ) 

405  1-1 

J-1 

L-1 

AI-SB(L) 

A2-ABA(L) 

A3-DSAB(L) 

IP-1 

A5-SB(L) 

GO  TO  301 

23  YL f L ) i ! S 2 ( *RT*S ) / SB (L ) L ^ } **2" (DSAB { L J ) **2 ) 7 ( 2 **SB ( L J *AB A ^ L ^ 

XL ( L )-ABA( L )*COSB 

PUNCH  108, ID 

TAU1-A5 

AI-SB(L) 

A2-ACB(L) 

A3-DSCB(L) 

IP-2 

A5-SB(L) 

GO  TO  301 

24  YU(L ) 1 \ 2B*R T*S ) /SB ( L ) ( L J ) **2_ ( DSCB ( L ) ^ **2 J 7 ( 2 »*SB ( L ^ *ACB ( L } } 
XU(L)-ACB(L)*COSB 

TAU2-A5 

L-2 

BX ( L )— SB ( L— 1 ) 

BY(L)-0.0 

L-1 

AI-DSAB(L) 


0-42 


R*v.  1. 


1 


( 


( 


( 


)) 


A2-SA(L) 

A3-ABA(L+1) 

IP-3 

A5-ABA(L+1) 

GO  TO  301 

25  GAM2-A5 

30  BET1«XL(L)/YL(L) 

BET1-ATANF(ABSF(BET1 ) ) 

36  BET2-XU(L)/YU(L) 

BET2-ATANF(ABSF(BET2) ) 

1 F ( XL ( L) ) 37 *35 ,35  ^ 

37  ALPH2-1 ,5708-(TAU1+GAM2+BET1 ) 

GO  TO  34 

35  ALPH2-1 ,5708-(TAU1-BET1+GAM2) 

34  YL(L+1)-YL(L)+(SA(L)*SINF  ALPH2) 

XL(L+1 )«XL(L)+(SA(L)*COSF( ALPH2) ) 

L-1 

AI-DSCB(L) 

A2-SC(L) 

A3-ACB(L+1) 

IP-4 

A5-ACB(L+1) 

GO  TO  301 

38  GAM3-A5 
IF(XU(L))26, 28,28 

26  ALPH3-TAU2+GAM3+BET2-1.5708 
GO  TO  29 

28  ALPH3-TAU2-BET2+GAM3- 1.5708 

29  YU(L+1)-YU(L)+(SC(L)*SIHF(ALPH3) 

XU(L+1 )-XU(L)+(SC(L)*C0SF(ALPH3) ) 

L-2 

N1-N-1 
D047  l ml  N1 

39  AC(L)»((XU(L)-XL(L) )*( XU(L)-XL(L) ) )+( (YU(L)-YL(L) )*(YU(L)-YL(l) 


AC(L)»SQRTF(ABSF(AC(L) )) 

I F(XL(L)-XU(L) )41 ,40,42 

40  XLAH-1 .5708 
GO  TO  43 

41  XLAM«ATANF( (YU(L)-YL(L) )/( XU (L )— XL ( L) ) ) 

GO  TO  43 

42  XLAM-3J416-ATANF((YU(L)-YL(L))/(XL(L)-XU(L))) 

43  AI-AC(L) 

A2-DSAB  L) 

A3-DSCB(L) 

IP-5 

A5-DSCB(L) 

GO  TO  301 
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0-43 


44 


45 


46 


47 


51 

54 


58 

52 

55 


RH0-A5 

XNU-XLAM-RHO 

BY(L+1)-YL(L)+(DSAB(L)*SINF(XNU)) 

BX(L+l)«XL(L)+(DSAB(L)*COSF(XNU)) 

S1-S 

RT1-RT 

AI-DSAB(L) 

A2-ABA(L+1 ) 

A3-SA(L) 

IP-6 

A5-ABA(L+1 ) 

GO  TO  301 
DELT-A5 
SIGH-XNU-OELT 

YL(L+1 )«YL(L j+(SA(L)*SI NF(S IGH 
l)-XL(L)+  ' 


xl(l+i )«xl(l)+(sa(l)*cosf(s igh) ) 

EPS-ATANF( ( (2.*RT1 )/(S1-DSAB(L) ) 
XNEW-3 • 1 41 6-XLAM-EPS 
Al-DSCB(L) 

A2-SC(L) 

A3-ACB(L+1) 

IP-7 

A5-ACB(L+1) 

GO  TO  301 
DELT-A5 

SIGH-DELT-XNEW 

YU(L+1 )-YU(L)+(SC(L)*SI NF(S IGH) ) 
XU(L+1)-XU(L)+(SC(L)*C0SF(SIGH)) 
L-L+1 

DO  48  L-1 ,N1 
XU(L)«XU(L)*1.2 
YU(L).YU(L)*1.2 
IF(XU(L))51,52,52 
XU(L)-ABSF( XU(L) ) 

IF  YU(L))54,58,58 
YU(L)-ABSF(YU(L) ) 

PUNCH  110,l,XU(L),L,YU(L) 

GO  TO  48 

PUNCH  113,L,XU(L),L,YU(L) 

GO  TO  48 

IF(YU(L))55,56,56 
YU(L)-ABSF( YU(L) ) 

PUNCH  114,L,XU(L),L,YU(L) 

GO  TO  48 

PUNCH  103,L,XU(L),L,YU(L) 
CONTINUE 
DO  49  L-1 ,N1 
XL(L)-XL(L)*1 .2 


)/( 


1  ,-(RT1/($1-DSAB(L) ) )**2) ) 


1 


YL(L)«YL(L)*1.2 
IF(XL(L))61 ,62,62 

61  XL(L)-ABSF(XL<L)) 

IF(YL(L))64,68,68 

64  YL(L)-ABSF(YL(L)) 

PUNCH  109,L,XL(L),L,YL(L) 

GO  TO  49 

68  PUNCH  111,L,XL(L),L,YL(L) 

GO  TO  49 

62  I F  C YL ( L ) ) 65 , 66 , 66 

65  YL(L)-ABSF(YL(L) )  ,  % 

PUNCH  112  ,L,XL(L) ,L,YL(L) 

GO  TO  49 

66  PUNCH  102,L,XL(L),L,YL(L) 

49  CONTINUE 

GO  TO  90 

300  A4-A3 

XX-2 .*A2*A1 +2 ,*A2*A3+2 ,*A3*A1-A2*A2-A3*A3-A1 *A1 
XX-(XX*A6)/(A3*A2*A1) 

IF(XX-1.)900,900,901 
901  TYPE  115 

900  A3-1 ,+XX/24.+3«*XX*XX/640.+5.*XX*XX*XX/7l68. 
A3-A3+35.*XX**4/294912.+63.*XX**5/2883584. 
A3-SQRTF(A6)*A3 
GO  TO  (2, 13.15, 17, 19, 21), IN 

301  S-(A1+A2+A3)/2.  %  %% 

RT-SQRTF(ABSF( ( (S-A1 )*(S-A2)*(S-A3) )/S) ) 

A5-ATANF ((2.*RT/( S-AS ) ) / ( 1 .- ( RT/ ( S-A5 ) ) **2 ) ) 

GO  T0(23, 24, 25, 38, 44,45, 46), IP 

101  FORMAT( 1 3, 12,11,12,13,12,11,12,13,12,11,12) 

100  FORMAT (13, F7.4)  ,  .  . 

102  FORMAT (4H  XL( , I 3,2H)«,F10.4,3HYL( , I3,2H)-,F10.4) 

103  FORMAT (4H  XU( , I3,2H)-,F10.4,3HYU( , I 3,2H)-,F10,4) 

105  F0RMAT(3HSC(.I3,2H)-,F9.4,4X,3HSA(, I3,2H)-.F9*4) 

106  FORMAT (4HACB(. 1 3,2H)-,F9.4,4X,4HABA( , 1 3 ,2H5-.F9*4) 

107  FORMAT(5HOSCB(.l3,2H)-,F9.4f4X,5HDSAB(, I3,2H)-,F9.4) 

108  FORMAT (2HA-,  3)  , 

109  FORMAT (4H  XL  ,  1 3,3H)— , F9.4,3HYL( ,  1 3 ,3H)— ,F9.4) 

110  F0RMAT(4H  XU(,I3,3H)— ,F9.4,3HYU(,l3,3Hj— ,F9.4j 

111  FORMAT (4H  XL(,I3,3H)— ,F9.4,3HYL(, I3,2H)-,F10.4) 

113  FORMAT (4H  XU( , 13, 3H)— ,F9.4,3HYU(, I3,2H)-,F10,4) 

112  FORMAT (4H  XL( ,  1 3,2H j-,F10.4,3HYL( ,  1 3 ,3H j— ,F9*4) 

114  FORMAT (4H  XU( , I3,2H)-,F10.4,3HYU( , 13, 3H)— ,F9.4) 

115  F0RMAT(20HS/R  GREATER  THAN  ONE) 

END 
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Appendix  H 


CALCULATION  OF  OFFSETS  FOE  FRAME  BENDING  TEMPLATES 


Bending  templetes  for  accurately  determining  the  shape  of  frames  for 
forming  are  one  form  of  Information  required  of  any  lofting  system. 
Assuming  an  equation  Is  at  hand  for  the  shape  of  the  hull  at  a  given 
frame,  this  program  calculates  the  proper  offsets  for  describing  the 
frame  bending  template. 

The  program  Is  written  in  FORTRAN  II  for  the  IBM- 1620  computer.  The 
method  of  solution,  operating  details,  a  sample  problem,  flow  diagram, 
and  program  listing  are  Included  in  this  Appendix. 


METHOD  OF  SOLUTION 

The  input  information  necessary  to  solve  the  problem  consists  of: 

— ■  The  function  that  describes  the  shape  of  the  hull,  Y  *  F  (Z) 

The  end  points  of  the  frame  segments,  Z^  and  Z ^ 

—  The  step  size  between  solution  points,  S 

The  offset  (V)  which  gives  depth  to  the  template,  or  allows 
for  the  width  of  the  frame. 


1.0. 0-2 


H-l 


This  problem  la  graphically  presented  in  Tig.  H-l.  The  solution  is  as  follows: 


1 


GIYEH:  Y  -  f(Z),  Z^  Z2>  V  ,  S 

TXHB:  OFF  SKIS  WL  at  intervals  S  as  shown  in  Fig.  H-l 


METHOD: 

1. 


■  f(Z)  at  Z1 
-  f(Z)  at  Z2 


2. 


M  - 


Z1  -  Z2 
Y1  -  Y2 


3. 


Y 


tan"1  |m| 


4.  H  -  S  sin  Y 
R  *  S  cos  # 
F  -  V  sin  P 
G  -  V  cos  P 

5.  $>  -  90°  -  X 

6.  YAj  -  Yj  -  G 


’■  l-W- 

8.  ZKX  -  Zjl+^F 

Z  EOT  «  Z2  +$T 


&  assumes  a  value  of  4-1  if  the 
slope  is  positive,  -1  if  the 
slope  is  negative 


and  have  now  been 

established.  The  offset  at  these  points  lsy  of  course 
equal  to  V  . 


9.  The  Coordinates  of  Points  PA1 
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A-3 


10.  Step  through  the  following  for  eech  required  offset  , 
L  -  1  ,  2  ,  3  ,  ....  until  Z^  +  H 


11-  ^  •  H 
aL  “  aL-l  ‘  H 

-  SR 

"t  ‘WL-1  -  <f* 

12.  The  solution  for  “l  involves  «  trial  end  error  solution. 
Newton's  Method  ves  selected  for  this: 


First : 

YC^  “  f(Z)  evaluated  at  ZC^ 

(1) 

YC^  -  AZc£  +  BZcJ  +  CZC^  +  D 

Also: 

(2) 

ycl  "  "l  "  M(s>  <“l  ■  “l> 

Subtracting 

(2)  from  (1)  gives 

(3)  0  -  AZCjJ  +  BZc£  +  CZC^  +  »  -  Y^  -  M^XM^-  Z^), 


ZC^  is  the  only  unknown  in  (3)  and  may  be 
solved  for  using  Newton's  Method. 


13.  YCl 
!4.  ^ 


f(z)za 


<ycl  ’  ^  >2  + 


(ZCl  -  z^)2 


H-4 
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A  definition  of  the  program  input  symbols  end  card  formats  are  given 
below: 


Input  Symbols 


A 

B 

C 

D 

IDSXA 


Z2 

DIS 


S 


Coefficients  of  a  standard  cubic  equation  describing 
the  molded  shape  of  the  frame.  The  equation  is  of  the 
form  AX3  +  BX2  +  O  +  D 

Four  digit  integer  frame  identification 

Upper  limit  of  the  frame  segment  to  be  analyzed 

Lower  limit  of  the  frame  segment  to  be  analyzed 

Minimum  offset  to  allow  for  depth  of  frame  (  F  in  Fig.  H-l) 

Distance  along  chord  of  curve  between  required  offsets. 


First  Card  Format: 

Variable 


foPE,  ISSSS& 


z. 


Format 


2F10.5 


Card  Columns 
15 


B 

4F15.8  j 

‘  16  -  30 

C 

(31-45 

D 

1 

^46  -  60 

BIS 

S 

2F7.3 

(61  -  67 
\68  -  74 

IDSIA 

15 

75  -  79 

f 1-10 

0  -  20 
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The  order  of  output  Is  as  follows: 

1.  The  frame  Identification  is  typed 

2.  T)w  coordinates  of  and  (Wg-  H-l)  are  typed 

3.  The  coordinates  of  each  successive  point  ?A^  and  the  offset 
at  that  point  are  typed 

gtfSpyt  Messages 

If  the  slope  of  the  chord  of  the  frame  segment  equals  zero,  a 
message  is  typed  and  the  angle  y  is  set  to  90° 


«•  St&iisu, 

Switch  1  -  OFF,  After  completion  of  the  calculations  for  a  given 
problem,  the  program  looks  for  an  entire  new  set 
of  data  in  the  card  reader. 

Switch  1  -  OH,  After  completion  of  a  problem,  the  program  looks 

for  a  new  limits  card  (Z^  ,  Z^)  in  the  card  reader. 

All  other  switches  are  Ignored 


1.0. 0-2 
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f  SAMP  LX  PKQBLXM  -  FUffi  BIND  IDG  PROGRAM 

Oat*  for  Sample  Problem  -  frame  No.  1234 

F(Z)  »  - .00013834Z3- .00210378Z2- .330263Z+36.0 

OIS  -  3.0  S  -  3.0 

Zj  -  14.0  Z2  -  0.0 

1620  Typewriter  Listing  while  executing  sample  problem: 


260000200003^)  •*  • —  . . . —Memory  Clearing  Instruction 

LOAO  SUBROUTINES 
ENTER  DATA 


FRAME  NO.  1 234 

Z(1)  Y ( 1 )  Z(2)  Y(2) 

14.00000  30.58437  0.00000  36.00000 


PT.  NO. 
1 
2 


3 

4 

5 

6 
7 


Z  Y 

12.91815  27.78623 

10.12000  28.86807 

7.32185  29.94991 

4.52371  31.03175 

1.72556  32.11359 

-1.07258  33.19543 

-1.08185  33.20185 


PERP.  01  ST.  TO  CURVE 
3.00000 
3.16161 
3.22862 
3.21436 
3.13377 
3.00320 
3.00000 


FLO*  CHART  - 

3JLAXX0K  COT  OFFSETS  FOB.  FUM  XSXBTMB  HMFUTSS 


■2 


H-8 
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RftOUM  LIST  ISC 


fwrT  |  i 

iNG-^7^53982+  • 995354* (PM/ PL  j-  . 2886  79* ( ( PM/PL )**3 ) 


OFFSETS  FOR  FRAME  BENDING 
DIMENSION  Z(3) »Y(3) 

Y(lj-O.0 
Y(2)-0,0 
AL  PA-0.0 

t  READ  100.A.B.C,D,DIS,S,  IDSTA 
PRINT  !05,IDSTA 

14  READ  101,Z(1).Z(2) 

CALL  CURVE  AfS,C,D,Y(n,Z(n) 
CALL  CURVE(A,B,C,D,Y(2),Z(2)) 
PRINT  106 

PRJNT  102.Z(1),Y(1),Z(2).Y(2) 
PRINT  104 

4  SLP-(Z(1 )-Z(2))/(Y(1 )— Y(2) ) 
PSLP-— SLP 
SP-ABSF(SLP) 

PH-SP-1 . 

PL-SP+1 . 

IF(SLP)26f21,26 

26  ANG-.  7853982+. 9 _ , _ - 

ANG-ANG+.079331*((PM/PL)**5 
23  H-S*SINF(ANG) 

R-S*COSF(ANG) 

BANG-(1 .5708- ANG) 
YA-Y(lj-DIS*COSF(BANG) 

YB-Y ( 1 ) 

Q-0 1 S*S I NF ( B  ANG ) 
IF(SLP)15,16,16 

15  ZA-Z(1 )-Q 
ZEND-Z(2)— Q 
GO  TO  17 
16  ZA-Z(1)+Q 
ZEND-Z(2)+Q 

17  ZB-Z(1) 

L-1 

,  PRINT  103,L,ZA,YAfDIS 
6  L-L+1 
ZA-ZA-H 
ZB-ZB-H 

I F(SLP) 25  *18,18 
25  YA-YA+R 
YB-YB+R 
GO  TO  19 

18  YA-YA-R 


t 


YB-YB-R 
19  ZC-Z8 

10  CALL  CURVE(A,B.C,D,ALPA,ZC) 
ALPA-ALPA-PSLP*( ZC-ZB  )-YB 

7  BETE-3.*A*ZC*ZC+2,*B*ZC+C-PSLP 
ZC1*ZC— (ALPA/BETE) 

I  F(ABSF(ZC1-ZC)-,0001 )5,5»1 1 

11  ZC-ZC1 
GO  TO  10 

5  ZC-ZC1 

CALL  CURVE(A,B,C,D,ALPA,ZC) 

8  DZ-(ZC-ZA)*(ZC-ZA) 
DY-(ALPA-YA)*(ALPA-YA) 
DIF-SQRTF(DZ+0Y) 

PRINT  103,L,ZA,YA,DIF 

I F ( ZA-ABSF ( ZEND+ABSF (H) ) ) 9 , 6 , 6 

9  L-L+1 

YA-Y( 2)-D I S*COSF(BANG) 

PRINT  107.L.ZEND.YA.DIS 
22  IF(SENSE  SWITCH  1)14,1 
21  PRINT' 108 
ANG-1.5708 
GO  TO  23 
100  FORMAT ( 


101  FORMAT 

102  FORMAT 

103  FORMAT 

104  FORMAT 

105  FORMAT 

106  FORMAT 

107  FORMAT 

108  FORMAT 
END 


4F15.8f2F7.3,l5) 

2F10.5) 

4F10.5/) 

I4,10X,F10.5,2X.F10.5,13X,F10.5) 

7HPT,  N0,,10X,1HZ,12X,1HY,10X,20HPERP,  DIST.  TO  CURVE) 
9HFRAME  NO., 15/) 

3X,4HZ(1).7X,4HY(1),7X.4HZ(2),7X.4HY(2)) 

I4,10X,F10.5,2X,F10.5,13X,F10.5/J 

20H  SLOPE  EQUALS  ZERO) 


SUBROUTINE  CURVE(A.B,C,D,FX,X) 

FX-A*X*X*X+B*X*X+C*X+0 

RETURN 

END 


l 
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Appendix  I 


DECK  OFFSETS  ROUTINE 


INTRODUCTION 

The  purpose  of  this  program  is  to  provide  vertical  deck  camber  offsets 
describing  some  standard  deck  forms  easily  expressed  mathematically. 

The  configurations  which  this  program  handles  are: 

—  Straight  line  sheer  and  straight  line  camber 

—  Straight  line  sheer  and  parabolic  camber 

—  Parabolic  sheer  and  parabolic  camber 

The  program  is  intended  to  be  as  general  as  possible  within  the  confines 
of  the  above  cases.  It  is  possible  to  vary  the  frame  spacing  along  the 
length  of  the  ship  and  to  vary  the  transverse  interval  between  offsets. 
It  is  also  possible  to  have  variations  in  slope  of  the  straight  sheer 
line  and  to  put  knuckles  in  the  straight  line  camber. 


FORMULAS 

The  offsets  calculated  by  this  program  are  vertical  offsets,  whose  base 
plane  is  a  horizontal  plane  at  a  height  above  the  base  line  of  the  ship, 
equal  to  the  height  of  the  lowest  point  on  the  aheer  curve.  The  formulas 
for  finding  the  height  of  the  deck  at  the  centerline  are  given  below: 

(1)  The  height  of  the  deck  at  the  centerplane  amidships 

h  -  .04  Y  -  2y  (.24)/ 12 
n  n  'n 

where  yn  is  the  molded  halfbreadth  amldship 

(2)  The  height  of  the  deck  at  the  centerplane  and  at  the 
forward  perpendicular 

hf  -  (.2Pl  +  20)/ 12 

where  F^  is  the  length  between  perpendiculars 
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(3)  The  height  at  the  aft  perpendicular 

h#  -  (.1PL+  10)/ 12 

Figure  1  gives  the  configuration  of  the  decks  in  the  longitudinal  direction. 
Figure  2  presents  the  cross  sections  of  the  decks.  The  definitions  of  the 
symbols  used  are  the  same  as  given  in  the  FORTRAN  Input  Symbol  Definitions. 

The  program  is  written  in  Fortran  II  for  the  IBM- 1620  computer. 


INPUT  DATA 


Input  Symbol  Definition 


Symbol 


Definition 


D 


PL 

P 

YM 

DIST 


CST 

LI 

L2 

L3 

L4 


Distance  from  the  midship  section  of  the  ship  to 
the  section  that  has  the  lowest  freeboard  (lowest 
point  on  the  sheer  profile).  D  Is  negative  if 
this  point  is  aft  of  midship. 

Ship's  length  between  perpendiculars 

Overall  length  of  the  ship 

The  maximum  molded  half-breadth  of  the  ship 

Difference  in  length  between  the  overall  length 
of  the  ship  measured  forward  from  midship  and 
the  length  between  perpendiculars  measured  for¬ 
ward  from  midship. 

Length  of  the  transverse  interval  between  two 
deck  offsets  (Fig.  2 -a). 

For  the  cases  having  straight  sheer,  this  is  the 
number  of  different  slopes  in  the  sheer  line 
forward  of  midship.  (LI  ■  3  in  Fig.  1-b). 

Same  as  LI,  except  aft  of  midship  (L2  -  2  in  Fig.  1-b). 

For  the  cases  having  straight  camber,  this  is  the 
number  of  different  slopes  in  the  camber  from  the 
center  plane  to  the  sheer  line  (L3  *  2  in  Fig.  2-b) . 

The  number  of  intervals  of  different  frame  spacing 
forward  of  the  section  with  minimum  freeboard.  The 
minimum  number  of  intervals  must  be  2.  That  is, 
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there  muet  be  at  least  one  change  In  frame  spacing. 

L5  -  The  number  of  Intervals  of  different  frame  spacing 
used  aft  of  the  section  of  minimum  freeboard.  The 
maximum  number  is  2  . 

DB  -  Length  of  each  different  slope  (corresponding 

to  LI)  from  midships  forward.  Each  DF  is 
presented  in  order  starting  at  midship  and  working 
forward . 

BA  -  Same  as  DF  (corresponds  to  12  )  starting  at 

midship  and  working  aft. 

8  -  Length  in  feet  of  each  longitudinal  section  having 

a  different  station  spacing  (corresponds  to  U  ) 
starting  at  the  point  of  minimum  freeboard  at 
working  forward. 

V  -  The  interval  between  frames  in  each  section  E  above. 

BA  -  Same  as  B  ,  except  aft  (corresponds  to  L5  ) 

VA  -  Same  as  V  ,  except  aft 

T  -  Length  in  feet  of  each  of  the  slopes  (corresponds 

to  L3  )  from  the  centerline  plane  to  the  sheer  line. 


Input  Data  Cards 

There  are  four  different  sets  of  data  specified  for  this  program, 
one  set  for  each  of  the  four  cases  given  in  the  Introduction. 

The  limitations  on  those  variables  which  haye  dimension  restrictions 
are  given  below. 


Variable  Minimum  Maximum 


E  ,  EA  2  20 

V  ,  VA  ?.  20 

DB  ,  DA  1  21 

T  1  21 


The  actual  FORTRAN  format  field  descriptions  have  been  used  to 
describe  the  data  cards.  The  field  descriptions  are  the  FORTRAN 
F  field,  which  contains  a  fixed  point  decimal  number,  and  the 
I  field,  which  conteins  an  integer  number  that  is  always  right 
justified. 
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A.  Parabolic  Sheer  and  Cambsr 


Card  1  (This  card  la  a  haader  card.  For  this  case  It 
contains  a  1  In  Column  5) 


Card  2 


Format 

F15.9 

F15.9 

F15.9 

F15.9 

F15.9 

Columns 

1-15 

16-30 

31-45 

46-60 

61-75 

Variable 

D 

PL 

P 

YK 

HIST 

Card  3 

Format 

F15.9 

15 

15 

Columns 

1-15 

16-20 

21-25 

Variable 

CST 

U 

L5 

Next  14  Cards 

Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

S 

V 

Next  15  Cards 

Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

IA 

VA 

B.  Straight  Sheer  and  Camber 

Card  1  (Header  card  contains  a  2  In  Column  5) 


Card  2 


Format 

F15.9 

F15.9 

F15.9 

F15.9 

F15.9 

Columns 

1-15 

16-30 

31-45 

46-60 

61-75 

Variable 

H 

PL 

P 

YM 

HIST 

Card  3 

Format 

F15.9 

15 

15 

15 

15 

15 

Columns 

1-15 

16-20 

21-25 

26-30 

31-35 

36-40 

Variable 

CST 

LI 

12 

L3 

U 

L5 

Next  LI  Cards 

Format  F15.9 
Columns  1-13 
Variable  OB 

Next  L2  Cards 

Format  F15.9 
Columns  1-15 
Variable  HA 

Next  L3 

Format  F15.9 
Columns  1-15 
Variable  T 
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Next  14  Cards 


Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

E 

V 

Next  15  Cards 

, 

Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

EA 

VA 

Straight  Sheer  > 

and  Farabollc  Camber 

Card  1  (Header  card ; 

contains  a  3 

in  Column  5) 

Card  2 

Format 

F15.9 

F15.9  F15.9 

F15.9 

F15.9 

Columns 

1-15 

16-30  31-45 

46-60 

61-75 

Variable 

0 

PL  P 

YM 

DIST 

Card  3 

Format 

F15.9 

15  15 

15 

15 

Columns 

1-15 

16-20  21-25 

26-30 

31-35 

Variable 

CST 

11  L2 

14 

15 

Next  LI  Cards 

Format 

F15.9 

Columns 

1-15 

Variable 

DB 

Next  12  Cards 

Format 

F15.9 

Columns 

1-15 

Variable 

DA 

Next  14  Cards 

Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

E 

V 

Next  L5  Cards 

Format 

F15.9 

F15.9 

Columns 

1-15 

16-30 

Variable 

EA 

VA 

D.  Parabolic  Sheer  and  Straight  Camber 


Card  1  (Header  card; 

contains 

a  4  in  Column  5) 

Card  2 

Format 

F15.9 

F15.9 

F15.9 

F15.9 

F15.9 

Columns 

1-15 

16-30 

31-45 

46-60 

61-75 

Variable 

D 

PL 

P 

YM 

DIST 

Card  3 

Format 

F15.9 

15 

15 

15 

Columns 

1-15 

16-20 

21-25 

26-30 

Variable 

CST 

L3 

L4 

L5 

Next  L4  Cards 

Format 

Columns 

Variable 

F15.9 

1-15 

E 

F15.9 

16-30 

V 

Next  15  Cards 

Format 

Columns 

Variable 

F15.9 

1-15 

EA 

F15.9 

16-30 

VA 

OUTPUT 

All  the  output  will  be  punched  on  cards.  The  deck  height  offsets  of 
Section  D  will  be  punched  first,  followed  by  those  of  the  sections  forward 
of  D  .  After  completing  these  sections  the  message  "Deck  Offsets  Aft" 
will  be  punched.  The  message  will  be  followed  by  the  offsets  from  the 
section  just  aft  of  D  to  the  stern. 

For  each  section  the  program  will  punch  offsets  out  to  the  halfbreadth  (YM). 
At  stations  forward  and  aft  of  the  midship  station  (or  parallel  middlebody) 
some  of  these  offsets  will  not  be  valid  since  they  will  lie  outside  the 
sheer  line. 

SENSE  SWITCH  SETTINGS 

All  the  sense  switches  are  normally  off,  with  one  exception.  The  value  of 
the  transverse  interval  between  offsets  (CST)  can  be  changed  during 
execution  by  turning  Switch  3  on.  The  program  will  halt  and  wait  for  a 
new  value  of  CST  to  be  typed  in  (Format  F15.9).  Switch  3  must  be 
turned  off  before  pushing  start  to  continue. 
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( 
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SAMFIX  na BUM 
Input  Pat* 

1 


-10.0 

50.0 

4.0 

2 

20.0 

10.0 

20.0 

10.0 

10.0 

10.0 

10.0 

5.0 

Output 

60.0 


12.0 


5.0 


DECK  OFFSETS  FORWARD. - 
FRAME  -10.0000 
.47999999 
.42666666 
.26666667 
.00000001 
FRAME  0.0000 

.64489794 
.59156461 
.43156462 
.16489796 
FRAME  10.0000 

1.13959170 
1 .08625840 
.92625847 
.65959181 
FRAME  20.0000 

1.96408150 
1.91074810 
1.75074810 
1.48408150 
FRAME  30.0000 

3.11836710 
3.06503380 
2.90503380 
2.63836720 
DECK  OFFSETS  AFT.- 
FRAME  -20.0000 

.82222223 
.76888893 
.60888893 
.34222233 
FRAME  -25 .0000 

1.24999990 
1.19666660 
1 .03666660 
.77000003 
FRAME  -30.0000 

1 .84888870 
1.79555540 

1 ifiiiilao 


1  60  OVERALL  1 

1— _ SO*  LBP _ J 

s' 

»o‘ 

25' 

_ 

EiCTtQW  ^ 
^  , 

a 

plane  a 

5  / 

7 

1 

1  1 

SIDE  VIEW  or  SAMPUE  PROBLEM 
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FLOW  DIAGKAM 


1.0. 0-2 


1-11 


1-12 


compute  offsets  for 

PARABOLIC  SHEER  {  CAM  BET? 
CONDITION. 


1.0. 0-2 


1-13 


1.0. 0-2 


1-15 


i 


DETERMINE  IF 
THIS  STA. 
COMPLETE 


6 


1 


2 


1-18 


OFF 


3 


3 


0 
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& 


85 


I 


C.sZj-AN^X-OaJ 

. II 


Y=0 


I 


(funchA 

FFAMtt  ) 

■£ 


ETsC,-AFY* 

I 


/punch\ 

(^uuV) - H>) 


1-21 


PUNCH 
l  VALUB 
^  IT  . 


I 


1.0. 0-2 


i-:3 


CpunchA 

FRr6J 


Cj*C,  +  A*TI* 


1.0. 0-2 


1-25 


< 


< 


( 


1.0. 0-2 


I- 


CpuwchN 

PRrJ 


00  73i|. 


Ct  =  C,  -  A*^1* 

am-  rc1.,-cx)/rri-71.j) 


73JT 


1-29 


•Wi 


2 


1-30 


ooooo 


FORTRAN  PROGRAM  LISTING 


*  *  *  DECK  OFFSETS  *  *  * 

IDNT-1  PARABOLIC  SHEER  AND  CAMBER 

IDNT- 2  STRAIGHT  SHEER  AND  CAMBER 

IDNT-3  STRAIGHT  SHEER  AND  PARABOLIC  CAMBER 

IDNT-4  PARABOLIC  SHEER  AND  STRAIGHT  CAMBER 
DIMENSION  DB(21),DA(21),E(20),EA(20),V(20),VA(20),T(21),C(21) 
DIMENSION  Z ( 2 1 ) 

LE-1 

GO  TO  600 
26  PAUSE  2 
600  READ  506, I DNT 

READ  500,0, PL, P,YM,DIST 
HF-(.2*PL+20.)/12. 

HA" ( , 1 *PL+1 0 , ) / 1 2 , 

HM-(2.*YM*.24)/12. 

PL2-PL/2. 

PL2SQ-PL2*PL2 

DSQ-D*D 

A«HM/(YM*YM) 

A1«(HF-HM)/(PL2SQ-2.*D*PL2+DSQ) 

A2-(HA-HM)/(PL2SQ+2.*D*PL2+DSQ) 

B1—  2.*D*A1 

B2— 2.*D*A2 

AB1-B1*B1/(4.*A1)+HM 

AB2-B2*B2/(4.*A2)+HM 

P1-PL2+DIST 

P2-P-P1 

GO  TO  (131,132.130,133), IDNT 
131  READ  502,CST,L4,L5 
34  DO  52  1-1, L4 

52  READ  500 , E ( I ) , V ( I ) 

DO  53  1-1, L5 

53  READ  500, EA( I ) , VA( I ) 

E( 1 )-E( 1 )+D 

EA( 1 )-EA( 1 )— D 
DO  41  1-2, L4 

41  E( I )-E( I )+£( 1-1 ) 

DO  42  1-2, L5 

42  EA( I )»EA( I )+EA( 1-1 ) 

LE-1 

X-D 

GO  TO  (43,44,701,702), IDNT 
43  PUNCH  503 

54  Y-0. 

PUNCH  505, X 

56  ZP— A*Y*Y+A1*X*X+B1*X+AB1 
PUNCH  500, ZP 
IF(SENSE  SWITCH  3)  30,31 
30  ACCEPT  500, CST 


55 

n 

59 


61 

63 


32 

33 


31  Y-Y+CST 

IF  YM-Y)55.56,56 
IF  E(LE)-X  58,57.59 
IF(U-LE)54,60,5& 

LE-LE+1 
X-X+V(LE) 

IF(P1-X)60,54,54 
60  PUNCH  504 
LE-1 

X-D-VA(I) 

Y-O. 

PUNCH  505.X 

ZP— A*Y*Y+A2*X*X+B  2*X+AB  2 
PUNCH  500, ZP 
I FCSENSE  SWITCH  3)  32,33 
ACCEPT  500, CST 
Y-Y+CST 

IF(YM-Y)62,63.63 
62  IF(-EA(LE)-X)66.64,65 

64  IF(L5-LE)61,26,65 

65  LE-LE+1 

66  X-X-VA(LE) 

I F(-P2-X)6l ,61 ,26 
132  READ  502,C1,L1,L2,L3,L4,L5 
L3-L3+1 
L1-L1+1 
L2-L2+1 
DO  1  1-2, LI 
1  READ  500, DB( I ) 

DO  3  1-2, L2 
3  READ  500, DA( I ) 

DO  5  1-2, L3 
5  READ  500 , T ( I ) 

GO  TO  34 
44  DA( 1 )«D 
DB ( 1 ) -D 
Z( 1 )-HM 
T(  1 )-0. 

PUNCH  503 
DO  27  1-2, L3 
27  T( I )-T( I )+T( 1-1 ) 

DO  6  J-2.L1 
X-DB(J-I) 

DB( J)-DB( J)+DB( J-1 ) 
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Z( J)-A1*DB(J)*DB(J)+B1*DB( J)+AB1 
AN«(Z(  J)-Z(  J-1 )  )/(r>B(J)-DB(J-1 )  ) 
C(1 )-Z( J-1 )+AN*(X-DB( J-1 ) ) 
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35 

36 


Y-0. 

PUNCH  505.X 
DO  7  1-2. L3 
C(l)-C(1)-A*T(n*T(l) 

AM-(C( l-1)-C(l))/(T(l)-T(l-1)) 
ZT-C( 1-1 )-AM*(Y-T( I  —  1 ) ) 

PUNCH  500, ZT 

IF(SENSE  SWITCH  3)  35,36 

ACCEPT  500, Cl 

Y-Y+C1 


I F(L3-I ) 1 1 , 10, 1 1 

10  IF(T(I )-Y)7,9,9 

11  IF(T( I )-Y)7,7,9 
7  CONTINUE 

IF(E(LE)-X) 12,28,13 
28  I F(L4-LE)6,6, 12 

12  LE-LE+1 

13  X-X+V(LE) 

I F(L1-J) 15, 14,15 

14  IF(DB(J)-X)6,8,8 

15  I F(DB(J)-X)6,6,8 
6  CONTINUE 

PUNCH  504 
LE-1 

X-D-VA(I) 

DO  16  J-2.L2 
DA(J)—  DA(J)+DA(J-1) 
Z(J)-A2*DA(J)*DA(J)+B2*DA(J)+AB2 
AN«(Z( J)-Z( J-1 ) ) / (DA( J— 1 )-DA( J) ) 
19  C(1)-Z(J)-AN*(X-0A(J)) 

Y-0. 

PUNCH  505.X 

DO  17  1-2, L3 

C( I )-C( 1 )-A*T( I )*T( I ) 

AM-(C( I —1 )— C ( I ))/(T{ I )— T( 1—1 )) 


18 

ZT-C( 1-1 )-AM*(Y-T( 1 

-i: 

PUNCH  500, ZT 

IF(SENSE  SWITCH  3) 

37 

37 

ACCEPT  500, Cl 

38 

Y-Y+C1 

IF(L3~I )21 ,20,21 

20 

IF(t(I)-Y)17,18,18 

21 

IF(T(I)-Y)17,17,18 

17 

CONTINUE 

1 F(-EA(LE)-X) 25,29, 

24 

29 

IF(L5-LE)16, 16,24 

24 

LE-LE+1 

25 

22 

23 

16 

130 


70 

72 

701 


75 


X-X-VA(LE) 

I F(L2-J)23.22,23 
IF  DA(J)-X)19J9,16 
IF(DA(J)-X)19J6,16 
CONTINUE 
GO  TO  26 

READ  502,CST,L1,L2,L4,L5 

L1-L1+1 

L2-L2+1 

DO  70  1-2, LI 

READ  500 , DB ( I ) 

DO  72  1-2, L2 
READ  500, DA( I ) 

GO  TO  34 

2(1 )»HM 

DB  ( 1 )-D 

D A ( 1 )-0 

PUNCH  503 

DO  74  1-2, LI 

DB ( I )-DB( I )+DB( 1-1 ) 

Z( I )-A1*DB( I )*DB( I )+B1*DB( I )+AB1 
AN-(Z(I)-Z(I-1))/(DB(I)-DB(I-1)) 
C(1)-Z(I-1)+AN*(X-DB(I-1)) 

Y-0. 


PUNCH  505 ,X 
76  ZT-C(1)-A*Y*Y 

PUNCH  500, ZT 
I FCSENSE  SWITCH  3)766,67 
766  ACCEPT  500, CST 
67  Y-Y+CST 

I F(YM-Y) 77,76,76 
77  IF(E(LE)-X)79,78,80 

78  IF(L4-LE)74,74,79 

79  LE-LE+1 

80  X-X+V(LE) 

81  1 F ( P 1  — X) 74.82,83 

82  I F  (DB  ( 1  )-X)74,75,75 

83  I F(DB( 1 5 — X) 74, 74,75 
74  CONTINUE 

PUNCH  504 
LE-1 

X-D-VA(LE) 

DO  84  1-2, L2 

DA(  I )— DA(J  )+DA(  1-1 ) 

Z( I )-A2*DA( I )*DA( I )+B2*DA( I )+AB2 
AN-(Z( I )-Z( l-1))/(DA( l-1)-DA( I )) 
85  C(1)-Z(I)-AN*(X-DA(I)) 

Y-0. 
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PUNCH  505.X 
86  ZT-C(1)-A*Y*Y 

PUNCH  500, ZT 

I F ( SENSE  SWITCH  3)  762,763 

762  ACCEPT  500, CST 

763  Y-Y+CST 
IF(YM-Y)87,86,86 

87  IF(-EA(LE)-X)90.88,89 

88  IF(L5-LE)84,84,89 

89  LE-LE+1 

90  X-X-VA(LE) 

91  I F(-P2-X)93,92.84 

92  I F(DA( I )-X)85,o5,84 

93  I F (DA( I )— X) 85 ,84, 84 
84  CONTINUE 

GO  TO  26 

133  READ  502,CST,L3,L4,L5 

T ( 1 )-0. 

L3-L3+1 
GO  TO  34 

702  DO  732  I-2.L3 

732  READ  500 , T ( I ) 

DO  150  l-2,L3 

150  T( I )-T( I )+T( 1-1 ) 

PUNCH  503 

733  C(1)-A1*X*X+B1*X+AB1 

Y-0. 


PUNCH  505.X 

DO  734  I-2.L3 

C( I )-C( 1 )-A*T( I )*T( I ) 

AM-(C( l-1)-C(l))/(T(l)-T(l-1)) 
735  ZP«C( 1-1 )-AM*(Y-T( 1-1 ) ) 

PUNCH  500, ZP 
I F(SENSE  SWITCH  3)68,69 

68  ACCEPT  500, CST 

69  Y-Y+CST 


I F  ( L3- 1 )  73  7 , 736 , 73  7 
IF(T( I )-Y) 734, 735,735 
I F( T( I )— Y) 734, 734, 735 

734  CONTINUE 

I F(E(LE)-X) 738-39,40 

39  I F ( L4-LE ) 733 , 742 , 738 
738  LE-LE+1 

40  X-X+V(LE) 

I F ( P 1 — X) 742,733,733 
742  PUNCH  504 
LE-1 

X-D-VA(I) 


736 

737 
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743  C(1)-A2*X*X+B2*X+AB2 

Y-0. 

PUNCH  505, X 
DO  744  1-2, 13 
C(I)-C(1)-A*T(I)*T(I) 
AM-(C(I-1)-C(I))/(T(I)-T(I-1)) 
45  7P-C( 1-1 )-AM*(Y-T( I  — 1 ) ) 


PUNCH  500, ZP 

I F{ SENSE  SWITCH  3)  764,765 

764  ACCEPT  500, CST 

765  Y-Y+CST 

IF  1.3-1)47,46,47 

46  I F(T( I )-Y >744,45,45 

47  I F(T( I )-Y) 744, 744,45 
744  CONTINUE 

IK-EA(LE)-X)48,49,50 

49  IF(L5-LE)743,741,50 

50  LE-LE+1 

48  X-X-VA(LE) 

51  I F(-P2-X) 743, 743,741 
741  GO  TO  26 

500  FORMAT(5F15.9) 

502  F0RMAT(F15.9,5 15) 

503  FORMAT  (22HDECK  OFFSETS  FORWARD.-) 

504  FORMAT  (18HDECK  OFFSETS  AFT.-) 

505  FORMAT (6H FRAME  .F10.4) 

506  FORMAT  (215) 

END 
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PROGRAM  FOR  INNER  EDGE  OF  WEB  FRAMES 


The  purpose  of  this  program  is  to  find  the  inner  contour  of  a  web 
frame  when  given  offsets  describing  the  mold  line  contour  of  the 
frame.  That  is,  given  a  set  of  points  which  lie  on  a  smooth  mold  line 
of  a  frame,  we  wish  to  apply  a  transformation  to  these  points.  This 
transformation  is  to  be  such  that  the  transformed  points  lie  on  a  curve 
which  is  similar  in  contour  but  which  now  passes  an  ordinate  distance 
from  the  upper  point  and  C2  from  the  lowest  point.  (See  Fig.  J-la.) 

To  accomplish  this  transformation,  the  curve  is  first  translated  the 
distance  •  By  examining  Fig.  J-la  it  can  be  seen  that  if  the  curve 
is  now  rotated  so  that  the  end  point  was  a  distance  C2  from  its  original 
position,  the  resulting  curve  would  be  too  short  and  in  a  distorted 
position.  To  overcome  this  condition  the  curve  is  extended  by  fitting  a 
parabola  to  the  three  lowest  data  points  and  extrapolating  a  new  end 
point. 

A  special  condition  which  may  sometimes  be  found  on  web  frames  is  the 
existence  of  a  knuckle  near  the  upper  point.  This  condition  is  illustrated 
in  Fig.  J-lb.  The  inner  contour  will  be  a  smooth  curve  from  the  beginning 
at  the  upper  edge.  The  true  distance  to  which  the  points  below  the  knuckle 
must  be  translated  is  not  given  by  .  Therefore,  a  new  must  be 
calculated  which  gives  a  closer  approximation  of  the  correct  distance. 

The  new  distance  is  approximated  by  passing  a  parabola  through  the 

uppermost  three  data  points  below  the  knuckle,  or  if  given  including  the 
point  at  the  knuckle  itself,  can  then  be  found  by  solving  for  the 
point  on  the  parabola  corresponding  to  the  upper  frame  edge. 

The  program  is  written  in  FORTRAN  for  the  IBM-1620  computer. 
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a.  Input  Bata 


FORTRAN  Input  Symbols 


N  -  Number  of  points,  must  be  an  integer  with  no 
decimal  point,  and  must  be  right  Justified  to 
Column  5 .  N  <  25 

Cl  m  Ordinate  distance  curve  is  to  be  moved  at  first 
point 

C2  -  Ordinate  distance  curve  is  to  be  moved  at 
last  point 

CC  ■  Abscissa  distance  from  first  point  to  place 
where  Cl  is  given.  This  applies  only  when 
the  curve  contains  a  knuckle 


X(I),Y(I) 


Coordinates  of  points  describing  the  curve 


Data  Cards 


There  are  a  total  of  (N+l)  cards 


First  Card: 


Variable 

Format 

Card  Columns 

N 

15 

1-5 

Cl 

Cb  -  15 

C2 

3F10.5 

<16  -  25 

CC 

(26  -  35 

If  Sense  Switch  1  is  OFF,  CC  will  be  ignored  and 
nothing  need  be  punched  in  Columns  26  -  35. 

Remaining  (N)  cards : 

Variable  Format  Card  Columns 


jl-15 
(16  -  30 
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b.  Sanaa  Switch  Settings 

Switch  1  ON  -  Frame  Curve  contains  a  knuckle 
Switch  1  OFF  -  Frame  curve  does  not  contain  a  knuckle 

c.  Output  (Via  Typewriter) 

Transformed  points  (X,Y) 
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Stiff!*  IMBUM 


7  “2.0 

-3.5 

8.0 

10.296 

7.0 

9.758 

6.0 

9.157 

5.0 

8.516 

4.0 

7.853 

3.0 

7.065 

2.0 

5.540 

260000  200003R5 
LOAD  DATA 


Output 

8.000000 

8.296000 

7.046353 

7.679557 

6.097758 

7.000316 

5.152369 

6.281204 

4.208745 

5.540163 

3.275142 

4.674525 

2.400629 

3.074258 

2.000000 

2.040000 

LAYOUT  OF  SAMPLE  PROBLEM. 


TUN  DIAGRAM 


FIRST  POINT 
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1 


I 


PftOGRAM  URHC 


*  * 


1 

50 

51 

52 

20 

40 

1 

i, 

41 

1 

1 

1 

T 

1  1 

42 

21 


30 

31 

2 

3 

4 

5 


CONTOUR  TRANSFORMATION  *  * 

DIMENSION  X( 25) , Y ( 25 ) 

I F ( SENSE  SWITCH  1)50,51 
RE AD 100, N, Cl ,C2,CC 
GO  TO  52 

READ  100, N, Cl ,C2 
DO  20  1-1  .N 
READ101 ,X( I ) ,Y( I ) 

I F(SENSE  SWITCH  1)40,42 

BB-(X(3)-X(2))*(Y(4)-Y(2))-(X(4)-X(2))*(Y(3)-Y(2)) 

BB-BB/((X(3)“X(2))*(X(4)-X(2))*(X(4)-X(3))) 

AA«(Y(3)-Y(2))/(X(3)-X(2))-BB*(X(3)-X(2)) 

Z-X(1)+CC 

WW-Y(2)+AA*(Z-X{2))+BB*(Z-X(2))**2 
C1-C1-WW+Y ( 1 )+CC*(Y( 1 )-Y(2))/(X(1)-X(2)) 

DO  41  1-1, N 

K-N+1-1 

X(K+1)-X(K) 

Y(K+1 )-Y(K) 

X(1)-X(1)+CC 
Y  ( 1 )  -WW 
N-N+1 
TX«X(  1 
TY-Y ( 1 ) 

SXX-X(N) 

SYY-Y(N)+C2 
DO  21  -1 ,N 

x( I )-X( I j-TX 
Y(I)-Y(I)-TY 

B-(Y(N-2)-Y(N) )*( X(N-1 )-X(N))-(Y(N-1)-Y(N))*( X(N-2) 
B-B/( ( X(N-2)-X(N) )*{ X(N-1 )-X(N) )*( X(N— 2)-X(N-1 ) ) ) 
A-(Y(N-1 )-Y(N) )/( X(N-1 )-X(N) )-B*( X(N-1 )-X(N) ) 
D-X(N)**2+(Y(N)+C2-C1 )**2 
Z-X(N) 

T-Z**2+Y ( N ) **2-D 

IF(T)2,9.31 

Z-.5*Z 

GO  TO  30 

I F(Z)4,9,3 

H-.5 

GO  TO  5 

H— .5 

CA-T 

Z-Z+H 

W«Y(N)+A*(Z-X(N))+B*(Z-X(N))**2 
T«Z**2+W**2-D 
I F(CA*T)6,9,5 
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-X(N)) 


I 


I 


9 


6  HA-ABS(H) 

I F ( HA- . 0001)7, 8,8 

7  Z-Z-H/2. 

GO  TO  9 

8  H— H/2. 

GO  TO  5 

9  C3«(Z-X(N))**2+(Y(N)+C2-C1-W)**2 
CS-1.0-.5*C3/D 

SN-SQRTO  .0-CS**2) 

DO  25  1-1,  N 
IF(C2*X(N})22,22,23 

22  SX-X( I )*CS+Y( I )*SN+TX 
SY— X(  I  )*SN+Y(  I  )*CS+TY+C1 
GO  TO  25 

23  SX«X(l)*CS-Y(l)*SN+TX 
SY«  X( I )*SN+Y( I )*CS+TY+C1 

25  PRINT  101 ,SX,SY 

PRINT  101.SXX.SYY 
GO  TO  1 

100  FORMAT( I5.3F10.5) 

101  FORMAT( 2F15 .6) 

END 
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ROUTINE  FOR  IDENTITY  MATRICES 


The  function  of  this  program  is  to  provide  an  identity  matrix  and 
associated  cards  for  providing  a  basic  feasible  solution  for  the 
matrices  produced  by  SMOG  1  and  SMOG  2. 

OPERATING  INSTRUCTIONS 


The  program  is  written  in  FORTRAN  II  for  the  IBM-1620. 

Fortran  Input  Variable  Definitions 

IOIM  -  The  number  of  the  first  column  in  the  "identity" 
matrix. 

IOWS  -  The  number  of  rows  in  the  non-basis  matrix  produced 
by  SMOG. 

COST  -  The  cost  on  the  columns  of  the  matrix 

Input  Format 

The  input  consists  of  one  card  with  the  following  format: 

Format  15  15  F10.4 

Variable  IOIM  IOWS  COST 

Output 

The  output  is  in  the  following  order: 

(1)  a  card  punched  "BASIS" 

(2)  A  set  of  basis  headings 

(3)  The  identity  matrix  including  the  cost  row  and  elements. 

Sense  Switches 
None  used 
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PROGRAM  LISTING 


1  READ  150, I OLM, IOWS, COST 
PUNCH  151 

I  OF- 1 OLM 
DO  2  1-1, IOWS 
PUNCH  1 52, 1  OF, I 

2  IOF-IOF+1 
I  OF- 1 OLM 

DO  3  1-1, IOWS 
PUNCH  153, IOF, COST 
PUNCH  1 54, 1  OF, I 

3  IOF-IOF+1 
GO  TO  1 

150  FORMAT ( 15, I  5 ,  F 1 0 . 4) 

151  F0RMAT(5HBAS IS) 

152  F0RMAT(7X,1HC, I4,2H  R,I4) 

153  F0RMAT(6X,2H  C,I4,6H  OFSET,F12.5) 

154  F0RMAT(7X,1HC, I4.2H  R,l4,8H  1.  ) 

END 


sample  input  data 
100  30  1000.') 
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SAMPLE  OUTPUT 


BASIS 


c 

100 

R 

1 

c 

101 

R 

2 

c 

102 

R 

3 

c 

103 

R 

4 

c 

104 

R 

5 

c 

105 

R 

6 

c 

106 

R 

“» 

/ 

c 

107 

R 

8 

r 

w 

108 

ft 

9 

c 

109 

R 

10 

c 

110 

R 

11 

c 

111 

R 

12 

c 

112 

R 

13 

c 

113 

R 

14 

c 

114 

R 

15 

c 

115 

K 

16 

c 

1 16 

R 

17 

c 

117 

R 

18 

c 

118 

R 

19 

c 

119 

k 

20 

c 

120 

R 

21 

c 

121 

R 

22 

c 

122 

R 

23 

c 

123 

R 

24 

c 

124 

R 

25 

c 

125 

R 

26 

c 

126 

R 

27 

r* 

U 

127 

R 

28 

r 

128 

R 

29 

c 

129 

R 

30 

c 

100 

OFSET 

1000.00000 

c 

100 

R 

1 

1  .0 

c 

101 

OFSET 

1000.00000 

c 

101 

R 

2 

1.0 

c 

102 

OFSET 

looo.onooo 

c 

102 

R 

3 

1.0 

c 

103 

OFSET 

1000.00000 

c 

103 

R 

4 

1 .0 

c 

104 

OFSET 

1000.  ocr.n.: 

n 

u 

104 

R 

5 

1 .0 

c 

105 

OFSET 

1000.00000 

c 

105 

R 

6 

1.0 

c 

106 

OFSET 

1000.00000 

c 

106 

R 

7 

1.0 
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c 

107 

OFSET 

c 

107 

R  8 

c 

108 

OFSET 

c 

108 

R  9 

c 

109 

OFSET 

c 

109 

R  10 

c 

110 

OFSET 

c 

110 

R  11 

c 

111 

OFSET 

c 

111 

R  12 

c 

112 

OFSET 

c 

112 

R  13 

c 

113 

OFSET 

c 

113 

R  14 

c 

114 

OFSET 

c 

114 

R  15 

c 

115 

OFSET 

c 

115 

R  16 

c 

116 

OFSET 

c 

116 

R  17 

c 

117 

OFSET 

c 

117 

R  18 

c 

118 

OFSET 

c 

118 

R  19 

c 

119 

OFSET 

c 

119 

R  20 

c 

120 

OFSET 

c 

120 

R  21 

c 

121 

OFSET 

c 

121 

R  22 

c 

122 

OFSET 

c 

122 

R  23 

c 

123 

OFSET 

c 

123 

R  24 

c 

124 

OFSET 

c 

124 

R  25 

c 

125 

OFSET 

c 

125 

R  26 

c 

126 

OFSET 

c 

126 

R  27 

c 

127 

OFSET 

c 

127 

R  28 

c 

128 

OFSET 

c 

128 

R  29 

c 

129 

OFSET 

c 

129 

R  30 

1000.00000 

1  .0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1  .0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1  .0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 

1000.00000 

1.0 
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